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Abstract 

A ubiquitous observation in cell biology is that the diffusive motion of macromolecules and organelles is anomalous, and a de- 
scription simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is 
commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarising their densely 
packed and heterogeneous structures. The most familiar phenomenon is a sublinear, power-law increase of the mean-square 
displacement as function of the lag time, but there are other manifestations like strongly reduced and time-dependent diffusion 
coefficients, persistent correlations in time, non-gaussian distributions of spatial displacements, heterogeneous diffusion, and 
a fraction of immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we 
summarise some widely used theoretical models: gaussian models like fractional Brownian motion and Langevin equations for 
visco-elastic media, the continuous-time random walk (CTRW) model, and the Lorentz model describing obstructed transport 
in a heterogeneous environment. Particular emphasis is put on the spatio-temporal properties of the transport in terms of two- 
point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even if the 
mean-square displacements are identical. Then, we review the theory underlying commonly applied experimental techniques in 
the presence of anomalous transport like single-particle tracking, fluorescence correlation spectroscopy (FCS), and fluorescence 
recovery after photobleaching (FRAP). We report on the large body of recent experimental evidence for anomalous transport in 
crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments 
where a variety of model systems mimic physiological crowding conditions. Finally, computer simulations are discussed which 
play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed 
by a synthesis of the theoretical and experimental progress identifying open questions for future investigation. 

PACS numbers: 05.40.-a 87.16.-b 
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1. Preface 

Transport of mesoscopic particles suspended in simple sol- 
vents is governed by Brownian motion and is one of the pillars 
of biological and soft condensed matter physics. The pioneer- 
ing works of Einstein and von Smoluchowski identified the 
random displacements as the essence of single-particle motion 
and — employing ideas of the central-limit theorem — suggest a 
gaussian probability distribution with a mean-square displace- 
ment that increases linearly in time. This paradigm constitutes 
the basis of numerous applications ranging from colloidal sus- 
pensions, emulsions, and simple polymeric solutions to many 
biological systems. 

From this perspective, the emergence of anomalous trans- 
port in more complex materials, typically characterised by a 
sublinear increase of the mean-square displacement, appears 
as exotic. Yet experiments from many different areas reveal 
that anomalous transport is ubiquitous in nature, signalling 
that slow transport may be generic for complex heterogeneous 
materials; examples are crowded biological media, polymeric 
networks, porous materials, or size-disparate mixtures. One 
prominent ingredient common to many of these systems is that 
they are densely packed, known as macromolecular crowding 
in the biophysics community. The presence of differently sized 
proteins, lipids, and sugars in the cell cytoplasm, as well as the 
filamentous networks permeating the cell, is believed to be at 
the origin of the observed slowing down of transport, render- 
ing diffusion constants meaningless that were measured in di- 
lute solution. Similarly in heterogeneous materials displaying 
pores of various sizes, the agents interact strongly with sur- 
faces of complex shape giving rise to an entire hierarchy of 
time scales, which implies a power-law increase of the mean- 
square displacement. 

The goal of this review is to provide a pedagogical 
overview of the existing theoretical framework on anomalous 
transport and to discuss distinguished experiments from the 
different fields, thus establishing a common language for the 
observations coming from various experimental techniques. 
The focus here is on crowded biological systems in vivo and 
in vitro both in the bulk systems such as the cytoplasm or bio- 
logical model solutions as well as in planar systems, for exam- 
ple, crowded cellular membranes. We also want to emphasise 
the similarities with complex heterogeneous materials as they 
naturally occur, e.g., porous rock, to technologically relevant 
materials such as molecular sieves and catalysts. 

In Section 2, we provide some elementary background to 
quantify the stochastic motion of a single particle in a complex 
disordered environment. We recall the standard arguments that 
underlie the theory of Brownian motion both on the level of 
an ensemble average (Fickian diffusion) as well as in terms 
of single-particle trajectories. In particular, we rederive the 



widely used mean-square displacement of a tracer and more 
generally the entire two-time conditional probability density. 
Thereby, we introduce some notation that will be employed 
throughout the review and exemplify the important concept 
of scaling. Then we contrast complex transport from normal 
transport by requiring that deviations from normal transport 
are persistent way beyond the natural time scales of the system. 
The most celebrated indicator is the subdiffasive increase of 
the mean-square displacement. A second example which has 
drawn significant interest in the statistical physics community, 
although almost entirely ignored by biophysicists, is given by 
the velocity autocorrelation function. We recall that even the 
case of a single particle in a simple solvent encodes non-trivial 
correlations and long-time anomalies due to the slow diffusion 
of transverse momentum. Further, we emphasise that anoma- 
lous transport can no longer be described by a single diffusion 
constant, rather the transport properties depend on the con- 
sidered length and time scales. Thus, it is essential to make 
testable predictions for the spatio-temporal behaviour, which 
we elucidate in terms of the nongaussian parameter and, more 
generally, the shape of the distribution of displacements. 

The theoretical description of the phenomenon of anoma- 
lous transport is reviewed in Section 3. We focus on the three 
most widely used concepts each of which has been corrobo- 
rated from some experiments and excluded from others. The 
first one consists of relaxing the white-noise assumption in a 
Langevin approach. Then to induce a sublinear increase of the 
mean-square displacement a power-law time correlation has 
to be imposed on the increments while the probability distri- 
bution is still assumed to be gaussian. In this simple model, 
one easily derives two-time correlation functions such as the 
van Hove correlation function or its spatial Fourier transform, 
the intermediate scattering function. Simple generalisations to 
models where the correlation functions encode only a power- 
law tail at long times then lead to the gaussian model. A sec- 
ond, complementary approach is given in terms of a distribu- 
tion of waiting times characteristic of the complex medium. 
The class of these models is known as the continuous-time 
random walk model and allows for a simple solution after a 
Fourier-Laplace transform to the complex frequency domain. 
While these two descriptions are phenomenological in nature, 
the third concept is a microscopic approach that also provides 
insight into the mechanism how subdiffusion can emerge from 
simple dynamic rules in a complex environment. In this class, 
which we refer to as Lorentz models, the origin of slow trans- 
port is found in an underlying geometric percolation transition 
of the void space to which the tracer is confined. We recall 
how self-similar, fractal structures emerge and why transport 
in these systems is naturally subdiffusive. We compare differ- 
ent microscopic dynamics and discuss the crossover to hetero- 
geneous diffusion at long time scales, scaling behaviour, and 
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the appearance of an immobile fraction of particles. 

The section is concluded by delimiting the scope of the 
review to subdiffusion that is related to molecular crowding or 
the excluded volume effect. Polymer physics and single-file 
diffusion, where subdiffusive motion is well known too, will 
not be covered here. In addition, while many experimental 
data are interpreted in terms of subdiffusion, some of them 
actually display only deviations from normal transport that are 
not persistent over large time windows. Hence one has to be 
careful whether one speaks of subdiffusion or simple crossover 
scenarios. Typical examples where such apparent subdiffusion 
is observed comprise diffusion of more than one species each 
moving with its own diffusion coefficient. In the same spirit a 
particle could change its conformation thus exhibiting several 
internal states again characterised by several diffusive regimes. 
Sometimes, the data only suggest power laws which could be 
interpreted as subdiffusion, yet the correlations are due to the 
measurement technique and have no relation to the underlying 
physical processes. 

Sections 4 and 5 review anomalous diffusion in crowded 
biological systems from an experimental point of view. We 
first introduce the experimental techniques that have proved 
themselves as useful tools in biophysics to measure transport 
properties in mesoscopic samples and on microscopic scales. 
An important technique that is both intuitive and powerful is 
provided by single-particle tracking. Here the trajectory of 
some agent in a complex environment is recorded over suf- 
ficiently long time, allowing for the evaluation of, in princi- 
ple, all correlation functions. The most widely used quantity 
is of course the mean-square displacement, since it is rather 
robust with respect to statistical fluctuations and appears easy 
to interpret. A second powerful method that has been widely 
applied in the biophysical context is fluorescence correlation 
spectroscopy. Here the basic idea is to label few molecules by 
a fluorescent dye and record the fluctuating fluorescent light 
upon illuminating a small part of the sample with a laser. We 
briefly introduce the theory underlying the measurement and 
discuss how anomalous transport manifests itself in the corre- 
sponding correlation function. Complementary to these two 
single-particle methods is fluorescence recovery after photo- 
bleaching, which detects the diffusion front of fluorophores 
after depleting a small spot by an intense laser pulse. The 
technique is apt to measure very slow transport and immobile 
particles, and we outline how anomalous transport becomes 
manifest in the recovery curve as function of time. 

Section 5 is devoted to the plethora of experiments on 
anomalous transport in the cell interior and related model sys- 
tems and their interpretation in terms of the three theoretical 
frameworks introduced before. We have compiled and dis- 
cuss the most significant experiments in the field, focusing on 
the past decade. Almost all experiments agree that transport is 
hindered and slowed down by molecular crowding, manifested 
in a suppression of the diffusion constant. A large subset of 



experiments provides clear evidence for subdiffusive motion, 
but there are notable experiments which report normal diffu- 
sion, and some findings appear to be specific to certain bio- 
logical conditions. Important progress in the understanding 
of anomalous transport was made by biologically motivated 
model systems, where in contrast to living cells key parame- 
ters are adjustable. The discussion is completed by addressing 
computer simulations of simple and complex model systems, 
which provide essential support for the interpretation of ex- 
periments. The presentation is divided into three-dimensional 
transport in cellular fluids, e.g., in the cytoplasm of living cells, 
and transport in cellular membranes, which may be approx- 
imated as a (curved) two-dimensional manifold. Finally, the 
consequences of anomalous transport on reaction kinetics are 
briefly sketched and some recent progress on this emerging 
topic is reported. 

In Section 6, we first summarise what the state of the art in 
the field of anomalous transport currently is and provide some 
general conclusions on what has been achieved so far both 
theoretically and experimentally. We point out which ques- 
tions are still under debate and give suggestions where future 
research may go and what crucial issues still need to be ad- 
dressed. 

2. Basics of Brownian motion 

The naming of the observation of anomalous transport fre- 
quently found in complex systems already suggests that the 
phenomena are fundamentally different from the standard case 
that therefore qualifies as normal transport. Before reviewing 
how anomalous transport can be addressed theoretically and 
measured experimentally, we discuss the framework for nor- 
mal transport connecting macroscopic diffusion with fluctua- 
tions at small scales. 

The erratic movement of a mesosized particle suspended 
in a simple solvent is referred to as Brownian motion, after 
the Scottish botanist Robert Brown, who observed the contin- 
uously agitated motion of minute particles ejected from cer- 
tain pollen grains under a light microscope. The first theo- 
retical description has been achieved independently by Ein- 
stein [1] and Smoluchowski [2] in terms of a probabilistic 
approach. These ideas have been rephrased shortly after by 
Langevin [3, 4] in terms of stochastic differential equations by 
separating the force balance into a deterministic and a random 
part. The characterisation of the random forces is largely due 
to Ornstein [5] thereby laying the foundations of the modern 
calculus of Langevin equations. 

The experimental demonstration of the probabilistic route 
to a macroscopic law has been achieved by Perrin [6] and his 
students around the same time by meticulously analysing sin- 
gle trajectories of colloidal particles observed under a micro- 
scope. His contribution was awarded with the Nobel prize 
in 1926 as a breakthrough in proving the physical reahty of 
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molecules. 

The impact of Einstein's theory on Brownian motion, i.e., 
the normal case, can hardly be underestimated and constitutes 
one of the milestones in physics. Recently, on the occasion of 
100* anniversary of Einstein's annus mirabilis 1905, a series 
of reviews have been published highlighting the new concepts 
and future directions in the field of Brownian motion [7-10]. 

2.1. Simple diffusion 

In the molecular kinetic approach advocated by Einstein and 
Smoluchowski the suspended particle experiences rapid colli- 
sions with the solvent molecules. These events occur at the 
time scale of the liquid dynamics, typically in the picosecond 
regime, and at each encounter a tiny amount of momentum is 
exchanged. Simultaneously, the collisions are responsible for 
the macroscopic friction force counteracting the random kicks. 
At the time scales where the suspended particle moves signifi- 
cantly, the increments AR(t) = R{t + t')-R(t') after an elapsed 
time t are considered as random variables that are identically 
and independently distributed. By the central-limit theorem, 
the total displacement, being the sum of many independent 
tiny increments, then is governed by a gaussian distribution 



P(r,t) = [27T6r^(t)/dY'"^exp[-r^d/2Sr^(t)) 



(1) 



where d is the dimension of the embedding space. The prob- 
ability distribution P(r, t) is known in general as the propaga- 
tor or the van Hove self-correlation function [11, 12]. For in- 
dependent increments the variance 5r^{t) :- ^[R(f) - R(0)]^^ 
grows linearly with the number of steps, implying a linear in- 
crease of the mean-square displacement, 6r^(t) - 2dDt. The 
only transport coefficient characterising the Brownian motion 
is then the diffusion constant D, completely specifying the 
propagator. 



(2) 



The self-similarity of free Brownian motion becomes evident 
by writing the propagator in a scale-free form. 



P(r,f) = r-"'!Pgauss(r), rocrf 



-1/2 



(3) 



introducing a dimensionless scaling variable, f :- (2Dt) '^^r, 
r - |r|, and a scaling function. 



nauss(r) = {2n)-''IH''txp{-f''l2). 



(4) 



Other scaling forms will be encountered in the course of this 
review. 

The connection to the macroscopic description arises when 
considering many particles performing Brownian motion such 
that the probability cloud P{r, t) displays the same space-time 
dynamics as the macroscopic concentration. In particular, one 



verifies that the gaussian propagator fulfils the diffusion equa- 
tion 



d,P{r, t) = DVV(r, f), 



(5) 



and satisfies the initial condition of a spatially localised distri- 
bution, P(r, f) = 6{r). 

For applications, a representation in terms of spatial 
Fourier modes is advantageous, 

and f (k, t) is known as the (self-)intermediate scattering func- 
tion [11, 12]. It can be measured directly by neutron scattering 
employing the spin-echo technique [13] or, on larger length 
scales, by photon correlation spectroscopy [14]. The momen- 
tum transferred from the sample to the photon or neutron is 
then simply /zk. For the diffusion propagator, one readily cal- 
culates 

P(k, f) = exp(-D/t2f) , (6) 

implying that density modulations decay with a rate 1/tj. = 
Dk^. Then, long-wavelength perturbations are long-lived since 
by particle conservation the relaxation involves transport of 
particles over large distances. 

For future reference, we also provide the van Hove corre- 
lation function in the complex frequency domain. 



P(r, u) 



f 

Jo 



e'"' P(r, t) dt. 



Im[w] > 0. 



(7) 



The one-sided Fourier transform reduces to the standard 
Laplace transform, provided one identifies s = -ito; the ad- 
vantage of introducing complex frequencies as above is that 
the transform can be readily inverted numerically by evaluat- 
ing the integral 



2 r°° 

P(r,f)=- Re[f(r,w)]cos(wf)dw. 
^ Jo 

For the diffusive propagator one finds 



P(r, oj) 



1 



( V-iwD 



i27rDyn 



d/2-l 



-lOJ 

~D~ 



(8) 



(9) 



where Ky(-) denotes the modified Bessel function of the second 
kind. The expression simplifies for the dimensions of interest. 



P(r, oj) 
Pir, w) 
P(r, = 



1 



2 V-iwD 



exp 



D 



2nD 
1 

4nDr 



exp 



D 



D 



(d = 1), (10a) 
{d = 2), (10b) 
(d = 3). (10c) 



In scattering techniques, where the energy transfer to the 
sample is also recorded, the central quantity is the frequency- 
and wavenumber-dependent scattering function P(k,cj) - 
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nOO ■ 

e"^ f (k, f) df . The scattering cross section corresponding 
to an energy transfer ho) and a momentum transfer is then 
essentially given by Re[P(k, w)], known as the (incoherent) 
dynamic structure factor [11]. In the case of simple diffusion, 
the dynamics is represented by a simple pole on the negative 
imaginary axis. 



An equivalent way to characterise the dynamics of a tracer 
is to give a prescription on how individual trajectories are 
generated as a result of the stochastic fluctuations in the 
medium [15, 16]. The propagator P{r,t) is then the result 
of a suitable average over the possible individual realisations 
of the jittery motion for the tracer. The equations of motion 
naturally become stochastic differential equations referred to 
as Langevin equations, which incorporate the randomness of 
the kicks by the medium as "noise". The modern formula- 
tion in terms of a Langevin equation is mostly due to Ornstein, 
who shaped the notion of what is now known as random gaus- 
sian white noise. A mathematical rigorous introduction to the 
stochastic differential equations and Brownian motion can be 
found in the excellent textbook by 0ksendal [17]. 

For overdamped motion, the displacements R(f) are as- 
sumed to obey the stochastic differential equation 

5,R(f) = /,(f), (12) 

with noise terms //(f) = {r]i(t), . . . ,Tid{t)) that are consid- 
ered as independent, random quantities on sufficiently coarse- 
grained time scales. In fact, these rjiit), i = l,...,d, rep- 
resent already averages over many independent processes oc- 
curring on even shorter time scales such that the central- 
limit theorem applies. The probability distribution then cor- 
responds to a multivariate gaussian, symbolically written as 
P[r](t)] cc D[Ti(t)] exp (- / df T](tf/4D), which is characterised 
completely by the only non- vanishing cumulant [15], 

{r]i(t) Jlj(t')) = 2DSij6(t - f') , i,j^l,...,d. (13) 

Such a noise displays only short-time correlations and corre- 
sponds to a power spectral density that is flat at the frequencies 
of interest, commonly referred to as white noise. We have im- 
posed that different Cartesian directions 77,(f) are uncorrected 
and, invoking isotropy, the strength of the noise, 2D, is iden- 
tical for all directions. The idea of coarse-graining and the 
seemingly innocent assumption of independence then neces- 
sarily leads to gaussian white noise as the universal law for the 
statistics of the displacements at small times. Any deviation 
from this law indicates the existence of non-trivial persistent 
correlations in the system. 

The displacement after a finite lag time follows from for- 
mally integrating the Langevin equation, AR(f) = R(f) - 
R(0) = ^df' j;(f'), and being a sum of gaussian variables, it 



obeys again a gaussian distribution. Thus it suffices to cal- 
culate the first two cumulants. Since the mean of the noise 
vanishes, one infers (AR(f)) = 0, and the correlation function 
of the displacements follows from the delta-correlated noise as 

{ARi(t) ARj(t)) = IDtSij . (14) 

In particular, one recovers 6r^(t) = ^AR(f)^^ - IdDt, and the 
probability distribution is determined by the diffusion propa- 
gator, Eq. (2). 

2.2. Anomalous and complex transport 

The probabilistic reasoning presented in the previous subsec- 
tion suggests that normal diffusion emerges as a statistical 
law essentially by the central-limit theorem. In particular, 
the mean-square displacement is expected to increase linearly 
in time for time scales much larger than microscopic ones. 
In simple systems such as normal liquids [11, 12] one ob- 
serves diffusion already at time scales exceeding the picosec- 
ond scale. The phenomena of anomalous or complex trans- 
port deal with dynamics where this diffusive regime is not vis- 
ible even on time scales that are by many orders of magnitude 
larger than picoseconds. Conventionally, a non-linear growth 
of the mean-square displacement 6r^{t) is taken as indicator 
of such unusual behaviour Typically, the mean-square dis- 
placement is proportional to a power law, 6r^{t) oc t" , with an 
exponent < a < 1. Hence the mean-square displacement 
increases slower than for normal diffusion, formally the dif- 
fusion coefficient becomes zero, nevertheless the tracer is not 
localised. This kind of behaviour is referred to as subdiffusion 
or anomalous transport^ . Theoretically, the phenomenon then 
calls for reasons why the central-limit theorem does not apply 
at the time scales of interest. Rephrasing the argument in terms 
of increments reveals that persistent correlations are hidden in 
the dynamics on meso- or macroscopic time scales. 

We would like to make a distinction between a simple vi- 
olation of the central-limit theorem in some intermediate time 
window and mechanisms leading to subdiffusive behaviour 
that can in principle persist forever. In the first case some dy- 
namic processes are unusually slow that spoil the central-limit 
theorem on these scales, yet ultimately normal transport sets 
in. This scenario of complex transport occurs generically by 
having constituents of the medium of different sizes or soft 
interactions, e.g., polymers. Then the mean-square displace- 
ment displays only a crossover from some short-time motion 
to long-time diffusion. Since the crossover can extend over 
several decades (due to a series of slow processes occurring 
in the medium), fits by power laws are often a satisfactory de- 
scription. In the second case, the correlations in the incre- 
ments decay slowly and upon tuning suitable control parame- 
ters the window of subdiffusion can become arbitrarily long. 

' In different contexts one finds also superdiffusive transport corresponding 
to ff > 1 , which is beyond the scope of thus review. 
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log(r) log(g 
log(time) 



Figure 1. Schematic mean-square displacement (MSD) for interme- 
diate subdiffusive transport. Free diffusion at microscopic scales is 
followed by subdiffusive transport at intermediate time scales. In a 
physical system, the subdiffusive growth ends typically at a second 
crossover, where the MSD grows linearly again with reduced diffu- 
sion constant, D <k Do, or where it saturates, e.g., due to boundaries 
like the cell membrane. 



Hence in a well-defined limit the subdiffusion persists forever 
and the central-limit limit theorem never applies. We reserve 
the term anomalous transport for the latter scenario. Typi- 
cally, the mean-square displacement is expected to display two 
crossover time scales, see Fig. 1 and Ref. 18, which is also 
found in experiments [19]. The long-time diffusion coefficient 
is then strongly suppressed compared to the microscopic mo- 
tion at short times, for example, subdiffusion with a - 0.6 over 
4 decades in time yields a reduction of D over Do by a factor 
of (ft/r)'"" = lo"*^""* ~ 40. We postpone the discussion what 
physical mechanisms can lead to such drastic changes and fo- 
cus here on the measurable quantities suited to reveal complex 
and anomalous transport. 

2.2.7. Van Hove self -correlation function 

The basic observable is the fluctuating single-particle density, 
p(r, f) = 5(r - R(f)), and the corresponding correlation func- 
tion, P(r - r', f - f') - V {p(r, f) p(r', f')), is referred to as van 
Hove (self-)correlation function [11]. Here V is the volume 
of the container and the thermodynamic limit V ^ oo is an- 
ticipated. Furthermore by translational invariance and for a 
stationary stochastic process, P depends only on the elapsed 
time t - t' and the accumulated displacement r - r'. Then 
the van-Hove correlation function can be cast into the form 
P(r, t) - {5{r - [R(f) - R(0)])) which is interpreted directly 
as probability density for an observed displacement r after a 
lag time /. Furthermore for isotropic systems, to which we re- 
strict the discussion, only the magnitude r - |r| enters the van 
Hove function. 



The first generalisation of simple diffusion consists of as- 
suming that the increments AR(f ) follow a gaussian probability 
distribution with zero mean. 



gauss 



s(r,f) 



1 



\ 



[27T5rHt)/df' ^'''^\2drHt)/dj 



(15) 



Here the mean-square displacement 6r^(t) characterises the 
width of the distribution; for ordinary diffusion, it holds 
6r^{t) = 2dDt, of course. The higher moments. 



5r«(f) := <|AR(f)|"> 



r |r|" P(r, t) , 



(16) 



are obtained by performing gaussian integrals, for example, 
one finds for the mean-quartic displacement, 5r'*(f) = [{d + 
2)/d][6r\t)f. 

Equivalently to the van Hove function, one can study 
the single-particle dynamics by monitoring the decay of den- 
sity fluctuations in the wavenumber representation p(k, f) - 
exp(ik-R(f)). The coiTesponding coiTelation function, P{k, t) - 
{p(k, f)*p(k, 0)), is the (self-)intermediate scattering function 
and merely the spatial Fourier transform of the van Hove func- 
tion. Again by isotropy, P{k, t) depends only on the magni- 
tude k - |k| of the wavenumber. The explicit representation 
P{k, t) - {exp[-ik • AR(f)]) permits an interpretation of P{k, t) 
as the characteristic function of the random variable AR(f), 
such that the lowest order moments can be obtained from a 
series expansion for small wavenumbers, k Q, 



P(k, t) = 1 



— 6r^(t) + 

2d Mid + 2) 



6r\t) + 0{k^). (17) 



Here we used that the orientational average over a d- 
dimensional sphere yields- cos^ = Ijd, cos* § - 3/d(d + 2). 

The logarithm of the characteristic function generates the 
cumulants. 



In P(k,t) 



2d 



6r^(t) ■ 



k*[6rHt)f 



6r\t) 



d + 2 



[6rHt)f 



+ Oik^), k^O. (18) 



In the case of gaussian transport, the Fourier transform of 
Eq. (15) yields P^^^ssikj) = exp(-k^6r^(t)/2d), and ah cu- 
mulants apart from the second one, 6r^{t), vanish identi- 
cally. Thus a simple, dimensionless indicator for transport be- 

' Let iij be the Cartesian components of a unit vector. Then by symmetry one 
argues that 

njn] = ASij 
niUjUkni = BiSijSki + SiuSji + SnSjk), 

where — indicates a spherical average. Contracting the indices i and j in 
the first relation reveals A = contracting in the second shows B = 
l/d{d + 2). Thus 4 = 3/d{d + 2). 
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yond the gaussian approximation is the non-gaussian parame 
ter [11, 12], 

d 6r\t) 



aiit) :-- 



d + 2 [6rHt)f 



1 



(19) 



The subscript indicates that there is an entire series of similarly 
defined quantities involving higher moments of the displace- 
ments. The inequality (X^) > (X)^ for the random variable 
X = |AR(f)|^ implies a lower bound on the non-gaussian pa- 
rameter, Qf2(f) > -2/(d + 2). 

The fact that probability theory imposes certain constraints 
on P{k, t) as a function of wavenumber k naturally poses the 
question if additional conditions apply if P(k, t) is consid- 
ered as a function of time t. More generally, what class of 
functions are permissible for correlation functions? Decom- 
pose the fluctuation density into Fourier modes p7(k, w) - 
j^j^2 6"^'p(k, t) for real frequencies oj and a long but fi- 
nite observation time T > 0. Then the corresponding power 
spectral density is obtained via the Wiener-Khinchin theo- 
rem [15, 16], 

^im 1 {\^pt(K ojf^ = 2Re[P(k, co)]. (20) 

Thus Re[P(k, oj)] > and inversion of the one-sided Fourier 
transform yields 



P(k. 



•"=-r 

Jo 



Re[P(k,w)] cos(wf)dw. 



(21) 



i.e., the propagator is the cosine transform of a non-negative 
function. Transforming again to complex frequencies by 
one-sided Fourier transform one derives relations of the 
Kramers-Kronig type [20, 21]; in particular, one observes that 
Re[P(k, oj)] > not only for real but for all complex frequen- 
cies in the upper half-plane, Im[w] > 0. 

Since the particle is to be found somewhere a parti- 
cle conservation law holds, and the intermediate scatter- 
ing function approaches unity in the long-wavelength limit, 
limk^o ^'(k, f) = 1 ■ For the one-sided Fourier transform, the 
particle conservation law suggests the representation P(k, oj) - 
1 /[-ia) + k^D(k, oj)], where D{k, a>) is known as the frequency- 
and wavenumber-dependent diffusion kernel [11, 12]. From 
Re[P(k, oj)] > for Im[(w] > the same property is inherited 
for the diffusion kernel, Re[D(k, a>)] > 0. 

Of particular interest is the long-wavelength limit Z{a)) - 
D{k — > 0,0)), which encodes the spatial second moment of 
the tracer motion. Note again Re[Z(tij)] > for Im[dLt] > 0. 
Expanding for small wavenumber. 



that Z(w) is the one-sided Fourier transform of the velocity 
autocorrelation function (VACF), 



Z(0 = i<v(0.v(0)> = l^..^(0. 



(23) 



Reversely, the mean-square displacement is obtained by inte- 
gration, 

6r^{t) = 2d f dt' (f - t')Z{t') . (24) 
Jo 

For stochastic processes where the derivative of the increments 
AR(f) does not exist, e.g., for a Brownian particle, the VACF 
may be defined via the mean-square displacement and can 
be shown to be a negative and completely monotone func- 
tion [22]. 

In the case of ordinary diffusion, 6r^(t) = 2dDt, the diffu- 
sion kernel simply assumes a constant, D(k, a>) i-> Z(cj) = D, 
and the velocity decorrelates instantaneously, Z(f) i-» D6(t — 
0). Furthermore, the non-gaussian parameter vanishes identi- 
cally, a'2(f) = 0, as do all higher cumulants. 



2. 2. 2. Distribution of squared displacements 

The measurements of the squared displacements AR(f)^ for 
a single particle along its trajectory often does not represent 
well the mean-square displacement 5r^{t), rather a significant 
scattering of the data is observed. This observation suggests 
to introduce the distribution function for the squared displace- 
ments u, 

p(u,t):^{6[u-AR(t)^)). (25) 

The probability distribution is obviously properly normalised, 
p(u,t)du - 1, its mean reproduces the mean-square dis- 
placement, and higher moments yield the even displacement 
moments, 

u''p(u, t) du . (26) 

The fluctuations in the squared displacements are given by the 
second cumulant and are via Eq. (19) already encoded in the 
non-gaussian parameter. 



Var[AR2(f)] = (|AR(f)|'*) - (AR(f)')' 



d + 2 2 
d d 



[Sr\t)]\ 



(27) 
(28) 



Rather than dealing with the probability distribution it is 
often favourable to work with a moment-generating function, 
which is the Laplace transform of the probability distribution, 



P{k, w) = 



1 



1 k^Z{cj) 
-ioj + k^D(k, cj) -icj {-ioj)^ 



Xoo 
due-"^"-'p(u,t) = (exp(-AR(f)Vw2)) , (29) 



and comparing to Eq. (17), one 
-(a//2d) J^dt e'"'(5r2(f). Integration 



finds Z{oj) - where the second representation follows from Eq. (25). The 
by parts reveals convention is chosen such that is the variable conjugate 
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to AR(f)2, and w carries the dimension of a length. Since 
the exponential is approximated by unity for small displace- 
ments, AR(f)^ w^, and rapidly approaches zero for large 
ones, AR(f)^ » w^, M(w, t) essentially constitutes the proba- 
bility for the particle to be still or again within a distance w of 
its initial position. In Section 4.2, it will be shown how this 
quantity can be directly measured by fluorescence correlation 
spectroscopy (FCS), where w corresponds to the beam waist 
of the illuminating laser [23]. 

The van Hove correlation P(r, f) is the probability distribu- 
tion of all vector displacements r, and the probability distribu- 
tion of the squared displacements p(u, t) follows by marginal- 
ising, 

p{u, t) = J d''r 6{u - r^) P(r, t) . (30) 

For the important case of statistically isotropic samples, P{r, t) 
depends only on the magnitude of the displacement r - |r| 
and the integral can be evaluated. In spherical coordinates, the 
angular integration yields the surface area of the li-dimensional 
unit sphere, Qj, as a factor The radial integral collapses due 
to the Dirac delta function, 6{u - P')- 6{y/u - r)/2u, and one 
obtains 

p(u, t) = ^u'''-^^'^P{r = V^, f) ■ (31) 

For gaussian transport, solely characterised by the time- 
dependent mean-square displacement, dr^it), Eq. (15) yields 

the gamma function evaluates to Y{dl2) - -^r, 1, V^/2 for 
d - 1,2,3. From these expressions one readily calculates also 

rU 

the cumulative distribution functions, J^^ p{u, f) du. 

The van Hove function P{r, t) and its spatial Fourier trans- 
form, the intermediate scattering function, P(k, f), the distribu- 
tion of the squared displacements p{u, f), and the correspond- 
ing moment generating function M{w,t) all encode spatio- 
temporal information on the motion of the tracer For isotropic 
systems they are all equivalent in principle, in practice they 
are sensitive to different aspects of transport occurring on dif- 
ferent length scales. Correlation functions that involve more 
than two times provide even more information on the dynam- 
ics, and may hold the key to distinguish different theoretical 
models that yield the same two-time correlation functions. 

3. Theoretical models 

The paradigm of anomalous transport is tantamount with a vi- 
olation of the central-limit theorem on arbitrarily long time 
scales. Modelling such processes requires including persis- 
tent correlations that manifest themselves as self-similar dy- 
namics in the mean-square displacements. Different models 
and theoretical approaches have been pursued that genericaUy 



lead to subdiffusion. Here we focus on the three most widely 
used frameworks. The perhaps simplest approach is based on 
stochastic differential equations where the noise term displays 
persistent correlations which then transfer to the increments. 
Since usually the statistics of the noise is still assumed to be 
gaussian, they differ essentially only in the form of the tem- 
poral correlations they incorporate and we summarise them as 
gaussian models. The second category, the continuous-time 
random walk (CTRW), consists of jump models where parti- 
cles undergo a series of displacements according to a distribu- 
tion with large tails. Here the central-limit theorem does not 
apply since the mean waiting time for the next jump event to 
occur becomes infinite. The last class of Lorentz models relies 
on spatially disordered environments where the tracer explores 
fractal-like structures that induce anomalous dynamics. 

3.1. Gaussian models 

Here we collect properties of a class of models that give a phe- 
nomenological description of complex and anomalous dynam- 
ics which all result in gaussian propagators. 

3.1.1. Fractional Brownian motion 

A simple model for subdiffusion is fractional Brownian mo- 
tion, introduced rigorously by Mandelbrot and van Ness [24] 
as superposition of Brownian processes with power-law mem- 
ory. Here we follow a heuristic approach [25-27] that sum- 
marises the essence of fractional Brownian motion. Assum- 
ing the stochastic differential equation (12), (9,R(f) = i/(f), we 
have already seen that if the noise ///(f) is delta-correlated in 
time, the mean-square displacements increase linearly. If ad- 
ditionally the noise obeys a gaussian statistics, this property 
is inherited also for the displacements AR(r) and transport is 
completely characterised by the diffusion propagator P(r, f), 
Eq. (2). The idea is now to incorporate persistent correlations 
in the noise such that transport is drastically slowed down with 
respect to normal diffusion. Since the noise plays the role of a 
fluctuating velocity, we use the same notation as in Section 2.2, 
to express the noise correlator in terms of the velocity autocor- 
relation function, 

{mit) rij(t')) = dZ(\t - f'l) 6;j , i,j^l,...,d. (33) 

Here the different Cartesian components are taken as uncor- 
related, which certainly holds for an isotropic system. In the 
Fourier domain, this implies 

(77,(w)* T]j(aj')) = 4Kddij 6(0) - oj') Re[Z(w)] , (34) 

/"CO ■ 

where Z(lo) - e Z(t)dt is again the one-sided Fourier 
transform of the velocity autocorrelation function. For ordi- 
nary diffusion, Z(a)) = D is constant and the noise corresponds 
to white noise. 
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The case of subdiffusion, 6r^(t) = IdKat" with an expo- 
nent < Q- < 1 and a generalised diffusion coefficient Ka > 0, 
then yields a spectral density Z(w) = (-iw)'"°'A'Q,r(l +a) [24]. 
Hence the strength of the noise approaches zero as the frequen- 
cies become smaller, which explains that transport slows down 
with increasing correlation time. In the temporal domain, Z(f) 
is represented by a pseudofunction ', 

Z{t)^a{a-l)KaPf\tr^. (35) 

Up to here only the mathematical frame has been set and no 
assumptions on the nature of the stochastic process has been 
made. 

In fractional Brownian motion, the statistics of the noise 
correlator is assumed to be characterised by the only non- 
vanishing cumulant Z(f), Eqs. (33) and (35), i.e., a stationary 
gaussian process although not white noise. Then the statistics 
of the increments AR(f) is again gaussian, and the propagator 
reduces to Pgauss(r, f)> Eq. (15). Its scaling form corresponds 
to that of simple diffusion, Eq. (3), 

PFBMir, t) = r-'r^,,,,{f) , r cc rr"/\ (36) 

sharing the gaussian scaling function, Eq. (4), but not the 
scaling variable, r. In particular, the non-gaussian parameter 
Q'2(f) = vanishes by the construction of fractional Brownian 
motion. 

In contrast to simple diffusion, fractional Brownian motion 
is not a Markov process; in particular, the van Hove correlation 
function is not sufficient to characterise the statistical proper- 
ties completely. Multiple-time correlation functions encode 
non-Markovian behaviour, for which fractional Brownian mo- 
tion makes detailed predictions. As an example, 

{[R(t)-R(0)f[R(t + T)-R(T)f) 

= AcfKy" + 2dKl{\t + T\" + \t- T\" - 2T"f , (37) 

which has been derived recently to study the ergodic properties 
of fractional Brownian motion [29]. For T -Q, this expression 
reproduces the quartic moment, 5r'*(f), and is compatible with 
a vanishing non-gaussian parameter, Eq. (19). 

3.1.2. Langevin equations for visco-elastic media 

The erratic motion of a spherical particle immersed in a com- 
plex medium can be described quite generally by a Langevin 
equation. Rather than directly addressing the displacements 
one may base the description on the velocity \(t) = R(r) and 
formulate a force balance equation. The paradigm has been 
given by Langevin [3] himself, 

OTv(f) = -^v(0 + m , (38) 

^ The noise correlator corresponds to a distribution, and pseudofunction 
means that integrals with test functions </>(?) extract only Hadamard's fi- 
nite part [28], / if(t) Pf \t\"-- df := / - i/3(0)]|r|"-2 df. In particular, one 
easily verifies that the one-sided Fourier transform of Z{t) yields Z{(jJ). 



where m denotes the mass of the particle, the deterministic 
friction force —(\{t) is merely the Stokes drag, and f(f) is a 
fluctuating force with zero mean, (fiit)) - 0. The friction 
constant constant, ^ = 67TT]a, is directly connected to the sol- 
vent viscosity 77 and the particle radius a. The statistics of 
the random forces f(f) is characterised completely by the only 
non-vanishing cumulant [5], 

{fi(t) fj(t')) = 2kj,T( 6ij 6(t -t'), (39) 

where T is the temperature of the environment and de- 
notes Boltzmann's constant. Thus the Cartesian components 
fi(t) of the forces are gaussian distributed and independent for 
different times. The variance at equal times is again dictated 
by the fluctuation-dissipation theorem, see, e.g., Ref 15 for 
details. The delta-correlation in the temporal domain for the 
forces translates to white noise for the corresponding power 
spectral density. The velocity autocorrelation then decays ex- 
ponentially [3], 

(v,(f) v/f')) = {k^TIm) 6ij exp(-|f - f'l/r^) , (40) 

where = m/^ = m/Gnria is the momentum relaxation time. 
Similarly, the mean-square displacement of the t/-dimensional 
motion is calculated to 

6r^(t) ^2dD[t + Tp (e"'^^" - l)] , (41) 

where D = k^T/^ is the diffusion constant according to the 
Stokes-Einstein relation. 

The description can be easily generalised for the case 
of visco-elastic media [30]. Here the response of the com- 
plex solvent to shear is encoded in the complex frequency- 
dependent viscosity, rjiai). In the conventions employed in this 
review, Re[77(w)] > corresponds to the dissipative part and 
Im[7;(a>)] encodes the reactive part. Equivalently, one may em- 
ploy the complex shear modulus G(cj) := -icorjiaj). For exam- 
ple, in the Maxwell model, G(w) - -ia)TMGoc/(l - iwr^f), 
the modulus is characterised by a high-frequency elastic re- 
sponse Goo and a crossover time scale tm- The Stokes drag 
in a visco-elastic medium then depends on the frequency and 
the Langevin equation is discussed conveniently in the Fourier 
domain [31], 

- iwmv((y) = -f(w)v(w) + f(w), (42) 

where ((oj) = 6nri(Lo)a replaces the Stokes friction coefficient. 
By the fluctuation-dissipation theorem, the force correlator 
has to be modified accordingly, 

{fAc^r fj(aj')) - 47r^BT-Re[^(c^)] 6ij 5(u - oj') . (43) 

Since the fluctuations arise in the surrounding solvent as a 
sum over uncorrelated regions, the forces are again assumed 
to be gaussian. Then it is clear that the van Hove correlation 
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function for the particle corresponds to a gaussian propagator 
and the dynamics is specified entirely by the mean-square dis- 
placement, 6r^(t). Rather than solving for 6r^(t), we solve for 
the one-sided Fourier transform of the velocity autocorrelation 
function [30,31], 



Z(w) 



-imoj + ({cl)) 



(44) 



Relying on the relation ^(w) = 67r?/(<y)fl, the local visco-elastic 
response of a complex medium is inferred from the motion of 
tracer particles in microrheology experiments [30]. 

Subdiffusion at long times is obtained if Z(w) = 
(-iajy^" KaTil + a) for a> ^ [25], i.e., the elastic modulus 
displays power-law behaviour, G(a)) ~ (-iw)", which appears 
to be generic in many biological materials and soft matter sys- 
tems for intermediate frequencies. This empirical observation 
is formulated in the soft glassy rheology model [32]. The mi- 
croscopic mechanism remains in general unspecified, yet for 
the case of a solution of semiflexible polymer networks, the 
bending rigidity of a single filament suggests an elastic power- 
law response, G{(jj) ~ (-iw)^^"* [33-36]. Similarly, by cou- 
pling to the elastic degrees of freedom of a membrane, effec- 
tive fractional friction kernels can be generated in the same 
way with various exponents depending on the level of descrip- 
tion of the membrane [37, 38] 

3.1.3. Long-time anomalies 

The assumption of an instantaneous friction term in the 
Langevin equation (38) is in fact incorrect even at long times, 
as has been pointed out already by Hendrik Antoon Lorentz. 
The reason is that the Stokes formula applies for steady mo- 
tion of the particle only and the friction is accompanied by 
a long-ranged vortex pattern in the velocity field of the en- 
trained fluid. For unsteady motion, the particle excites inces- 
santly new vortices diffusing slowly through the fluid. As a 
consequence the friction force depends on the entire history 
of the particle's trajectory, an effect known as hydrodynamic 
memory. The theoretical description is achieved most conve- 
niently in the frequency domain. The drag force for a sphere 
performing small-amplitude oscillations of angular frequency 
0) has already been calculated by Stokes [39] and leads to a 
frequency-dependent friction coefficient [15], 

f(w) - 6nr]a (l + -\/-ia>Tf^ - icomf/l. (45) 

For steady motion, co = 0, the formula reduces to the Stokes 
drag. The last term appears as an acceleration force for half 
of the displaced fluid mass, mf - Anpfa^ /3, and it is natural 
to absorb this contribution by introducing an effective mass 
for the particle, meff = m + ;«f/2. The second modification 
is a non-analytic contribution due to the slow vortex diffusion 
in the liquid as the particle performs unsteady motion. The 
characteristic time scale, Tf = pia^/j], is the time needed for a 



vortex to diffuse over the distance of the radius of the particle. 
By the fluctuation-dissipation theorem, Eq. (43), the spectrum 
of the random forces is no longer white but displays a coloured 
component that increases as a square root with frequency. Re- 
cently, the power spectral density of the thermal noise has been 
measured experimentally for a single bead by ultra-sensitive 
high-bandwidth optical trapping [40] in excellent agreement 
with theoretical predictions. 

The velocity autocorrelation function in the frequency do- 
main, Eq. (44), acts as an admittance or frequency-dependent 
mobility and displays a non-analytic low-frequency expansion. 



Z(w) = D 



1 



V-it^Tf + (9(w)j. An explicit expression in the 
temporal domain is achieved in terms of error functions [41], 
here we focus on the long-time anomaly. 



Z(f) 



D 



(46) 



which is a direct consequence of the non-analytic terms in 
((ol>). The most striking feature is that Z(t) encodes persistent 
correlations manifested by a self-similar tail in strong contrast 
to the exponential decay of Langevin's original theory. These 
long-time tails have been discovered first in computer simula- 
tions for fluids [42, 43] and, only recently, have directly been 
observed for colloidal particles in suspension [44-47]. The 
mean-square displacement follows directly by integration, 

6r^(t)^6Dt[l-2^fi^t + 0{t-^)]. (47) 

The algebraic tail in the velocity autocorrelation manifests it- 
self in a slow approach to normal diffusive transport. 

The persistent correlations in the mean-square displace- 
ment, buried under the leading linear increase, show that the 
assumption of independent increments is not satisfied and that 
the regime of truly overdamped motion is never reached due 
to the hydrodynamic memory, even at long time scales. Nev- 
ertheless the central-limit theorem remains valid, although the 
convergence is slow due to the persistent power-law correla- 
tions induced by vortex diffusion. 

3.2. Continuous-time random walks ( CTRW) 

A different class of models that is widely discussed is the 
continuous-time random walk (CTRW) [48-51], originally in- 
troduced by Montroll and Weiss [52] for hopping transport 
on a disordered lattice. The particles spend most of the time 
bound to a trap with an escape time that depends sensitively 
on the depth of the trap. Rather than dealing explicitly with 
the quenched disorder on the lattice, the medium is treated 
as homogeneous with the new ingredient of a waiting-time 
distribution for the next hopping event to occur Anomalous 
transport can be generated within this framework by assum- 
ing waiting time distributions such that the mean waiting time 
becomes infinite. The central-limit theorem does not apply 
since longer and longer waiting times are sampled. It turns 
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out that the CTRW as mean-field approximation to hopping in 
the quenched trap model gives the same result in dimensions 
2 and higher [49]. The trapping in biology seems rather natu- 
ral, due to chemical attachments of molecules in the cell, such 
that binding on broadly distributed time scales may lead to 
CTRW dynamics. First we provide the general description of 
the CTRW model, then we show how subdiffusion can emerge 
and discuss scaling properties of the propagator. 

3.2.1. Model definition 

In the CTRW model, the particle is assumed to traverse the 
(/-dimensional space by a series of jumps such that the dis- 
placement, {, and the waiting time to perform the next jump, 
f, are drawn from a given distribution, ifr{{,t). For simplic- 
ity, we assume that the observation starts when the process 
is initialised. Hence, a CTRW belongs to the broad class of 
renewal-reward processes. In particular, the propagator P{r, t) 
fulfils a renewal equation [53], which follows by conditioning 
on the event that the particle has accumulated a displacement 
r at time t within the first step. If the first step occurs later than 
at time t the propagator is simply 6{r), otherwise the process 
is renewed [48, 52, 54, 55], 

P{r,t) = ^(r) df' J d''eifr(e,t') 

+ J^dt' Jd''£ P{r -tj-t') ifr{e, t') . (48) 

The first term corresponds to the propagator provided that no 
jump occurred and can be rewritten as 



Po(r,0 = 5(r) 



1 



r df'(A(f') 

Jo 



(49) 



where i]/{t) - J if/if, t) d''( is the jump probability density irre- 
spective of the size of the jump. The solution of the renewal 
equation is most easily achieved after a spatial Fourier trans- 
form, f (k, t) - J d'^r e P{r, t), and a subsequent temporal 
one-sided Fourier transform, P(k, w) = J^°°df e"^'P(k, f) for 
Im[(jj] > 0. By the convolution theorem the renewal equation 
simplifies to f (k, w) = Po(k, w) + w)f (k, a>). The spatio- 
temporal Fourier-Laplace transform of Eq. (49) can be be per- 
formed directly, P()(k, w) = [l-tfr(a>)]/(-ia)). Combining both 
results yields the Montroll-Weiss relation [48, 54, 56-58] for 
the propagator in terms of the jump probability distribution. 



P(k,w) 



1 1 - (A(w) 
-io) 1 - i/'(k, cj) 



(50) 



In many applications the jump distribution is not known, 
and additional assumptions are necessary to define the model. 
First, we assume that the jumps exhibit no preferred direction 
such that if/(f, t) depends only on the magnitude \(\, Hence, the 
dynamics becomes isotropic and P(k, to) is a function of the 



wavenumber, k - |k|, only. Second, the waiting-time distribu- 
tion is often taken to be independent of the jump size, and the 
corresponding jump distribution factorises, t) - A(t) il/{t) 
such that jA({)d''£ = 1. Then the propagator assumes the 
simple form 



P(k, oj) 



1 1 - i//(oj) 
-ioj 1 - A(k) (//(a)) ' 



(51) 



which will be the starting point for the discussion. 

We shall assume that the jump size distribution A(() is 
well-behaved, in particular it decays rapidly for large dis- 
tances. Then the characteristic function of the jump sizes, A(k), 
encodes all moments. 



{t'}+0{lf)^ (52) 



For the mean-square displacement, only the second moment 
of the spatial distribution is relevant, {£^)/d - d^/l(A:)/dA:^|t=o, 
and, by Eq. (22), the velocity autocorrelation function is com- 
puted to 

-iwiA(w)<^^> 



Z(w) 



1 



(53) 



(//{cj) 2d 

By normalisation of the waiting-time distribution, it holds 
tl/(a> — > 0) = 1, and for a well-behaved distribution the 
low-frequency expansion provides the moments of the waiting 
time, (t) :- ttfi{t)dt. 



>^(lo) = 1 -I- iw <r) + 0(a/) . 



(54) 



In this case, the diffusion coefficient D — Z(a) — > 0) is finite 
with value D = {e^)/2d{T}. 

3.2.2. Anomalous transport 

Anomalous transport is obtained if the jump rate distribution 
is non-analytic at zero frequency, e.g.. 



(A(w) - 1 - i-icjT)" + h.o.t., 



0, 



(55) 



with non-integer exponent, < a < 1, some time scale t > 0, 
and neglecting higher order terms. Then the diffusion coeffi- 
cient vanishes, D - Z(w — > 0) = 0, and the VACF inherits a 
leading non-analytic contribution for — > 0, 



Z(w) ^ {-lajy-" K^Til +a). 



(56) 



with KaT(l + a) - {fi)T^'' I2d. By means of a Tauber theo- 
rem [53, 59], the non- analytic ity at zero frequency corresponds 
to a negative long-time tail in the VACF, 



Z(f) ^ -a{\-a)Kj"- 



(57) 



Thus, persistent anticorrelations become manifest in the 
VACF. Integration over time shows that the time-dependent 
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diffusion coefficient, D{t), approaches zero at long times only 
algebraically, D{t) ^ aKat"^^ for f — > oo. Eventually, one 
finds for the mean-square displacement. 



6r^{t) ^ IdKaf, t 



(58) 



i.e., a subdiffusive increase that persists for arbitrarily long lag 
times. Reversely, the low-frequency singularity in tff(a>) is con- 
nected to a power-law tail in the waiting-time distribution, 



(59) 



such that even the mean waiting time is infinite. 

Let us briefly discuss the next higher moment, the mean- 
quartic displacement. The expansion of the incoherent dy- 
namic structure factor, Eq. (51), in powers of the wavenum- 
ber yields the one-sided Fourier transform of the mean-quartic 
displacement, cf. Eqs. (17) and (52), 



Jo 



'dte^'6r\t) = md + 2)^^ + 2/^^^. (60) 



(-iw)3 <^2)(-iw)2 



In particular, the leading low-frequency singularity is again 
completely governed by the tail of the waiting-time distribu- 
tion. 



Jo 



df e'"'(j/(f) ^ Mid + 2)Klr{l + afi-ico)-^-^" ; (61) 



in particular, the 4th moment of the jump size, {^^), does not 
enter the leading order Application of a Tauber theorem [53, 
59] yields the asymptotic long-time behaviour. 



6r\t) ^ Mid + 2) l!;y"\ Ky'^, t^^. 
1(1-1- 2a) 



(62) 



One concludes that for a CTRW with power-law distributed 
waiting times, the non-gaussian parameter approaches a con- 
stant value. 



2r(l+af 



(63) 



irrespective of the jump size distribution A{(). The method 
can be easily extended to higher moments, and one can show 
that all non-gaussian parameters assume non-vanishing long- 
time limits, with values that depend only on the exponent of 
subdiffusion, a. 

Equations (51) to (53) imply a small-wavenumber approx- 
imation of the propagator. 



Pik,co) ^ 



1 



-ia> + k^Z((L)) ' 



(64) 



Propagators of this form will be referred to as generalised hy- 
drodynamics approximation. It allows for an exphcit solution 



of the spatial inverse Fourier transform; actually, all that needs 
to be done is to replace the diffusion coefficient D in Eqs. (10) 
by its frequency-dependent generalisation, Z(a)). For example 
in three dimensions, one finds 



1 



47rZ(w)r 



exp 



id^3). (65) 



3.2.3. Scaling limit and the fractional Fokker-Planck 
equation 

Let us discuss a simple, possibly the simplest distribution for 
the waiting time that displays a tail, the Cole-Cole distribu- 
tion,"^ 



1 



1 + i-icoT)" 



(66) 



In the temporal domain, this corresponds to the waiting-time 
distribution [55, 60] 



m = --£„(-(f/T)"), 



(67) 



where E„(-) denotes the Mittag-Leffler function. Inserting 
in Eq. (53), the velocity autocorrelation function assumes a 
power law for all frequencies. 



Z(a;) = ^(-iwT)'-^ 
2dT 



(68) 



and the corresponding mean-square displacement is subdiffu- 
sive for all times, 6r^(t) - 2dKaf'. The non-gaussian parame- 
ter can be evaluated exactly, 

= -F?TTTT-l + 7T^7^r(l + a) - . (69) 
r(l-l-2a) d + 2/£2y \tJ 

Within the generalised hydrodynamics approximation, the 
Fourier-Laplace transform of the propagator, P{k, oj), is given 
by Eq. (64) for long wavelengths, A; — > 0. Together with 
Z{oj) given by Eq. (68), this form corresponds to the scaling 
limit [55] after coarse-graining on large length scales and long 
times. The corresponding propagator in real space and time 
then is the solution of the so-called fractional Fokker-Planck 
equation [48], which is the scale-free limit of CTRW [55, 60]. 
The approach presented here devoids the introduction of frac- 
tional derivatives which naturally emerge if the equations of 
motion are introduced in the temporal domain. We refer the 
mathematically inclined reader to the excellent review by Met- 
zler and Klafter [48]. 

Such a frequency-dependence was introduced empirically and has since 
been widely used to describe the stretched dielectric response of polymeric 
liquids and glassforming materials, see Kenneth S. Cole and Robert H. 
Cole, J. Chem. Phys. 9, 341 (1941). 
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Figure 2. Scaling functions of the propagators for CTRW, frac- 
tional Brownian motion (FBM), and lattice percolation on the in- 
cipient infinite cluster in rf = 2 dimensions for subdiffusive motion 
with a = 0.695 on double-logarithmic (main panel) and linear (inset) 
scales. Scaling functions and scaling variables are normalised such 
that fij £'dxP{x)/x = 1 and Slj f^dixxPix) = d. 

The spatial backtransform is obtained from Eq. (9) upon 
replacing D by the above result for Z(a>). It is conve- 
nient to introduce a dimensionless, rescaled frequency, cb :- 

{{(^)l2dy^^"T r^l" oj. Then, 

P{v, c.) = (-ia>)«<''-2)/4 ^ (i-ia>rl^) . (70) 

The evaluation of the complex modified Bessel function of the 
second kind, Ky{-), as well as the temporal inverse Fourier 
transform, Eq. (8), can be achieved numerically, e.g., with 
Mathematica. In the fractional Fokker-Planck limit, the 
time- and space-dependent propagator exhibits a scaling prop- 
erty that holds for all times and distances, 

P{r, t) = r^'Pair) , f oc rr"!^ ; (71) 

the omitted prefactor, {fiy^^^r"!'^ , renders r a dimensionless 
scaUng variable. Although the scaling property is of the same 
form as for fractional Brownian motion, Eq. (36), the scaling 
functions are different. Here, Va depends on the space di- 
mension d and on the exponent a. The scaling functions are 
displayed in Figs. 2 and 3 for d - 2,3> and are compared to 
the scaling functions of fractional Brownian motion and ob- 
structed transport on the percolating cluster for the same sub- 
diffusion exponent. Note that Eq. (70) was derived previously 
on a lattice [61, Eq. (20)], the subsequent inverse Laplace 
transform to the time domain, however, relied on asymptotic 
approximations of the Bessel function, yielding a scaling func- 
tion, Paix) = exp(-jt;2/(2-«')/2) [50, 61], which strongly dis- 
agrees with our numerical findings at small arguments. 

Since CTRW is not a Markov process, the propagator 
P(r, t) does not characterise the dynamics completely. Higher 
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Figure 3. Scaling functions of the propagators for CTRW and frac- 
tional Brownian motion in = 3 dimensions for subdiffusive motion 
with a = 0.515 on double logarithmic (main panel) and linear (inset) 
scales. The dotted, black line indicates the anticipated asymptote for 
transport on the percolating cluster, 'P„{f) ~ f''. 



correlation functions involving multiple times have been cal- 
culated [62] and the procedure has been generalised to arbi- 
trary order [63]. In the fractional Fokker-Planck limit, the 
dynamics can be written as a subordinated Brownian motion, 
i.e., the individual trajectories have the very same shape as 
for Brownian motion, yet the time the particle needs to follow 
these trajectories is provided by a different clock [64]. Hence 
to distinguish the CTRW model, respectively the fractional 
Fokker-Planck description, from fractional Brownian motion, 
one can use the statistical properties of the individual trajecto- 
ries and derive suitable measures highlighting the difference, 
see e.g. Refs. 65-67. 

The propagator describes the probability for a tracer to 
have moved a distance r after a lag time t. If the stochas- 
tic process is in equilibrium this ensemble average can be ob- 
tained also as a moving-time average, i.e., by following a sin- 
gle trajectory of a tracer For CTRW such that the mean wait- 
ing is strictly infinite, the process never reaches equilibrium, 
no matter how remote the process has been initiated in the 
past [53, 68-70]. Then, time correlation functions exhibit age- 
ing, i.e., they explicitly depend on the points in time of the 
measurements and not just the time lag inbetween. A closely 
related phenomenon is referred to as weak ergodicity break- 
ing and was introduced by Bouchaud [71] in the context of 
spin glasses. The reason is at the very heart of the concept 
of CTRW with power-law distributed waiting times, since as 
time proceeds, longer and longer waiting times are drawn from 
the distribution before the particle actually moves. In conse- 
quence, individual trajectories can differ strongly from the av- 
erage, or in other words, the average is not dominated by the 
typical reahsations of the trajectories. Hence, one has to be 
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careful when comparing predictions for individual trajectories 
to ensemble averages and one should calculate probabilities 
for the time averages of observable quantities [72-76]. 

3.3. Obstructed motion: Lorentz models 

The models discussed so far implicitly assume that the parti- 
cle motion occurs in a homogeneous medium; for example, the 
waiting time in a CTRW does not depend on the position. In 
contrast, the Lorentz model paradigmatically describes trans- 
port in a spatially heterogeneous medium and displays many 
facets of anomalous transport like subdiffusion, crossover phe- 
nomena, immobilised particles, and long-time tails. In the 
original version, H. A. Lorentz [77] discussed the motion of 
a ballistic particle which is elastically scattered off randomly 
placed, hard spherical obstacles to lay a microscopic basis 
for Drude's electric conductivity of metals. The medium is 
reminiscent of a Swiss cheese: it consists of a homogeneous 
phase supporting free particle transport, punched by the pos- 
sibly overlapping obstacles. Thus at high obstacle density, a 
tracer faces a heterogeneous environment characterised by a 
significant excluded volume and a highly ramified remaining 
space. Variants of the model include Brownian tracer parti- 
cles [22, 78-80] or correlated obstacles, e.g., by simply forbid- 
ding obstacle overlaps [78, 81-86]. Saxton [87] and later Sung 
and Yethiraj [83, 84] recognised the biophysical relevance of 
such obstacle models for protein transport in cellular mem- 
branes. 

The model is described by a number of parameters. The 
obstacle radius cr fixes the unit of length, and the unit of 
time T is specified in terms of the velocity v of a ballistic 
tracer or the microscopic diffusion constant Do in case of a 
Brownian tracer Then, the dimensionless obstacle density, 
n* - ncr', remains as the only interesting control parameter; 
n is the number of obstacles per volume. Alternatively, one 
may use the porosity or fraction of unoccupied, void volume, 
e.g., (f - exp(-ni;) for independently inserted obstacles of vol- 
ume V. 

3.3.1. Microscopic theories 

The ballistic motion of the tracer is randomised by the scat- 
tering and becomes diffusive at long times. Assuming un- 
corrected collisions, Lorentz [77] computed the diffusion co- 
efficient about a century ago to Do - va-j^nn* for d - 3; 
for d - 2, Lorentz-Boltzmann theory yields Do = 3vcr/16n* 
[88]. These results are valid only in the dilute limit, n* — > 0. 
Actually, subsequent collisions with the obstacles are persis- 
tently correlated due to the frozen disorder, and a proper treat- 
ment involves the resummation of so-called ring collisions. A 
systematic low-density expansion of 1/D within kinetic the- 
ory [88, 89] revealed non-analytic terms at the leading {d - 2) 
or next-to-leading {d - 3) order, supported by pioneering com- 
puter simulations [90, 91]. For a Brownian tracer, the short- 



time motion is diffusive already by construction and the long- 
time tail of the VACF is obtained at leading order in a low- 
density expansion [80], employing a scattering formalism for 
the Smoluchowski operator. 

The microscopic theoretical treatment beyond the regime 
of low densities could be extended using self-consistent ap- 
proximations. Gotze et al. [92, 93, 94] closed the Zwanzig- 
Mori equations for the incoherent dynamic structure factor 
by introducing a set of approximations for the frequency- 
dependent memory kernel; an approach that formed the ba- 
sis of the successful microscopic mode-coupling theory for 
the description of the glass transition [21]. For the Lorentz 
model, the theory predicts a localisation transition, i.e., a crit- 
ical obstacle density n*, where long-range transport ceases. 
Din*) - for n* > n*. As the critical density is approached, a 
temporal window emerges with a subdiffusive increase of the 
mean-square displacement; concomitantly, the long-time dif- 
fusion constant is predicted to vanish according to a power law, 
D(«*) ~ \n* - «*!'' for n* t n*- A similar picture was obtained 
within a self-consistent repeated-ring kinetic theory [95], al- 
though both theories have deficiencies and differ in their de- 
tailed predictions. 

Recently, transport of a dense fluid in a porous host struc- 
ture has attracted new interest. Significant progress in this 
direction has been achieved by Krakoviack [96, 97, 98] by 
generalising the mode-coupling theory of the glass transition. 
There, a series of non-equilibrium phase transitions was pre- 
dicted, which roughly speaking correspond to glass-glass tran- 
sitions in size-disparate mixtures [99]. Some of the predicted 
phenomena have already been observed in computer simula- 
tions [100-103]. In the limit of small fluid density, the model 
degenerates to the Lorentz model and the physics should be 
dominated by long- wavelength phenomena [104]; however, 
these features appear not to be correctly reproduced by the cur- 
rent theories. A theoretical framework combining all aspects 
of glassy dynamics and localisation is yet to be developed. 

3.3.2. Percolation and random-resistor networks 

A phenomenological approach to the transport properties of 
the Lorentz model relies on a mapping to random-resistor net- 
works. The medium is thought of as a set of voids connected 
by channels, and transport occurs via hopping from void to 
void. The geometry of the channels, formed by close by ob- 
stacles, imposes a distribution of transition probabilities be- 
tween the voids. For many voids, although being neighbours 
in space, a direct connection is blocked by obstacles and the 
transition probability is zero. The anticipated link to random- 
resistor networks is provided by representing the voids as the 
nodes of a regular lattice and interpreting the transition prob- 
abilities as electric conductances between lattice nodes. The 
macroscopic diffusion constant is then identified as the to- 
tal conductance of the network. For a given obstacle con- 
figuration, the network may be constructed rigorously from 
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a Voronoi tessellation [83, 84, 105, 106]. Increasing the ob- 
stacle densities corresponds to diluting the conductive bonds 
of the network, precisely as in the bond percolation prob- 
lem [50, 107]. If the bonds are sufficiently diluted, the network 
falls apart into differently sized clusters, and above a critical 
threshold and in the thermodynamic limit, no spanning cluster 
exists that supports macroscopic transport. 

Let us briefly summarise the essence of percolation theory; 
we refer the reader to the excellent textbooks by ben Avra- 
ham and Havlin [50] and Stauffer and Aharony [107] for a pro- 
found introduction. Directly at the percolation transition, the 
incipient infinite cluster is a fractal in a statistical sense. It is of 
inextensive weight and occupies a volume Scx,(L) ~ if' within 
a ball of radius L, defining the fractal dimension, c/f < d. A 
self-similar hierarchy of finite clusters coexists with the infi- 
nite cluster, whereby the distribution of cluster sizes s follows 
a power law, i ' "'^"'f. Away from the transition, the medium 
is no longer scale-free and self-similarity holds only on length 
scales below the correlation length, ^. In particular, the proba- 
bility to encounter finite clusters of linear extent larger than ^ 
is at least exponentially suppressed. On the percolating side 
of the transition, the infinite cluster looks homogeneous at 
scales larger than ^. The correlation length exhibits a non- 
analytic singularity at the transition in form of a power-law 
divergence, ^ ~ \n* - n'T'', which introduces a second crit- 
ical exponent, v. The macroscopic conductivity E of a per- 
colating random-resistor network vanishes at the percolation 
threshold with a power-law singularity as well, "E ~ \n* - n*f. 
The conductivity exponent /i constitutes a third critical expo- 
nent, which describes transport and dynamic phenomena and 
genuinely extends the set of geometric exponents (t/f, v); in 
general, it does not follow by means of a scaling relation from 
the other two. 

The percolation transition shares many aspects of a contin- 
uous thermodynamic phase transition and is tractable by renor- 
malisation group methods [107, 108]. An important lesson is 
that the critical exponents df, v, and fi are determined by a 
non-trivial fixed point of the renormalisation group flow, i.e., 
a systematic and consistent coarse-graining of the medium. 
As a consequence, the microscopic details of a specific model 
become irrelevant for the leading singular behaviour near the 
transition, and in this sense, the critical exponents are univer- 
sal numbers. The percolation transitions of many different lat- 
tice types belong to a single universality class and are charac- 
terised by the same set of exponents, which depends only on 
the dimension of space. A compilation of recent values for the 
critical exponents of random percolation is given in Table I. 

For transport on continuum percolation clusters, a pecu- 
liarity arises: the dynamic universality class may be different 
from that of lattice models. In the context of random-resistor 
networks it was shown that a sufficiently singular power-law 
distribution of weak bond conductances, no-(crbond) ~ ""bond 
for (Tbond — » 0, can dominate the renormalisation flow such 
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Table I. Static and dynamic percolation exponents for the leading and 
sub-leading critical behaviour. Numbers in parentheses indicate the 
uncertainty in the last digit. Sources: (a) Ref. 109, (b) Refs. 110 and 
111, (c) from exponent t = 1 + djdf = 2.186(2) [112], (d) Ref. 113, 
(e) continuum percolation theory yields d„ = d[ + 2/v for a ballistic 
tracer [79, 114, 115]. The values for y and z were calculated from 
exponent relations, Ww = i^^f [116] andz = 2d^/(2+d{-d) [50]. The 
dynamic universality class does not split for d = 2 [1 14, 1 17, 1 18]. 



that the conductivity exponent /u deviates from its universal 
value on lattices, /j. = maxjyu'"", v((i - 2) + 1/(1 - a)} [119- 
121]. In continuum percolation, the weak conductances arise 
from narrow gaps or channels connecting the voids. For uncor- 
rected spheres or discs, indeed a power-law distribution was 
derived, which is singular enough in three dimensions to mod- 
ify the exponent; the chain of arguments has been summarised 
in Ref. 79. Halperin et al. [117] predicted jj. - v -H 3/2 for the 
conductivity problem, and Machta et al. [120] found jj. - v + 2 
for the diffusion constant of ballistic tracers. 

The structural aspects of the Lorentz model with overlap- 
ping obstacles have been studied extensively in the context of 
continuum percolation. It shares the phenomenology of lattice 
percolation, and simulation results are consistently described 
by the same universality class [79, 106], even for obstacle 
mixtures with two different radii [122, 123]. In the light of 
renormaUsation group theory, this is no surprise and corrob- 
orates the mapping to random-resistor networks. Percolation 
thresholds were obtained with high precision yielding a criti- 
cal void porosity ipc = 0.0301(3) for spheres [106, 122, 123] 
and ifc - 0.323 652 5(6) for discs [124, 125]; numbers in 
parentheses indicate the uncertainty in the last digit. Simu- 
lations for the Lorentz model [22, 79, 115, 118, 126, 127] 
confirm the picture that the localisation transition is indeed 
driven by the percolation transition of the medium. For point 
tracers, the critical density of the localisation transition is de- 
fined by the percolation threshold. The critical porosity for 
localisation of a tracer with finite radius cr, follows trivially in 
the case of hard obstacles as (pdcr,) - exp(-n*(cr, -i- a-y'v), 
e.g., the localisation transition occurs at a void porosity of 
fcict - cr) - (fc^l^ ~ 65% for tracer and obstacles being 
spheres of the same size. A similarly high sensitivity of the 
percolation threshold to the tracer radius was found numeri- 
cally for obstacles modelled by soft repulsive discs [128]. 
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3.3.3. The ant in the labyrinth 

Consider a random walker ("the ant") on percolation clus- 
ters that has to explore the ramified and self-similar structure 
of the clusters ("the labyrinth"). The problem was posed by 
de Gennes [129] and is amenable to scaling arguments [129- 
132] corroborated by Monte-Carlo simulations [133]. If the 
walker is restricted to the incipient infinite cluster at critical- 
ity, one expects subdiffusive motion, 5r^(f) ~ f^^"'* , at all time 
scales beyond the microscopic regime, f » r := ct^/Dq. The 
exponent d^/ is known as walk dimension and one may prefer 
to consider it the fundamental critical exponent for the dynam- 
ics rather than the conductivity exponent, fi. Off criticality, 
where the infinite cluster becomes homogeneous on distances 
much larger than ^, tracer transport crosses over to normal dif- 
fusion, i5r^(f) - IdDcat, at long times, t » t^. The crossover 
time scale may be defined via 6rl^(tx) - and diverges as 
f.r ~ ■ Crossover matching yields t^J'^^ ~ D^^t^, which im- 
plies that the macroscopic diffusion constant vanishes at the 
critical density in a singular way, Dac(n*) ~ \n* - n*|'^", where 
jt/oo = v(dv, - 2). Transport on a finite cluster is subdiffusive as 
well with the same exponent 2/d^ as long as the tracer has not 
fully explored the cluster [116, 132]. 

If tracers on all clusters are included, the dynamics is non- 
ergodic since the time average over a single trajectory differs 
from an ensemble average. The average over a self-similar 
hierarchy of differently sized clusters reduces the exponent 
of subdiffusion, 6r^(t) ~ f^^' for t <s f <k t^, with the dy- 
namic exponent z given by z = dvi/[l - (d - df)/2] > d^,. 
In the absence of a percolating cluster, the mean-square dis- 
placement saturates for long times, 6r^(t) ^ for f » f^, 
and measures the mean-cluster size, { ~ t]J~ ~ ^■"(''"''r)/2^ 
which is distinct from the size of the largest finite clusters, ^. 
At the percolating side of the transition, the mean-square dis- 
placements for an all-cluster average increase linearly for long 
times, 6r^(t) ^ 2d Dt for f » f We shall refer to this be- 
haviour also as heterogeneous diffusion although only the mo- 
tion on the infinite cluster is diffusive, the all-cluster averaged 
propagator being different from a gaussian even for long times. 
Since the finite clusters do not contribute, the diffusion con- 
stant is further suppressed by the small weight of the infinite 
cluster yielding D(n*) ~ \n* - n*f with fi - v{d^ - 2 + d -df) 
for the all-cluster-averaged motion. By virtue of an Einstein 
relation, I. ~ D, the exponent /u corresponds to the conduc- 
tivity exponent. The previous relation connecting fi and d„ 
emphasises that the leading singularities of transport-related 
observables can be described by a single exponent, e.g., the 
walk dimension, d^,. 

The dynamic scaling hypothesis [134] suggests that the full 
time-dependence in the scaling limit is encoded in universal 
scaUng functions, which extends the notion of universality for 
the critical exponents. More specifically, a scaling form of the 
mean-square displacement is expected to hold at sufficiently 
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Figure 4. Ensemble-averaged mean-square displacements for differ- 
ent obstacle densities in the three-dimensional Lorentz model with a 
ballistic point tracer. Data taken from Ref. 115. 
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Figure 5. Time-dependent, ensemble-averaged diffusion coefficients 
for different obstacle densities in the two-dimensional Lorentz model 
with a Brownian point tracer. Data taken from Ref. 22. 



long times, f » r, 

6r^(t;n) = t^''6'Rl(i), 



i := f/ff ~ t{ 



(72) 



dli^ and (5!Rj denote scaling functions which describe the 
crossovers from subdiffusion to heterogeneous diffusion or lo- 
calisation, respectively. They approach the same constant for 
small arguments, and the correct long-time behaviour of 6r^(t) 
is reproduced by 6'Rl{t) ~ f'^'^ and 6'Rlit) ~ r^'^ for f ^ oo. 
Numerical results for these functions can be found in Refs. 22 
and 116 for d - 2 and in Ref. 115 for the Lorentz model in 
d = 3. The crossover from intermediate subdiffusion to the 
long-time behaviour occurs slowly and extends over several 
decades in time. 
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3.3.4. Simulation results 

The ensemble-averaged mean-square displacements of a bal- 
listic tracer in the three-dimensional Lorentz model were de- 
termined by computer simulations over the full range of obsta- 
cle densities [115], reproduced in Fig. 4. A temporal window 
of subdiffusive motion emerges as the critical obstacle density 
is approached, and at the critical density, subdiffusion is ob- 
served over 6 decades in time. The subdiffusion exponent has 
been determined as a = 2/z =s 0.32 and is consistent with the 
predictions for continuum percolation, yu - v + 2 x 2.88 [120]. 
Likewise, a tracer on the incipient infinite cluster obeys sub- 
diffusion with 2/t/w ~ 0.42 [127]. Away from criticality, the 
subdiffusive motion crosses over to either heterogeneous dif- 
fusion or localisation beyond the time scale f j. If particles on 
all clusters are considered, the non-gaussian parameter slowly 
diverges in time, aiCO ~ t^d-d,)/d^.^ reflecting the distinct length 
scales ( and ^ [79, 115]; it approaches a constant for tracers 
only on the infinite cluster [127]. 

The same phenomenology was found within simulations 
for the two-dimensional Lorentz model with a Brownian point 
tracer [22] and earlier for an obstructed random walker on 
a percolating square lattice [87]. The time-dependent diffu- 
sion coefficient, D(t) = d,6r^{t)/4, in Fig. 5 nicely exhibits a 
double-crossover scenario from free diffusion at microscopic 
time scales to a growing subdiffusive window, r <K f « f^, 
and back to heterogeneous diffusion with a suppressed diffu- 
sion constant for long times. 



D(f)- 



4Dof, f<KT, 

4Dt , t » tx; 



(73) 



the data are consistent with the exponent a - 2/z ~ 0.659 
from lattice percolation in two dimensions. 

In the context of ion conductors, the electric conductivity 
and susceptibility are relevant material parameters. Simulation 
data for the Lorentz model reveal an anomalous frequency- 
dependence inherited from the time-dependent diffusion coef- 
ficient, D{t), by a Fourier-Laplace transform [126, 127, 132]. 

Closer inspection of Fig. 4 suggests an apparent density 
dependence of the subdiffusive motion, which contrasts the 
notion of one universal exponent governing the subdiffusive 
regime and which seems to violate scaling, Eq. (72). The dis- 
crepancy is resolved by noting that the scaling property only 
holds asymptotically and that correction terms arise if the time 
scales under considerations are not large enough. It was ar- 
gued that the leading correction term is non-analytic and well 
approximated by a power law, 6r^{t; n*) = t'^^~S'Rl{t){l + Ct'-'), 
where C is a non-universal constant [115]. The correction 
exponent y, on the other hand, is universal and related to 
the correction of the cluster-size distribution at criticality. 



ing theory for transport on percolation clusters [116]. If cor- 
rections of this form are taken into account, the scaling prop- 
erty, Eq. (72), is corroborated by extensive simulation data for 
two- and three-dimensional Lorentz models [22, 79, 1 15, 127]. 
Both, significant corrections to scaling and slow crossovers 
imply that a power-law fit to experimental data over a limited 
time window is likely to produce apparent subdiffusion expo- 
nents that deviate from the universal value. Further, deviations 
from the idealisation of hard, immobile obstacles potentially 
wash out the sharp localisation transition. 

The mean-square displacement obtained in the computer 
simulations are ensemble averages where the initial position 
of the tracers have been chosen according to their equilib- 
rium distriubtion, i.e., anywhere in the void space. Further- 
more, an average over different realisations of the disordered 
environment is performed to ensure that all sample-to-sample 
fluctuations are averaged out. Nevertheless, it would be inter- 
esting to investigate the convergence of the averaging proce- 
dure. If the tracer is confined to the infinite cluster only, one 
should measure how long single trajectories have to be in or- 
der that the time-averaged mean-square displacement reflects 
the ensemble-averaged counterpart. 

3.3.5. Spatially resolved transport 

The scaling form of the propagator is more involved than for 
the models discussed so far due to the presence of finite clus- 
ters and the finite correlation length. For transport restricted to 
the infinite cluster, one expects [116, 135] 



P„(r,f;«*) = r"'!Po<,(r/^,f^-''») 



(74) 



P{s-n*) 



-\-d/d, 



in the scaling limit, r, ^ » cr and f » r; the density depen- 
dence is encoded on the rh.s. via the divergent correlation 
length, ^ ~ \n* - n'T". At the critical density, this simplifies 
to the one-parameter scaling form familiar from the preceding 
sections, 

P^rJ-nD^r-^'p^irt-'"-). (75) 

We have determined the scaling function !Poo( ) for a random 
walker on the incipient infinite cluster of a square lattice (site 
percolation) using Monte-Carlo simulations, see Fig. 2. For 
smalljirguments, the fractal dimension enters the scaling func- 
tion, Poo(x) ~ x''' for X <K 1, with consequences for the return 
probability, which is effectively probed in fluorescence cor- 
relation spectroscopy [23]. The scaling of the intermediate 
scattering function was studied numerically [136] and it was 
found that certain aspects are surprisingly well predicted by a 
mode-coupling approach. 

For transport phenomena including all clusters, a dynamic 
scaling ansatz for the van Hove function has been conjec- 
tured [115, 135], 



(1 + Cs ) for 5 — > oo. The exponent re- 



lation ydv, - ndf was derived within a cluster-resolved scal- 



P(r,t-n*)^^''-^'p^ir/^,t^-'^); 



(76) 
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the modified prefactor follows from averaging over the clus- 
ter size distribution [116] and the index ± indicates either 
the localised phase (+) or the regime of heterogeneous dif- 
fusion (-). The scaling form has implications for the mean- 
square displacement, Eq. (72), and the non-gaussian param- 
eter, both of which were corroborated by extensive simula- 
tions [22, 79, 115, 133]. 

A fraction of particles is generically trapped or "immobile" 
at long time scales by the presence of finite clusters at all ob- 
stacle densities. Immobile particles become manifest for all n* 
in a finite Lamb-MoBbauer factor or non-ergodicity parame- 
ter [135], 

f(k\n ):^\im P{k,t)>0, (11) 

in agreement with simulations [126, 137]. The long- 
wavelength limit, limt^o f{k), quantifies the immobile fraction 
and displays a singularity at the percolation transition [135]; 
more generally, a scaling form has been suggested [126], 

\ - f{k-n*) ^ ^''-'T^ik^) (78) 

for small wavenumbers, k 2n/cr, and close to the transition, 
^ » cr; the asymptotics of the scaling functions are f_{x — > 
0) = const, T+{x ^ 0) = 0, and T±{x) ~ x''^* for x » 1. 

3.3.6. Velocity autocorrelation function 

The obstacles impose an excluded volume on the tracer parti- 
cles, which forces the latter to reverse the direction of motion 
at some point. As a result, the particle displacements carry 
persistent anticorrelations which become manifest in an alge- 
braic long-time tail of the VACF, -Z(f) ~ f^''^^"', first derived 
within kinetic theory for low obstacle density [89]. More gen- 
erally, the power-law decay at long times emerges for diffu- 
sion in a frozen, disordered environment [138, 139], and it is 
akeady produced by repeated encounters of a Brownian tracer 
with a single obstacle [80], or equivalently, two colloidal par- 
ticles [140]. The phenomenon is universal and occurs inde- 
pendently of the obstacle density with the same exponent at 
sufficiently long time scales; close to the localisation transi- 
tion, the long-time tail competes with an intermediate, critical 
power law inherited from the subdiffusive motion [92] . These 
predictions were confirmed by simulations on a lattice [141] 
and in the continuum [22, 118]. Molecular dynamics simu- 
lations for supercooled mixtures of hard spheres [142, 143] 
and Lennard- Jones particles [144], both exhibiting slow glassy 
dynamics, reveal the same power-law signature in the VACF, 
which indicates that the underlying mechanism may also ap- 
ply for glassforming materials. A related long-time tail was 
predicted within a mode-coupling theory for the force-force 
correlation function in dense colloidal suspensions [145]. 

In the localised regime, n* > «*, where diffusive trans- 
port is absent, Machta [146] argued that another long-time tail 
becomes dominant if there are cul-de-sacs with a broad dis- 
tribution of exit rates. The singular distribution of channel 



width in the overlapping Lorentz model implies a long-time 
tail, -Z(f) ~ rl2+i/(d-i)]^ fQj. ^ ballistic tracer [1 14], which has 
been detected in recent simulations [118]. 

3.4. Other sources of subdiffusion 

To conclude the section, let us mention further examples of 
mechanisms, which are known to lead to subdiffusion. The 
review by Bouchaud and Georges [49] introduces a series of 
lattice models with quenched disorder, for example random 
trap and random barrier models or comb-like structures, which 
under certain conditions yield subdiffusive transport. Next, 
single-file diffusion describes the motion of strongly interact- 
ing particles aligned in a tube-like structure such that the ex- 
cluded volume impedes passing. As a consequence, the mean- 
square displacement grows as 6r^{t) ~ f'^^ [147-149]. 

For flexible polymers, the chain connectivity induces 
strongly correlated motion between the individual chain seg- 
ments. Transport slows down as the polymer weight in- 
creases. For intermediate time scales, a labelled monomer dis- 
plays a subdiffusive mean-square displacement until the slow- 
est Rouse mode has relaxed [150-153]. Self-avoidance due to 
excluded volume and entanglement effects from topological 
constraints add to the complexity of the dynamics. The theo- 
retical implications from the interplay of polymer physics and 
crowding (excluded volume) remain largely unexplored to this 
date. 

Apparent subdiffusion may result from incorporating sev- 
eral processes that occur on different time scales, for exam- 
ple in multi-component mixtures, due to a polymeric deple- 
tion layer around tracers [154, 155], or due to internal states. 
In these cases, the processes are intrinsic to the constituents or 
their interactions and hence cover finite time windows only. In 
particular, one can usually not manipulate these processes to 
generate a self-similar distribution of time scales required for 
anomalous transport in the sense of Section 2.2. 

4. Experimental techniques 

Before we will review the biophysical experiments addressing 
anomalous transport, let us introduce some widespread exper- 
imental techniques that have been developed during the past 
decades and that have proved themselves as useful tools for 
the measurement of molecular transport in mesoscopic sam- 
ples. 

4.1. Single-particle tracking 

An important technique that is both intuitive and powerful is 
provided by single-particle tracking, for reviews in the con- 
text of biophysics see Refs. 19 (giving an extensive historical 
account in the supplement) and 156. Here, a nanoscopic re- 
porter particle or a fluorescent dye is introduced in the probe 
and followed by confocal video microscopy and digital image 
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processing. The trajectory is recorded over sufficiently long 
time, which provides direct access to the full statistics of the 
spatial displacements and which, in principle, allows for the 
evaluation of all correlation functions. The most widely used 
quantity is, of course, the mean-square displacement as it is 
rather robust with respect to experimental noise and appears 
easy to interpret. 

The spatial resolution of the displacements is typically on 
the order of a few nanometres, while the temporal resolution is 
limited by the image capture rate to about 10 ms [156]. High- 
speed tracking with a resolution of 25 ps could be implemented 
by using colloidal gold nanoparticles as tracers, which yielded 
a sufficiently high signal-to-noise ratio [19]. The length of 
the recorded trajectories and thus the longest accessible time 
scale is limited by photobleaching of the fluorophores and by 
residual drifts of the sample baseline. Single-particle tracking 
is well suited to study transport in cellular membranes and to 
elucidate their structural details [19, 157]. Limitations occur 
for tracking inside the cytoplasm since the spatial resolution 
transverse to the focal plane is naturally an order of magni- 
tude lower [158]; progress was made recently using tailored 
computer algorithms that address the noise and image inho- 
mogeneities specific to the cytoplasm [159]. Single-particle 
tracking has also become a widespread tool to probe the mi- 
crorheology of the cytoplasm, see Refs. 160-162 for recent 
reviews and references therein. 

Three-dimensional tracking at the nanoscale has seen sig- 
nificant advances recently. Truly three-dimensional trajecto- 
ries in real-time at a resolution of 32 ms are provided by or- 
bital tracking of a single particle, which relies on a feed-back 
loop coupled to the laser-scanning microscope [163]. The 
"blindness" of this technique to the neighbourhood is over- 
come by simultaneous wide-field imaging, which is partic- 
ularly useful to track particles with nanometre resolution in 
a heterogeneous environment [164]. Off-focus imaging of 
trapped particles provides even subnanometre resolution in 
all three directions; it requires, however, information about 
the optical properties of the sample [165]. An alternative 
to strongly trapped probes is interferometric detection, where 
the motion can be resolved with sub-Angstr0m resolution at 
75 MHz [166], yielding significant time-correlation functions 
over 4 decades in time [40, 46] . 

Central for the analysis of single-particle tracking experi- 
ments is the evaluation of time-correlation functions as they 
contain the proper statistical description of the motion, see 
Section 2.2. The mere inspection of individual stochastic tra- 
jectories may easily lead to misinterpretations [157, 167]. The 
discussion and computation of time-correlation functions and 
related quantities is simplified in the case of time-translational 
invariance, i.e., for samples in equilibrium or in a steady state. 
Then, the correlation functions depend on the lag time only 
and provide information about the motion at a certain time 
scale (and not at a certain point in time); they are naturally 



displayed on a logarithmic time axis. For example, the mean- 
square displacement of a single trajectory x{t) of finite length 
T given on an equidistant time grid is efficiently evaluated as 
a time average employing a fast (discrete) Fourier transforma- 
tion T, 



C„(f; T) := (x(f) ^(0))^ = ^ f x(t' + f) x(t')At' 

i Jo 



\TAx(t)]\ 
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where oj,, = Inn/T and 
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x(t) e"'"' dt ; 



(82) 



the mean-square displacement then follows as 6x^{t) - 
2[C„(0; T) - CAf, T)] for t <k T. For very long trajecto- 
ries, a semi-logarithmic sampling of the trajectory data and a 
direct evaluation of the correlations has proved a useful alter- 
native [168]. 

Automated experimental assays can measure several thou- 
sand trajectories and thus permit the evaluation of displace- 
ment histograms, i.e., the propagator P(r, f), or of the squared 
displacements, see Section 2.2.2. Analysis of these distribu- 
tions helps to disentangle different transport processes, for ex- 
ample in mixtures of slow and fast particles [169, 170] or when 
stationary and non-stationary processes coexist [171]. Second, 
the propagator gives insight into the spatial aspects of a trans- 
port process that would remain undisclosed by a mere investi- 
gation of the mean-square displacement, as has been discussed 
extensively in Section 3. For example, tracking experiments 
demonstrated that transport characterised by a diffusive mean- 
square displacement may be non-gaussian in space [172]. Fi- 
nally, the ensemble of time-averaged mean-square displace- 
ments may display large fluctuations, in particular if transport 
is anomalous. Comparison of ensemble and time averages 
tests the ergodic hypothesis and may shed light on the under- 
lying processes in nanoporous structures [173, 174] as well as 
in cellular membranes [171]. 

4.2. Fluorescence correlation spectroscopy (FCS) 

A second powerful technique that has been widely applied 
in the biophysical context is fluorescence correlation spec- 
troscopy (FCS). Here the basic idea is to label few molecules 
by a fluorescent dye and record the fluctuating fluorescent light 
upon illuminating a small part of the sample. We briefly in- 
troduce the theory underlying the measurement and discuss 
how anomalous transport manifests itself in the corresponding 
correlation function, comprehensive information about the ad- 
vantages and hmitations of FCS may be found in the pertinent 
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literature, for reviews see Refs. 175-179. Compared to single- 
particle tracking, an advantage of FCS is the high temporal 
resolution in the microsecond regime and the use of tracers 
as small as a single fluorophore. On the other hand, the spatial 
sensitivity is typically on the order of some 100 nm (exceptions 
are discussed below), orders of magnitude larger than that of 
high-precision tracking experiments. 

A typical FCS setup consists of an illumination laser and 
a confocal microscope with a photon detector.^ The laser 
beam illuminates the detection volume with the intensity pro- 
file W(r) and excites the fluorophores in the focal volume. 
The emitted fluorescent light is collected in the detector, it 
depends on the fluctuating, local concentration c(r, t) of la- 
belled molecules. Thus, the detected light intensity is a spa- 
tially weighted average [Ref. 14, section 6.6], 



I(t) 



r W{r) c(r, f) , 



(83) 



where the prefactor e accounts for the total quantum efficiency 
of absorption, fluorescence, and detection. The output of the 
FCS experiment is the time-autocorrelation function Git) of 
the intensity fluctuation 51{t) - I(t) - (I) around the mean 
intensity; it is conventionally normaUsed as 



G{t) ^{61(1)61(0)} I {if . 

Introducing a spatial Fourier transform, W(k) 
Jd'^r e'*'"'' W(r) , and the intermediate scattering function. 



(84) 



d''re"'''-'<&(r,f)5c(0,0)). 



one arrives at the expression 



G(t) 



1 Jd'^k\W(li)\^S(li,t) 
N /d''fe|W(k)p 



(85) 



(86) 



where :- {if / {si^") - {c)Vef{ is interpreted as the number 
of fluorophores in the effective illumination volume Veff. In 
the experiments, the fluorophores are highly diluted and V^n 
usually contains only few molecules or even less than 1, turn- 
ing FCS essentially into a single-molecule fluorescence tech- 
nique [178]. Thus, 5(k, f) reduces to the incoherent interme- 
diate scattering function. 



5(k,f)^P(k,f) = (e""'"^'*^'^) 



(87) 



Inserting this in Eq. (86) and interchanging the k-integration 
and the statistical average provides the fundamental expression 
for the FCS correlation function [23], 

G(f) '^i^^^''^ l^(k)l^ (88) 

^ In practice, the detector consists of two cross-coiTelated avalanclie photodi- 
odes to reduce tlie detector noise. 



It depends solely on the approximation of dilute labelling, it 
holds for arbitrary illumination profile and makes no assump- 
tions on the fluorophore transport. 

For practical purposes, however, Eq. (88) shall be evalu- 
ated further. A conventional laser emits a gaussian beam pro- 
file, which together with the usual confocal setup leads to an 
illumination profile that is often approximated by an elongated 
ellipsoid. 



W{r) oc exp (-2(x2 + y^)lw^ - 2z^/zl) 



(89) 



with beam waist w and longitudinal extension zo- 
This implies a gaussian filter function, |H^(k)| oc 
exp(^-{kl + k^. + rf-k^)w^ , introducing the anisotropy 
parameter rj - zo/w. Then, the k-integration can be carried 
out and gives 



. 1 / / AR(0 

G(r;w)=- exp -— ^ 



' + ^-^AZitf 



(90) 



For planar motion, the displacement along the beam axis van- 
ishes, AZ(f) - 0, and a compact statistical expression for the 
FCS correlation follows. It reveals the close relationship to the 
probability distribution of the squared displacements, Eq. (29), 
if G(f; w) is interpreted as the corresponding characteristic 
function. The similarity of the representations in Eqs. (87) 
and (90) suggests that the FCS correlation encodes important 
spatial information analogous to scattering methods as soon as 
the beam waist w is considered an experimentally adjustable 
parameter [23]. Several FCS setups with spatio-temporal res- 
olution have been implemented recently by introducing vari- 
able beam expanders [180, 181], z-scan FCS [182, 183], sub- 
wavelength apertures [184, 185], near-field scanning optical 
microscopy [186], or stimulated emission depletion (STED) 
permitting spot sizes as small as 20 nm [187-189]. Finally, 
Eq. (90) provides a simple scheme for the efficient evalua- 
tion of autocorrelated FCS data in computer simulations from 
a given trajectory. 

For the analysis of a specific FCS experiment, the statisti- 
cal average needs to be performed and some knowledge about 
the statistical nature of the fluorophore displacements is re- 
quired. A common, but in the context of anomalous transport 
strong assumption is that of a gaussian and isotropic distri- 
bution of the displacements AR(f) after a fixed time lag, i.e., 
{AR(0) - and only the second cumulant 6r'^(t) = ^|AR(f)P) 

IS non-zero. It foflows that f(k,f) = exp[-k'^dr'^{t)/2d) being 
the characteristic function of the random displacements, and 
the FCS correlation in c/ = 3 is given by [190] 



1 + 2^ 



1 + 2^ 



-1/2 



which simplifies in two dimensions to [191] 
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Both expressions allow for a (numerical'') inversion at fixed 
time t and conveniently provide direct access to the mean- 
square displacement from an FCS experiment [192, 193], 
without resorting to fitting algorithms. This will be partic- 
ularly useful in the generic situation where the transport ex- 
hibits different regimes depending on the time scale, including 
crossovers from subdiffusion to normal diffusion (see, e.g., 
Ref 19). The quality of the gaussian approximation can be 
tested experimentally by resolving the spatio-temporal proper- 
ties of the molecular motion [23], e.g., by systematic variation 
of the confocal volume [180-189]. Then, the mean-square dis- 
placements obtained from Eqs. (91) for different values of w 
should coincide if the gaussian approximation is applicable. 

Specialising to free diffusion, it holds P(k, t) - 
exp (^-Dk^t^, and G(t) attains the well-known form 



Gfree(f) = ^(l + t/TDy\l + TJ-^t/ToJ 



1/2 



(92) 



with the dwell time td = w^/4D. Imposing subdiffusion at 
all times, 6r^(t) ~ f , and assuming gaussian spatial displace- 
ments as in fractional Brownian motion, Eqs. (91) yield 



^FBM 



(0 



^ [l + (tlTofX' [l + il-\tlTof\'^\ (93) 



This form has been widely employed in FCS-based studies 
of anomalous transport, maybe because it simply supplements 
the usual model, Eq. (92), by an additional fitting parameter a. 
Care must be taken that the additional fitting parameter a is not 
abused for disclosing experimental deficiencies (e.g., drift) or 
apparent subdiffusion where the physical picture would actu- 
ally favour several species with different diffusion coefficients. 

For FCS experiments with variable detection volume, it is 
beneficial to develop some scaling properties of the FCS corre- 
lation [23]. Starting from a fairly general scaling ansatz for the 
transport propagator, P{r,t) - r^'''P{rr^^^^), Eq. (90) implies 
the existence of a scaling function such that 



G(f; w)/G(0;w) = Qitw'''^) . 



(94) 



A representation of the l.h.s. as function of the scaling variable 
p^-d^ should collapse the FCS correlation data for different 
beam waists w onto a single master curve [23]. In particular, 
the half- value times ri/2, defined via G(ti/2)/G(0) = 1/2, are 
expected to obey 



Tl/2(w) ~ W 



(95) 



Wawrezinieck et al. [181] suggested to study the phe- 
nomenological relation t^(w) = fo + w^/4Deff, termed FCS 

' In three dimensions, one may alternatively consider to expand the last fac- 
tor of Eq. (91a) for 2> 1, yielding the approximation 

^ , , \-NG(t) [, ^ 3 \ -NG(t) _4 



2NG(t)-r]-^- [ 8 [/VGWP 
which is useful as long as NG{t) > rf^. 



diffusion law, where Td(w) is obtained from fitting the sim- 
ple diffusion model to the FCS data for different beam waists. 
Based on simulations, they proposed that the transport mecha- 
nism can be inferred from the sign of the axis intercept Iq: neg- 
ative values would indicate transport hindered by barriers and 
?() > would hint at transport in the presence of microdomains 
with slower diffusion. For a meshwork model, where the tracer 
diffuses freely in square domains separated by barriers of finite 
probability, the relation ti/2(w) = (w^ - fl^/12)/4D has been 
derived analytically for large w, with a being the mesh size and 
D the long-time diffusion constant [194]. This is consistent 
with the empirical findings of Ref. 181. A thorough theoreti- 
cal analysis for other transport models, however, remains still 
to be done [195]. 

A general representation of a correlation function for pure 
relaxation dynamics, e.g., for the incoherent intermediate scat- 
tering function, is given by [196] 



P(k, t) 
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with a set of positive and A:-dependent probability distribu- 
tions Ilk for the relaxation times r. From Eqs. (87) and (88), 
one finds for the FCS correlation function 



G(t) 



/-*oo 

e'/^Hv 

Jo 



iT)dT, 



(97) 



where the distribution 11^ is given by 11^ and depends on the 
illumination profile, TIw(t) oc Jd''k\W(k)f'Ilk(T). An interpre- 
tation in terms of many diffusing components may be obtained 
by transformation to 



(98) 



G(0 = f Gfree(f;To)n^o(TB)dT£,, 

Jo 

using the representation 

(l+f/r^)-"/2= fe-'/- ^"°^"C/r°'^ dT. (99) 
Jo Td r(ii/2) 

The specialisation to two freely diffusing components is 
widely employed to fit experimental data for transport in 
crowded media. 



^ToiTo) = f\S(TD - Tl) -I- flSiTD - T2) . 



(100) 



For normal as well as for subdiffusive motion, G(f) dis- 
plays a power-law decay at long times, G(f) ~ r^, imply- 
ing a power law also for the distribution of relaxation times, 
niy(T) ~ T"^"' as T — > CO, by means of a Tauber theo- 
rem [53, 59]. Iterating the application of the Tauber theorem, 
the distribution of residence times displays a power-law tail. 



n.„(TD) 



for Td — > only in the case of subdiffu- 



sion, y6 < d/2. Note that for spatially non-gaussian transport, 
the exponent /3 and the subdiffusion exponent may not be con- 
nected by a simple relation; in case of the Lorentz model, it 
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involves the fractal space dimension, fi - dfjd^ [23]. The nu- 
merical determination of the distributions of relaxation times 
from a given data set for a slowly decaying correlation func- 
tion is a formidable task [196]; it amounts to performing an in- 
verse Laplace transform, which is a mathematically ill-posed 
problem [178, 197]. Imposing the additional constraint that 
the "entropy" of the distribution shall be maximal, Sengupta 
et al. [198] nevertheless introduced a numerical procedure to 
determine the distribution of diffusion times from a given 
FCS correlation function G{t), which was successfully applied 
to the interpretation of experiments [199, 200]. 

For the measurement of slow, anomalous transport, it is 
crucial that the FCS correlation covers several decades in time 
and that the data at long times are still significant. In par- 
ticular, the choice of an appropriate fit model requires suffi- 
cient knowledge about the physics of the sample, and physical 
constraints on fit parameters or fit windows need considera- 
tion [201]. The application of FCS in complex cellular envi- 
ronments has several limitations which may artificially induce 
weakly subdiffusive behaviour [202, 203] or modify the prop- 
erties of true anomalous transport [201]; some of these limita- 
tions are addressed by advanced derivatives of the FCS tech- 
nique [178, 204]. We agree that "the physical insight gained 
from the empirical application of the concept of anomalous 
diffusion to experimental FCS data may be quite limited, es- 
pecially if no convincing microscopic origin of the deviation 
from the normal diffusion law is provided" [178]. This con- 
clusion emphasises the need for more refined theoretical mod- 
els and for the possibility of their experimental discrimination. 
Long measurement times, a high signal-to-noise ratio, and the 
combination of spatial and temporal information seem essen- 
tial for this task. 

4.3. Fluorescence recovery after photobleaching (FRAP) 

Simultaneously to the FCS technique, fluorescence recovery 
after photobleaching (FRAP) was developed [205, 206] and 
readily applied to membrane proteins on cells [205, 207-209], 
see Refs. 179, 210-214 for reviews. The method is based 
on a similar experimental setup as FCS, but initially the flu- 
orophores in the observation region are bleached by a brief, 
intense laser pulse. Afterwards the fluorescent light intensity 
emitted from this region is monitored while it recovers due 
to the diffusion of unbleached fluorophores from outside into 
the observation region. The method is in some sense comple- 
mentary to FCS: it uses a high fluorophore concentration and 
measures the collective transport in form of a diffusion front. 
It is applicable to very slow processes (even on the scale of 
several seconds [214]) and appears more robust than FCS if a 
significant fraction of molecules is immobile [179]. FRAP has 
become a valuable tool to measure protein motion and activity 
in drug delivery in pharmaceutical research [211] as well as in 
living cells [213]. 

A transparent presentation of the theoretical background 



was given by Elson [210] for normal diffusion in two dimen- 
sions, which we will briefly summarise. The experiment mea- 
sures the response of the local fluorophore concentration to 
an initial quench out of equilibrium. Similar to FCS, the de- 
tected fluorescent light intensity I(t) is given by Eq. (83) as 
integral over the time-dependent, local concentration c(r, f) 
weighted with the intensity profile of the detection laser W(r). 
Let Ceq(r) s Co be the homogeneous background concentration 
of fluorophores before bleaching, c(r, 0) the imprinted profile 
directly after the bleaching pulse, and <5c(r, f) :- c(r, f) - cq 
the partially recovered disturbance after a time t. Repeated 
measurements of the intensity evolution /(f) result in a non- 
equilibrium average, {/(O)ne' fluctuations across different 
measurements are not considered. The prebleach intensity and 
the recovery curve are then given by 

:= </(f ^ cx,))^^ = SCO J d''r W(r) , (101) 

- <'f(0>„.e. = /-co + e J d'r W(r) {6c(r, ?)>„,. ■ (102) 

The postbleach intensity is defined as /q lim,^o/,, and we 
will give results for the reduced fluorescence recovery function 
in the following, 

I(t) := . (103) 

Assuming an irreversible conversion of the bleached 
molecules to a non-fluorescent state, the fluorescence recovery 
is solely governed by fluorophore transport. Arguing along the 
lines of Onsager's regression hypothesis, the macroscopic re- 
laxation to the equilibrium distribution follows the same laws 
as the regression of a microscopic disturbance of the local flu- 
orophore concentration induced by thermal fluctuations [Ref. 
11, section 7.6]. (While the fluctuation-dissipation theorem 
supports this procedure for small concentration gradients, even 
for complex dynamics, the validity of this assumption beyond 
linear response remains to be corroborated.) Specifically, the 
time evolution of a disturbance 6c{r, 0) is governed within lin- 
ear response theory by the equilibrium correlation function of 
the concentration fluctuations, namely the van Hove function, 

(5c(r, f))„.e. = J"dV5(r-r',f)(5c(r',0), (104) 

where 

5(r-r',f) = — {6c(r,t)6c(r',0)) (105) 
Co 

for spatially and temporally homogeneous transport. In 
Fourier space, perturbations with different wavevector decou- 
ple and their relaxation is dictated by the intermediate scatter- 
ing function, Eq. (85), 

{6c(K f)>„.e. = 5(k, f) (5c(k, 0) . (106) 
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For the recovery of detected fluorescence intensity, one finds 



J d'^rdV W(r)S(r - r',t)Sc(r',0) (107) 



= r ^W(-k)6c(k,0)S(k,t) 



(108) 



The reduced fluorescence recovery function foflows as 

J(f) = Jd''k^V(k)S(k,t), (109) 



introducing a filter function. 
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W(-k)6c(k,0) 
Jd''kW(-k)6c(K0) 



(110) 



to describe the specific experimental setup. Note that bleach- 
ing and detection enter independently and may be imple- 
mented with different beam profiles. 5(k, f) is solely deter- 
mined by the physics of fluorophore transport in the investi- 
gated sample. If the molecules diffuse freely. 



5free(k,0 = exp(-Dcrf), 



(111) 



where the collective diffusion constant Dc differs in general 
from the diffusion constant of a single tracer molecule. 

Modelling the fluorophore bleaching as an irreversible 
first-order reaction, the concentration of unbleached fluo- 
rophore immediately after the bleaching pulse is given by 
c(r, 0) = CQe^'^^'^\ where K(r) is proportional to the beam pro- 
file of the bleaching laser and to the duration of the pulse. For 
free diffusion in two dimensions, bleaching with a gaussian 
beam, Eq. (89), and detection with the same, but attenuated 
beam, the recovery of fluorescence intensity can be expressed 
in terms of the incomplete gamma function or as an infinite 
series [206], 
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1 -Hn(l +2t/Td) 
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it involves the bleaching parameter K = K(0) and the diffusion 
time Td = w^/4Dc. The postbleach intensity is calculated as 
/o = /-oo(l - e-'^)/K. 

If the beam profiles of both the bleaching and the detec- 
tion laser are approximated by a step-like disc of radius w, the 
filter function ^(k) assumes a relatively simple form and is 
independent of the bleaching parameter. 



^(k) 



Jiikwf 
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where J\{-) denotes the Bessel function of the first kind; the 
postbleach intensity is Iq - /_ooe"^. In the case of free dif- 
fusion, Eqs. (109) and (ill) yield the fluorescence recovery 
curve after performing the k-integration [215], 



Jfree(f) = 1 - c'^' ' ' ' [hilT d 1 1) + /l(2Trf/f)] ; 
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/o( ) and /i( ) denote modified Bessel functions of the first 
kind. 

Bleaching patterns other than a spot were implemented ex- 
perimentally. In fluorescence pattern photobleaching recov- 
ery (FPPR), a fringe pattern of periodic stripes is bleached to 
probe the anisotropic dynamics at finite, non-zero wavenum- 
ber k [216-219]. Using a laser-scanning confocal microscope 
allows for bleaching and probing an arbitrary geometry, e.g., 
a line segment, with diffraction-limited resolution [220]. Fur- 
thermore, the size of a circular bleaching spot can easily be 
varied, which permits the intrinsic determination of the instru- 
mental resolution parameters [221] and which can elucidate 
the spatio-temporal properties of anomalous transport analo- 
gous to FCS experiments with variable observation volume. 

In some experiments, the initial fluorescence is only par- 
tially recovered, even after long measurement time, which is 
commonly attributed to further photobleaching, chemical re- 
actions, or immobile fluorophores [211]. The unrecovered flu- 
orescence intensity is quantified by the ratio 



R 



I-co — Iq 
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^). 
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In this case, fluorophore transport is, strictly speaking, not er- 
godic and one infers a finite long-time limit of S (k, t) from 
Eq. (109), 

F(k) := lim5(k,f), (116) 



determining the incomplete recovery. 



R 



/feW(k)F(k). 



(117) 



In the context of glassy dynamics, F(k.) is known as non- 
ergodicity parameter or Debye-Waller factor [21], and it plays 
the same role as the Lamb-MoBbauer factor in the case of the 
incoherent dynamics, c.f. Eq. (77). It follows that the ratio 
R for the same sample depends via W(k) on the geometries 
of the bleaching and detection spots. In particular, FRAP ex- 
periments with variable beam waist w of the laser probe the 
non-ergodicity at different scales. Salome et al. [222] sug- 
gested that R is an affine function of the inverse radius, 1/w, 
which is motivated by simulations and experiments on com- 
partmentalised membranes. Let us assume that a fraction Ri of 
the fluorophores is immobile and homogeneously distributed 
and that the mobile part is ergodic despite the presence of im- 
mobile molecules. Then, f (k) = Ri and the unrecovered in- 
tensity is identified with the immobile fraction, R - which 
is regularly encountered in the literature. More generally, Rj, 
is obtained as small-wavenumber limit of the non-ergodicity 
parameter, 

Ri ^ lira F(k). (118) 

For the study of anomalous transport, Feder et al. [223] 
generahsed the FRAP recovery curve, Eq. (112), by assuming 
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a gaussian subdiffusion model, which effectively amounts to 
replacing f/r^ by (tlTd)" ■ For the intermediate scattering func- 
tion, this implies a gaussian ansatz in space and a stretched 
exponential decay in time, 

5(k, f) = exp(-(f /T,n, (X r^/^ (119) 

Such a stretched fluorescence recovery occurs quite slowly and 
long measurements are required to distinguish the subdiffu- 
sion model from incomplete recovery [223]. More general 
considerations on FRAP for anomalous transport are still to 
be worked out; in particular, an extension to spatially non- 
gaussian transport would be useful. To this end, Eq. (109) may 
serve as starting point for a similar scaling analysis as above 
for the FCS correlation function. 

5. Anomalous transport in crowded biological media 

Application of the discussed techniques to crowded biologi- 
cal media have led to a plethora of experimental results on 
anomalous transport in the cell interior and in related model 
systems. In the main part of this section, we provide a compi- 
lation of the most significant biophysical experiments address- 
ing anomalous transport and report on the progress made dur- 
ing the past decade. Since a proper characterisation of anoma- 
lous transport requires many time and length scales, emphasis 
is put on the scales investigated by the different experiments. 
The presentation is divided into three-dimensional transport 
in cellular fluids (Table II) and crowded model solutions (Ta- 
ble III), and quasi-two-dimensional transport in cellular mem- 
branes (Table V) and Upid bilayer models (Table VI). Each 
experimental subsection is supplemented by a discussion of 
related computer simulations (Tables IV and VII). The section 
closes by briefly accounting for recent insights in the implica- 
tions of crowding for biochemical reactions. 

The topic is rather controversially discussed in the bio- 
physics community. Shortly after the turn of the millennium, it 
was recognised that the heterogeneous structure of intracellu- 
lar environments can not be neglected in the modelling and de- 
scription of macromolecular transport and biochemical reac- 
tions, with impact on basically all intracellular processes [224- 
227]. While a number of sound experiments on intracellu- 
lar transport and related model systems have observed pro- 
nounced subdiffusion, other experimental findings are well ex- 
plained in terms of normal diffusion [158, 228]. The situa- 
tion becomes even more delicate since "the sub-optimal ex- 
perimental conditions often encountered in cellular measure- 
ments do not allow ruling out simple Brownian diffusion mod- 
els" [228]. Further, the microscopic mechanisms underlying 
the subdiffusive motion have not been identified unambigu- 
ously until today. We agree with Elcock [229] that there is 
clearly a need for theory and simulation of microscopic mod- 
els that can make qualitative and quantitative predictions of the 
transport behaviour in crowded environments, at least in vitro. 



5.1. Crowded cellular fluids 

5.1.1. Cytoplasm and nucleoplasm (in vivoj 

Measurement of molecular transport in living cells faces a se- 
ries of complications typical for in vivo experiments: the cell 
size of 1 to lOOpm puts a natural upper Umit on the accessi- 
ble length scales, cells may be in different internal states of 
their cell cycle, and artificial probes may trigger a specific cell 
response. Further, transport may depend on various cellular 
processes like directed and active motion by motor proteins 
along the microtubule network, or cytoplasmic flows induced 
by the cytoskeleton pushing and pulling organelles around and 
locally liquefying the cytoplasm [251-253]. For example, the 
transfection pathway of viruses or gene carriers, visualised by 
in vivo single-particle tracking, crucially involves an interme- 
diate phase of slow passive transport, although cargo trans- 
port over large distances is mainly driven by molecular mo- 
tors [254, 255]. Such enzymatically driven processes may be 
suppressed or switched off by reduction of temperature or by 
depletion of ATP [256]. In the following, we will focus on 
passive transport induced by thermal fluctuations. 

An interesting question is to which extent may the cell be 
considered in a stationary state? The cell is certainly not in an 
equilibrium state, rather a plethora of biochemical processes 
is continuously taking place. The cell's life time is finite, it 
ages, and the cellular structures and processes change during 
the cell cycle. These changes, however, occur usually orders 
of magnitude slower than the typical experimental time scales 
of milliseconds or seconds at which intracellular transport is 
studied. We anticipate that a description of the cell as a sta- 
tionary state is a useful approximation, which of course has its 
limitations. The issue is partially addressed by some of the ex- 
periments covered in this section, a definite answer, however, 
requires future experimental work. 

The fission yeast Schizosaccharomyces pombe is a eukary- 
ote with a stiff cell wall that maintains a stable cylindrical 
shape, 12 [im in length and 4 |jm in width. It naturally contains 
lipid granules with diameters of about 300 nm, which can serve 
as endogenous, nearly spherical probe particles [230-233]. In 
their prominent work, Tolic-N0rrelykke et al. [230] combined 
the advantages of 2-3 min long trajectories from video mi- 
croscopy (at 25 Hz) and the high temporal resolution of laser- 
based particle tracking (22 kHz) to obtain mean-square dis- 
placements (MSDs) spanning a time window of more than 4 
decades. The data sets from the two techniques can be con- 
sistently extrapolated, bridging the gap between 1 and 40 ms. 
The data provide clear evidence of subdiffusive motion with 
an exponent between 0.70 and 0.74 in essentially the whole ac- 
cessible time window. The optical trapping experiments were 
refined five years later by the same laboratory [231], focus- 
ing on the different phases of the cell cycle. From the posi- 
tional power spectrum, the subdiffusive motion was confirmed 
at time scales between 0. 1 and 10 ms in all phases with a small. 
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cell type 


probe (size) 


experimental 
technique 


temporal and spatial 
scales investigated 


subdiffusion exponent 
or diffusion constant D 


year 


reference 


S. pombe 


lipid granules 
(300 nm) 


laser tracking, 
video microscopy 


50 |is-l ms, 
40 ms-3 min 


0.70-0.74 


2004 


ref. 230 






laser tracking 


0.1 ms— iUms (rt)U) 


0.81 (interphase) 
0.84 (cell division) 


2U09 


let. lil 






video microscopy 


lOms-lUs 


0.4 (time average) 


2010 


ref. 232 






video microscopy, 
laser tracking 


> 10 ms, 
0.1 ms-1 s 


0.6, 

normal for / < 3 ms, 
then y8 = 1 - a a 0.2 


on 1 1 


ret. Hi 


E. coli 


mRNA(lOOnm) 


video microscopy 


1 s-30 s, 

10^^-1 Hz (PSD) 


0.70 ± 0.07, 
0.77 ± 0.03 


2006 


ref. 234 






video microscopy 


I S— 10 S 


O./i + O.IU 


on 1 n 
2U1U 


let. liJ 


E. coli, 

C. crescentus 


chromosomal loci 
(GFP labelled) 


video microscopy 


1 S-IO^ S 


0.39 + 0.04 (ensemble and 
time averages coincide) 


2010 


ref. 235 


human ostero- 
sarcoma cells 
(nucleus) 


telomeres 
(GFP labelled) 


video microscopy 


lOms-lh 


subdiffusion with a varying 

with lag time, 

a\t). v.jZ — ^ u.jl — ^ 1. 


2009 


ref. 236 


SV-80 cells 


PS beads (3 |jm) 


video microscopy 


40ms-50s 


0.5-1 forf > 10 s, 
ill tor r < 3 s (motor 
proteins) 


2002 


ref. 237 


mammalian & 
plant cells 


rhodamine dye 


FCS 


0.1 ms < Ti/2 < 1 s 


0.6, or 2 components: 

Di =Daq/5,£)|/Z)2 =40 


1999 


ref. 238 


COS-7 and AT-1 
cells 


EGFP proteins 


FCS 


T|/2 « 1 ms 


0.7-1, or 2 components: 
£»! = A,q/5,£'i/£>2 = 10 


2000 


ref. 239 


HeLa cells 


FITC-dextran 
(1.8-14.4nm) 




0.4 ms < Tif2 < 16 ms 


U. / 1— U.o4 (depending on 
tracer size) 


zUU4 


let. Z4U 


HeLa and liver 
cells 


gold beads (5 nm) 


FCS 


10 [is-1 s, T|/2 ~ 0.3 ms 


0.53, with sucrose added: 
0.66 


2007 


ref. 241 


mammalian 
cells 


gold beads (5 nm) 


FCS 


0.1 ms < Ti/2 < 0.9 ms 


0.52 (cytoplasm), 
0.58 (nucleoplasm) 


2007 


ref. 242 


HeLa cells 


DNA 

(20bp-4.5kbp) 


FCS 


5 ms < Tfli < 20ms, 
To2 up to 500 ms (3 kbp) 


normal for small tracers, 
2 components above 250 bp: 
Di = Oaq/40 (size 
dependent) 


2005 


ref. 243 


D. discoideum 


GFP-actin, 
free GFP 


FCS 


170 |js (globular actin), 
240 |js (tree GFP) 


0.83 (both tracers) 


2010 


ref. 244 


3T3 fibroblasts 


dextran, FicoU 
n-58 nm"! 

-J O 11111 J 


FRAP 


50 [jm spot size 


normal with incomplete 

rppnvprv I~) c]f^ir\RT\f]^ on 

tracer size, D-^^/D up to 30 


1986, 
1987 


refs. 245, 
246 


3T3 fibroblasts 
MDCK cells 


dextrin Ficoll 
(4-30 nm) 


FRAP 


1 |js — 5 s 


normal, D-^q/ ^ 4, no 
dependence on tracer size 


1997 


ref. 247 


myotubes 


globular proteins 
(1. 3-7.2 nm) 


fringe pattern 
photobleaching 


Ti/2 ~ Is and 208 s, 
10-17 ^im wide stripes 


normal, D^q/D: 3-7, but 
10* for largest protein 


2000 


ref. 248 




FITC-dextran 
(2.9-12.6 nm) 


fringe pattern 
photobleaching 




normal, D^q/D: 2.5-7.7 


1996 


ref. 249 


E. coli 


GFP 


FRAP 


0.1-1 s 


normal, D a; Daq/ 1 1 


1999 


ref. 250 



Table II. Overview of in vivo experiments on crowded cellular fluids as guide to the discussion in Section 5.1.1. Empty fields repeat the entry 
above. PSD refers to the power spectral density in the optical trap, and Daq to the diffusion constant in aqueous solution. 



but significant variation of a from 0.81 in the interphase to 
0.85 during cell division. The larger exponents compared to 
the earlier work were attributed to possibly disrupted intracel- 
lular actin filaments there, and the difference in exponents dur- 
ing the cell cycle was interpreted as a sensitivity to changes of 
the cytoskeleton. 



From the corresponding video microscopy data in the in- 
terphase of S. pombe, time-averaged MSDs were computed 
from a small set of 8 long trajectories (up to 200 s) [232]. 
The MSDs cover time scales from 10 ms to 10 s and indicate 
subdiffusion with a relatively low exponent of or » 0.4. The 
ensemble- and time-averaged MSDs support this low exponent 
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for 0.1 s < f < 1 s, but the statistical noise permits a possible 
crossover to the previously found larger exponent at shorter 
times. Additional statistical quantities like higher cumulants 
or the mean maximal excursion are used to test the shape of the 
probability distribution of the displacements. It was concluded 
that the investigated data set is neither compatible with the 
CTRW model nor with a percolation scenario, but shares some 
features with fractional Brownian motion. Further improve- 
ments of the laser-tracking setup yielded very long trajectories 
with 22 kHz resolution of lipid granules during cell division 
of S. pombe [233]. The time-averaged MSDs at time scales 
between 10""^ to 1 s, clearly exhibit a crossover around 3 ms 
from normal diffusion to subdiffusion with exponent fi x 0.2, 
while the short-time MSDs obtained earlier [230, 231] display 
subdiffusion with the exponent ax 1-/3. Such a peculiar be- 
haviour is recreated in terms of a CTRW model with truncated 
power-law waiting-time distribution and has been related to 
the confinement by the optical trap. We refer the reader to the 
original literature for a detailed discussion [231-233]. Track- 
ing the granules by video microscopy at time scales above 
10 ms yielded subdiffusive motion with a x 0.8, similarly as 
in the earlier experiments. 

Summarising, the transport of granules in fission yeast has 
been studied with great precision over almost a decade and 
displays a complex behaviour varying over the range of ob- 
servation time scales. The experimental data are partially de- 
scribed in terms of FBM and a properly adapted CTRW model, 
but a consistent phenomenological description valid at all time 
scales is yet to be found. The microscopic mechanism leading 
to subdiffusive motion has not yet been resolved. 

In a highly recognised experiment addressing the physi- 
cal nature of bacterial cytoplasm, Golding and Cox [234] used 
video microscopy at 1 Hz to follow the motion of fluores- 
cently labelled mRNA macromolecules (radius on the order 
of 100 nm) over up to 30 minutes in E. coli, a rod-shaped bac- 
terium about 2 |jm long and half a micron in diameter. The 
analysis of 70 trajectories yielded subdiffusion at time scales 
between 1 to 30 s with varying prefactor, but rather robust 
exponent a - 0.70 + 0.07. Concatenating all trajectories, 
a power spectrum was obtained over three decades in time, 
which clearly indicates subdiffusive behaviour with exponent 
a = 0.77 ± 0.03. These findings agree with the limiting value 
of a in artificial crowding conditions using dextran [199], see 
Section 5.1.2. 

Subdiffusion with a significantly smaller value a - 0.39 + 
0.04 was obtained from time-lapse fluorescence microscopy 
tracking of GFP-labelled chromosomal loci in E. coli and 
Caulobacter crescentus over 10^ s [235]. For RNA proteins 
in E. coli, subdiffusion with a - 0.71 +0.10 was observed 
in agreement with Ref. 234. These findings were rationalised 
by supplementing the Rouse model for polymer dynamics by 
a fractional gaussian noise term [257] and arguing that the 
monomer dynamics at short times is then characterised by the 



exponent a' - a/2. Further, a CTRW scenario was ruled out 
as time- and ensemble-averaged data coincide. Analysis of the 
VACF showed pronounced anti-correlations and a description 
in terms of a Langevin equation with power-law correlated, 
spatially gaussian noise was favoured. 

Tracking the motion of telomeres in the nucleus of human 
osterosarcoma cells in a broad time range of almost 6 orders 
of magnitude displayed subdiffusive motion with varying ex- 
ponent depending on lag time [236]. Between 0.1 ms and 1 s, 
a small value of a ~ 0.32 + 0.12 was found, changing to 
0.5 1 + 0.20 in the intermediate regime up to 200 s, and ap- 
proaches normal diffusion for longer times; the amplitude of 
the power laws varied over an order of magnitude from one 
telomer to another An interpretion in the framework of the 
CTRW was found to not fully explain the experimental results, 
mainly since ageing behaviour could not be corroborated. The 
dynamics is, however, compatible with the dynamics of entan- 
gled polymers: the reptation model for a Rouse chain predicts 
an exponent of 1/4 for the short-time relaxation, crossing over 
to 1/2, and finally free diffusion. 

The interaction with cytoskeletal fibres and motor proteins 
was probed by engulfing 3 |jm polystyrene beads into vari- 
ous eukaryotic cells and monitoring their motion by video mi- 
croscopy [237]. The obtained MSD was found to grow super- 
diffusively with a ^ 3/2 at time scales less than 3 s, while 
a crossover to subdiffusive motion was observed for longer 
times (0.5 < a < 1). The latter finding, however, is subject 
to statistical noise and could be followed over less than half 
a decade only. The superdiffusive transport, dP'it) ~ P^^, 
was explained in terms of a generalised Langevin equation 
that accounts for the stochastic driving by motor proteins and 
the power-law memory due to the subdiffusive thermal fluc- 
tuations of semiflexible polymers, 6r^(t) ~ P^*. Further, the 
crossover to subdiffusion with a - 1/2 was rationalised by the 
decorrelation of motor activity at long time scales. 

Owing to technological advances in the 1990s, the FCS 
technique has quickly become an established tool for dynamic 
studies at the mesoscale, permitting also intracellular appli- 
cations. In vivo measurements benefit from the high tempo- 
ral resolution as well as from the choice of tracers ranging 
from tiny fluorescent proteins to labelled polymer coils and 
inert nanoparticles. Such tracers can be incorporated into the 
cell with only mild physiological side effects. In a pioneer- 
ing work, Schwille et al. [238] established the application 
of FCS to the cytoplasm by studying the diffusion of tetra- 
methyl-rhodamine dye with single molecule sensitivity in var- 
ious mammalian and plant cells. They found subdiffusion in 
the different cell types with a » 0.6, but the same data may 
equally well be rationalised by fitting a mixture of two nor- 
mally diffusing components, the faster one being 5-fold slower 
than in aqueous solution. The slow component showed diffu- 
sion constants up to 40 times smaller than the fast one and 
comprised 35-60% of the molecules; the slow component was 
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attributed to membrane-bound dye as the diffusion coefficients 
are of similar magnitude. 

The variation of anomalous transport with the spatial posi- 
tion in the cell was addressed using genetically modified COS- 
7 and AT-1 cells that express fluorescent EGFP proteins [239]. 
The obtained FCS correlations displayed anomalous diffusion 
depending on the position in the cell, the largest anomalies 
were observed in the nuclei. The data were analysed in terms 
of subdiffusion and two diffusive components. In the first case, 
the exponent a varied between 0.7 and 1 with half-value times 
around 1 ms. The two-component fits yielded a ratio between 
the diffusion constants of about 10, where the larger value was 
5-fold reduced compared to free diffusion of EGFP in aqueous 
solution. 

Weiss et al. [240] introduced differently sized FITC- 
labelled dextrans in HeLa cells and characterised their mo- 
tion with FCS. For all molecular weights studied, the corre- 
lation functions displayed anomalous transport and were fit- 
ted with the gaussian subdiffusion model, Eq. (93). The ob- 
tained values for the exponent a varied between 0.71 and 0.84 
non-monotonically depending on the size of the dextran poly- 
mers, which covered hydrodynamic radii in buffer solution be- 
tween 1 .8 and 14.4 nm. The dwell times To, on the other hand, 
increased systematically from 0.4 to 16 ms. Complementing 
their study by in vitro experiments with unlabelled dextran as 
crowding agent, the authors found a systematic decrease of a 
with the concentration of dextran, which suggests to quantify 
the degree of crowdedness in terms of the subdiffusion expo- 
nent a. 

Probing the cytoplasm and the nucleus of living HeLa cells 
as well as healthy and cancerous liver cells at the nanoscale 
by introducing fluorescently tagged 5 nm-sized gold beads re- 
vealed pronounced anomalous transport over five decades in 
time [241]. One advantage of using gold particles instead of 
(branched) polymers as probe is to reduce interactions with 
the cytoplasm, like polymer entanglement, potentially modi- 
fying the tracer transport. The FCS data were fitted with an 
empiric model that interpolates between subdiffusion at short 
times and normal diffusion at long times, see Eq. (91a) with 
6P-{t) - 6w'^[t/Td + (t/T,)"]. It consistently yielded a * 0.53, 
Tj X 0.3 ms, and tj ss 90 ms across all cell types, with 
a being slightly larger in the nucleus; these parameters im- 
ply a crossover time scale between the two transport regimes, 
tx - {TdT'^^y^^^ "^ ~ 56 s, which is beyond the scope of cur- 
rent experiments. When the cells were osmotically stressed by 
adding high, but non-aptotic concentrations of sucrose, raffi- 
nose, or NaCl, the subdiffusive transport changed markedly: 
a attained values around 0.66, while the dwell time t, in- 
creased by a factor of 3 in the cytoplasm and up to 2.4 ms 
in the nucleus. Considering the cyto- or nucleoplasm a con- 
centrated polymer solution, the observations were rationalised 
by a change of solvent conditions qualitatively modifying the 
visco-elastic properties of the medium. 



An extension of this study [242] covered a representative 
collection of mammalian cells, with different origin (from Chi- 
nese hamster ovary (CHO) cells to human osteosarcoma cells) 
and state of health (immortalised Thle-3 from a healthy liver 
vs. highly aggressive human glioma cells). Again, the gold 
nanobead displayed subdiffusive motion with a ^ 0.52 in the 
cytoplasm and a x 0.58 inside the nucleus, with only small 
variations across cell lines; the dwell times varied between 
0.1 and 0.3 ms with some exceptions reaching up to 0.9 ms 
in the nucleoplasm. In conclusion, very different mammalian 
cells show a similar degree of crowding at the nanoscale with 
the nucleus being somewhat less crowded than the cytoplasm. 
In particular, differences in the macroscopic visco-elastic re- 
sponse depending on a cell's development and disease state 
appear to be not reflected in the transport at the nanoscale. 

Dauty and Verkman [243] investigated the size-dependent 
transport of DNA in the cytoplasm of living HeLa cells with 
DNA molecules sized between 20 to 4,500 bp and labelled 
with a single fluorophore. The motion of DNA was followed 
by FCS, the fitting of correlation functions for sizes above 
250 bp required a two-component model with a short diffu- 
sion time of 5-20 ms independent of DNA size; the diffusion 
time of the slow component ranged up to ^ 500 ms for 3 kbp. 
Various sources were suggested as origin of the fast compo- 
nent, yet it could not be identified unambiguously. DNA dif- 
fusion was found to be significantly reduced in comparison 
to the free diffusion in saline by factors of up to 40, with a 
pronounced dependence on molecular weight above 500 bp. 
This dependence, however, almost disappeared by disrupting 
the actin cytoskeleton: the FCS data exhibited simple diffusion 
with an approximately 5-fold reduction. The authors have cor- 
roborated their findings by in vitro experiments using crowded 
solutions, cytosol extracts, and reconstituted actin networks. 
Only for the actin networks, the suppression of the DNA dif- 
fusion constant over its free value was sensitive to the molec- 
ular weight, recreating the behaviour in intact HeLa cells. In 
all other environments, simple diffusion was observed with re- 
duction factors not exceeding 5. The authors concluded that 
mobile obstacles can not explain the strongly reduced mobility 
of DNA in living cells and that the actin cytoskeleton presents 
a major restriction to cytoplasmic transport. Further, the sen- 
sitivity of diffusion to the molecular weight may be explained 
by entanglement effects with the actin mesh and reptation dy- 
namics of the elongated DNA molecules. 

Transport of globular actin molecules, a building block 
of the cytoskeleton, was investigated by FCS in living Dic- 
tyostelium discoideum cells [244] using mutants that express 
either GFP-labelled actin or free GFP. Measurements in the 
cytoplasm displayed anomalous transport for both probes, fits 
of the gaussian subdiffusion model yielded in both cases a ~ 
0.83 with half-value times of 170 [is for free GFP and 240 |js 
for G-actin. In the lysate of the actin mutant, transport was 
even slower with a 0.72, but with a half-value time reduced 
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by a factor of 2.5; here, a two-component fit was not able to 
detect an anticipated slow diffusive component. Further, a low 
degree of actin polymerisation was found in the cytoplasm, 
and thus, the slow transport was mainly attributed to crowding 
effects. Transport in the presence of an highly polymerised 
actin network was monitored by means of fluorescent LIM 
proteins marking the actin polymerisation front in the cell cor- 
tex. While transport was normal in the corresponding lysate, 
subdiffusive motion with a = 0.79 + 0.02 was observed in the 
living cell. Given that an exponent 3/4 is expected for in vitro 
actin networks [193], the observed anomalous transport may 
be caused by polymer fluctuations in this case. 

Luby-Phelps et al. applied the FRAP technique to fluo- 
rescently labelled dextran [245] and Ficoll [246] molecules 
microinjected into the cytoplasm of living Swiss 3T3 fibrob- 
lasts. The tracer radii ranged from 3 nm to 25 nm (Ficoll) and 
58 nm (dextran). The obtained diffusion constants were sub- 
stantially reduced relative to their values in aqueous solution, 
depending on the tracer size; a 30-fold reduction was observed 
for the largest Ficolls. The initial fluorescence was fully re- 
covered for molecules smaller than 14 nm (dextran) or about 
20 nm (Ficoll), indicating that all tracers are mobile. For larger 
molecules, the fluorescence recovery was increasingly incom- 
plete, and it was inferred that a fraction of the molecules is 
trapped. These observations are consistent with the picture of 
a significant excluded volume hindering cytoplasmic transport. 

Targeting the same experimental system, Seksek et al. 
[247] made different observations using an improved FRAP 
apparatus with microsecond resolution and measurement 
times of several seconds. The diffusion constants of labelled 
dextran and Ficoll molecules in the cytoplasm of T3T fibrob- 
lasts as well as MDCK epithelial cells were only 4-fold re- 
duced relative to aqueous solution. In particular, no size de- 
pendence of this factor could be detected for a comparable 
range of tracer radii, in contrast to the results in Refs. 245 
and 246. As in these works, the FRAP analysis was based 
on the half-value times of fluorescence recovery rather than 
on fits of the full recovery curve. Nevertheless, investigation 
of the recovery dynamics at long times did not display devi- 
ations form normal diffusion. In agreement with the earlier 
experiments, fluorescence was only partially recovered for in- 
creasing molecular weight, even at the scale of seconds. This 
is again consistent with the presence of spatial heterogeneity in 
the form of persistent structures or obstacles that immobilise a 
significant fraction of the larger macromolecules. 

A derivative of the FRAP technique is modulated fringe 
pattern photobleaching, which was employed by Arrio-Dupont 
et al. [248] to investigate the mobility of different globular pro- 
teins in cultured myotubes. The studied myotubes are elon- 
gated cells, 10-40 |jm wide and up to 1 mm in length. Apply- 
ing fringe patterns of 10-17 pm wide stripes yielded recovery 
time constants of several seconds for labelled proteins with 
hydrodynamic radii between 1.3 and 5.4 nm. These time con- 



stants translate into a 3 to 7-fold suppression of protein mo- 
bility relative to aqueous solution. For yS-galactosidase, which 
has a hydrodynamic radius of ai7.2 nm, a drastic suppression 
by 4 orders of magnitude compared to aqueous solution was 
found from the extremely slow fluorescence recovery with a 
time constant of 208 s. The results clearly display a monotone 
dependence of the suppression factor on protein size. Similar 
results were obtained for a series of differently sized dextran 
molecules in the same cell type [249], but the mobility is more 
size-sensitive for globular proteins. The apparent discrepancy 
with the findings of Ref. 247 may be explained by the very dif- 
ferent time and length scales probed by both experiments, and 
in this sense they complement each other. 

Combining FRAP and a photoactivation technique allowed 
for a non-invasive measurement of the motion of endoge- 
nously expressed GFP in the cytoplasm of E. coli bacte- 
ria [250]. The apparent diffusion constant was obtained as 
7.7 + 2.5 |jm^/s at the time scale of some 0.1 s, which is about 
1 1 times smaller than for free diffusion in water and signifi- 
cantly lower than in eukaryotic cells. 

5.7.2. Crowded model solutions fin vitro) 

Important progress in the understanding of anomalous trans- 
port was made by biologically motivated model systems, 
where in contrast to living cells key parameters are adjustable. 
In many cases, the crowded solutions are made from linear 
or branched (bio-)polymers of different rigidity, but proteins 
were used as crowding agents just as well. 

A component of the cytoskeleton is filamentous actin, an 
elongated, semiflexible polymer forming dense, entangled net- 
works, e.g., in the cell cortex. Solutions of reconstituted 
F-actin offer a model system to investigate the mechanical 
properties of these networks in vitro. The thermal motion 
of colloidal tracer particles (radius 0.25 pm) in semidilute F- 
actin solutions was tracked with video microscopy at 30 Hz 
over 20min [258]. The obtained MSDs display subdiffusion 
on time scales between 0.1 and 10 s with an exponent vary- 
ing between and 1 upon varying the actin concentration 
and thereby the mesh size of the networks between 0.17 and 
0.75 |jm. Normal diffusion was observed for mesh sizes larger 
than the tracer radius. The results can not be interpreted in 
terms of the macroscopic frequency-dependent shear viscos- 
ity. Analysis of the trajectories suggests that the tracers are 
caged by the network and the motion resembles a series of in- 
frequent and large jumps with power-law distributed waiting 
times. 

Diffusing wave spectroscopy experiments [259] on 
0.48 |jm-sized colloids embedded in F-actin solutions explored 
the large temporal window between 1 ps and 100 s. The MSD 
shows normal diffusion at short times up to 10 ms, where 
caging effects appear. The time-dependence of the diffusion 
coefficient is clearly observed and displays a suppression by 6 
orders of magnitude. 
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F-actin 


colloid (0.25 |jm) 


video microscopy 


0.1-lOs 


< Q- < 1 depending on 

ttipqIi *;i7p of aptiti tiPtworV 


2004 


ref. 258 




colloid (0.48 urn) 


diffusing wave 
spectroscopy 


l|js-100s 


crossover normal — > caging 
near 10 ms, D{t) spans 10* 


1999 


ref. 259 




nanosphere 
(25-250 nm) 


video microscopy 


50 ms- 10 s 


normal MSD, but 
anomalous, exponential 
propagator 


2009 


ref. 172 


fd virus 


protein (3.5 nm), 
silica spheres 
(35-500 nm) 


FCS, dynamic 
light scattering, 
video microscopy 


sensitivity: 10|js-1s 


D.JD: 1^10as[/a]T 


2005 


ref. 260 


dextran 


FITC-dextran 
( 10-500 kDa) 




T[/2. 0.3 —> jOms 


1 — > 0.6 monotonically 
as [dextran] t 


2UU4 


let. 240 




apoferritin 


FCS 


Ti/2 ~ 100 ms 
(675 [IS in saline) 


0.7-0.9 (mean 0.82), no 
ageing, stationary process 


2009 


ref. 261 




streptavidin, 
EGFP 


FCS 


1 ms-1 s 


1 — > 0.74 monotonically 
as [dextran] t 


2005 


ref. 199 


dextran 
(30 wt%) 


beads (50 nm) 


laser tracking 


50-500 ms 


subdiffusion {a = 0.82), no 
signs of ergodicity breaking, 
FBM most likely 


2012 


ref. 67 


Ficoll-70 


dye, proteins, 
polymers, beads 
(^1—1 uu nm ) 


FCS 


T1/2: 0.3 100 ms 
(rhodamine green) 


normal, D.^^jD > 100, 
exponential suppression: 

U Q.Xp\—(X [rlCOUJj 


2004, 
2008 


refs. 158, 
262 


micellar solution 
(C12E6) 


dye, proteins, 
nanoparticles 
(1.7-190 nm) 


FCS 


- 


normal, 2 components, 
i^aq/i^i: 1 — » 10 (tree 
tracers), D1/D2 « 10 


2006 


refs. 263, 
264 


dextran 


ribonuclease A 


FCS 


— 


normal, free and bound 

unjiciii. Ly{ j LJ2, ' 


2011 


ref. 265 


BSA protein 


BSA protein 
(3.6 nm) 


neutron scattering 


0.3-5 ns, 

wavelength: 0.5-3.3 nm 


normal, D.,q/Df,hoTt' 1 — > 10, 
hydrodynamic effects 


2010, 
2011 


refs. 266, 
267 


glycerol, PVP, 
Ficoll, proteins, 
E. coli lysate 


CI2 protein 
(7.4 kDa) 


nuclear magnetic 
resonance 




test of Stokes-Einstein 
relation for translational and 
rotational motion 


2010 


ref. 268 



Table III. Overview of in vitro experiments on crowded model fluids as guide to the discussion in Section 5.1.2. Empty fields repeat the 
entry above, right arrows indicate a systematic change upon variation of some parameter, up/down arrows indicate increase/decrease of some 
quantity. Abbreviations: FBM. . . fractional Brownian motion, [X]. . . concentration of X. 



In a related experiment, Wang et al. [172] tracked the mo- 
tion of nanospheres in semidilute F-actin solution. Choos- 
ing the tracer radius (25-250 nm) much smaller than the mesh 
size, normal diffusion was expected. The measured MSDs in- 
deed grow linearly in the window between 50 ms and 10 s. Yet, 
transport is anomalous: the van Hove function P(r, t) decays 
exponentially with r rather than following a gaussian. Like- 
wise, the propagator exhibits scaling upon rescaling distances 
as in normal diffusion, f oc rr^^^, but the resulting scaling 
function is anomalous. Similar findings were made for the 
quasi-one-dimensional diffusion of beads along phospholipid 
bilayer tubes. 

Another model system for a semiflexible biopolymer is 
the fd virus, being a stiff filamentous phage of contour length 
880 nm, diameter 6.6 nm, and persistence length 2.2 pm. Kang 
et al. [260] studied transport in semidilute suspensions of fd 
virus using differently sized tracers: a protein of radius 3.5 nm 



and labelled silica spheres with radii from 35 to 500 nm. Mea- 
surement of the long-time self-diffusion constant required suit- 
able techniques to cover the large range of particle sizes: 
fluorescence correlation spectroscopy for the protein and the 
small colloids, dynamic light scattering for intermediate-sized 
spheres, and video microscopy was found most appropriate for 
the large spheres. The obtained diffusion constant are system- 
atically suppressed by up to an order of magnitude when the 
concentration of fd increases; the effect is more pronounced 
for the larger tracers. Above the overlap concentration, i.e., 
in the semidilute regime, the fd rods span a network and an 
average mesh size can be assigned. The diffusion constants 
as function of the tracer radius divided by the mesh size, how- 
ever, display strong deviations from scaling — in contrast to the 
findings for cross-linked networks. This may be explained by 
the dynamically changing structure of the fd network, whose 
constituents themselves are subject to Brownian motion. 
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Aqueous solution of dextran molecules may also serve as 
a model system for a crowded fluid [240]. FCS measurements 
on a small fraction of labelled dextran molecules (10, 40, and 
500 kDa in weight) revealed a systematic decrease of the ex- 
ponent a with a concomitant increase of the half-value time as 
the crowding agent is more concentrated. The half-value times 
varied over an order of magnitude as a function of unlabelled 
dextran concentration, and the "anomaly parameter" a reached 
values below 0.6 for the 40kDa-dextran tracers. 

A benefit of using well-defined model systems is the pos- 
sibility to perform a detailed statistical analysis of anomalous 
transport at a given point in parameter space. FCS measure- 
ments on labelled apoferritin in a concentrated dextran solu- 
tion were repeated many times at the same spot to generate 
histograms of the residence times and of the exponent a [261]. 
The exponent covered values between 0.7 and 0.9 with mean 
0.82; a systematic variation of a, indicating ageing of the sam- 
ple, was not observed. These distributions were then critically 
compared to computer simulations of the FCS experiment for 
various models of anomalous transport: CTRW, FBM, and ob- 
structed diffusion on a percolating cluster (see Section 3), all 
of them generating the same subdiffusive MSD with (appar- 
ent) exponent a - 0.82. While FBM and obstructed diffusion 
provide a good description of the experimental data, the re- 
sults from the CTRW model are in qualitative and quantitative 
disagreement. This indicates that crowding-induced subdiffu- 
sion is more appropriately described by a stochastic motion 
with stationary increments than by a CTRW with untruncated 
waiting-time distribution. 

In the same lab, single-particle tracking experiments of flu- 
orescent 50 run beads in a solution of 30% 500 kDa-dextran 
were performed with temporal and spatial resolution of 4 ms 
and 10 nm, respectively [67]. For lag times between 50 ms 
and 500 ms, the time-averaged MSDs exhibited subdiffusion 
with an average exponent a » 0.82 as before. At about 1 s, a 
crossover towards normal diffusion was observed. The anal- 
ysis of 21 trajectories showed no signs of ergodicity break- 
ing, and comparing their geometric properties to the models 
of Section 3, fractional Brownian motion was found to match 
quantitatively. 

The monotone dependence of the subdiffusion exponent 
a with the concentration of crowding agent was corroborated 
by Banks and Fradin [199], who used the globular proteins 
streptavidin and EGFP as tracers in dextran solution. The con- 
ventional ^^-fit of the FCS data with the gaussian subdiffusion 
model was substantiated further: On a double-logarithmic rep- 
resentation, the long time decay of the FCS correlation is com- 
patible with the subdiffusion model down to the noise level of 
about 5 X 10"^ and clearly deviates from the power law r^^^, 
expected for normal diffusion. Second, an analysis assuming 
a large number of diffusive species, Eq. (98), yields a con- 
tinuous distribution of diffusion times between 0.1 and 10 ms, 
rather than, e.g., a bimodal distribution. The resulting expo- 



nents of subdiffusion decrease rapidly for increasing concen- 
trations of dextran obstacles and approach a = 0.74 + 0.02 
for high concentrations. This value appears to be related to 
the exponent 3/4 characterising the internal polymer dynam- 
ics of the dextran chains. Dextran aggregates are ruled out as 
source of anomalous diffusion, and the subdiffusive motion is 
remarkably robust against temperature variations, merely the 
half-value times reflect the changing viscosity of water Fur- 
ther, the motion of the small fluorescein molecules as well as 
of the dextran crowders itself shows simple diffusion even at 
high dextran concentration. If a globular protein is used as 
crowding agent, streptavidin motion is only slightly anoma- 
lous with a X 0.91 at high concentrations. 

Ficoll-70 served as crowding agent as well, being an in- 
ert, highly branched polysaccharide of approximately spheri- 
cal shape with a hydrodynamic radius of 5.5 nm. Verkman et 
al. [158, 262] characterised the transport of differently sized 
tracers in Ficoll-crowded solution using FCS with an illumi- 
nation region of 0.5 |jm in diameter. The tracers covered about 
two decades in size and ranged from rhodamine green over la- 
belled proteins, dextrans, and DNA fragments to fluorescent 
polystyrene beads with up to lOOnm in diameter, the latter 
mimicking cellular organelles. Although the FCS correlations 
displayed simple diffusion in all cases, tracer transport slows 
down drastically upon systematically crowding the solution. 
The obtained diffusion coefficients are suppressed by two to 
three orders of magnitude as Ficoll concentration increases up 
to 60 wt% and follow an exponential law. All tracers showed 
qualitatively the same behaviour, independently of their size. 
Using the smaller glycerol as crowding agent, the reduction of 
diffusion is smaller, but still exponential. For the large trac- 
ers, the slow diffusion can to a large extent be explained by 
the change of the macroscopic viscosity of the crowded fluid. 
This correlation is less pronounced for the small rhodamine 
green molecule, which appears to sense the microviscosity of 
its local environment. 

Another model for a crowded environment is given by 
aqueous micellar solution made of the non-ionic surfactant 
CiaEft. Szymahski et al. [263] studied diffusion constants 
of TAMRA-labelled lysozyme proteins with FCS for a wide 
range of C12E6 concentrations and found systematic deviations 
from free diffusion for concentrations above 3 wt%. A fit of 
the FCS correlation data, however, required a two-component 
diffusion model with only 10% of the proteins in the slow com- 
ponent; a possible explanation for the slow diffusion is the for- 
mation of protein-micelle complexes. Slow and fast diffusion 
constants are both suppressed by up to an order of magnitude 
at the highest concentration studied (37 wt%) and are separated 
by another order of magnitude. The gaussian subdiffusion 
model fits the data equally well yielding a systematic decrease 
of the exponent a down to 0.88, but the interpretation in terms 
of two components was favoured due to its clear physical inter- 
pretation. Focusing on the unbound fraction, the exponential 
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suppression of diffusion with concentration was consistently 
reproduced for a variety of fluorescent tracers covering sizes 
from 1.7 to 190 nm [264]. For the larger tracers, the diffusion 
constants are determined by the macroscopic viscosity of the 
micellar solution via the Stokes-Einstein relation. In contrast, 
the diffusion constant of small tracers show a high sensitivity 
to the tracer size: in the concentrated solution, apoferritin dif- 
fuses about 100-fold slower than lysozyme although the radii 
of the two proteins differ only by a factor of 3.6. The crossover 
length scale of 17 nm was identified with the persistence length 
of the micelles. 

The interplay of binding and crowding was addressed 
by means of a model system of negatively charged dextran 
molecules (500 kDa) as crowders and the positively charged 
protein ribonuclease A as tracer [265]. The tracer protein 
binds reversibly and non-specifically to a dextran molecule. 
Already at dilute concentrations of dextran (1 |jM), the FCS 
correlation signals the presence of a large fraction of bound 
protein with about 7-fold suppressed diffusion constant, re- 
flecting the increased hydrodynamic radius of the compound. 
Protein diffusion in solution of positively charged and neutral 
dextrans is unaffected at these low concentrations and reduced 
by merely a factor of 2 at the highest polymer concentration 
studied (lOOpM), crowding effects are still small. The data 
are fully explained by the coexistence of free and bound pro- 
tein, both of them exhibiting simple diffusion. 

Neutron scattering is a non-invasive technique to access 
protein solution samples at high protein concentrations at 
nanosecond and nanometre scales. Using quasi-elastic neu- 
tron backscattering, Roosen-Runge et al. [266, 267] probed 
the self-diffusion of bovine serum albumin (BSA) proteins in 
crowded aqueous solutions, where the same protein served as 
crowding agent. The protein motion was inferred from the 
incoherent dynamic structure factor P{k, co) of the hydrogen 
atoms, and the data are compatible with simple diffusion for 
wavelengths between 0.5 and 3.3 nm and time scales between 
0.3 and 5 ns. Increasing the protein content up to 30% volume 
fraction, a 10-fold reduction of the translational short-time dif- 
fusion coefficient over its value in dilute solutions was ob- 
served. Previous theoretical and experimental studies for hard- 
sphere colloids showed that diffusion is considerably slowed 
down already at short times when hydrodynamic interactions 
are taken into account, followed by a further slowing down 
at long times for high volume fractions [269, 270], see also 
Ref 271 discussed in Section 5.1.3. The reported BSA diffu- 
sion constants compare well with the findings for short-time 
diffusion of hard spheres if an effective hydrodynamic protein 
radius is used, which accounts for the hydration shell and the 
oblate ellipsoidal protein shape; the effective radius was de- 
termined to 3.6 nm from small-angle X-ray scattering. Noting 
further that at time scales of a few nanoseconds the caging 
by surrounding proteins is expected to be not yet effective, 
this suggests that the observed slow-down of protein diffusion 



is mainly caused by hydrodynamic interactions and that pro- 
tein diffusion in a crowded environment cannot be understood 
merely by excluded volume and confined motion. 

Nuclear magnetic resonance spectroscopy was employed 
to quantify both the rotational and translational diffusion of 
the protein chymotrypsin inhibitor 2 (CI2) in a variety of 
crowded solutions as function of crowder concentration [268]. 
The crowding agents comprised glycerol, synthetic polymers 
(PVP, Ficoll), globular proteins (BSA, ovalbumin, lysozyme), 
and E. coli cell lysates. The macroscopic shear viscosity of 
the solution, 77, increases with crowder concentration. Then 
the Stokes-Einstein relation already suggests a reduced dif- 
fusivity, D oc 1//;, which describes the measured CI2 diffu- 
sion constants in glycerol solution and also translational dif- 
fusion in ovalbumin, BSA, and cell lysate. In solutions of 
synthetic polymers, transport is affected less than expected 
from the increase of viscosity, and the translational motion 
of CI2 is impeded more than its rotational motion. Surpris- 
ingly, the opposite effect was found in protein-crowded solu- 
tions and in the cell lysate: rotational diffusion constants were 
suppressed stronger than translational diffusion and stronger 
than the Stokes-Einstein relation would imply. The findings 
were attributed to weak non-specific, non-covalent chemical 
interactions between proteins, while synthetic polymers tend 
to form a loose mesh work. This suggests that crowded pro- 
tein solutions are preferred for modelling the effects of the in- 
tracellular environment on protein transport. 

5.1.3. Computer simulations (in silicoj 

Anomalous transport and subdiffusion were unambiguously 
identified in computer simulations for minimalist, generic 
models as well as for elaborate cytoplasm models targeting 
a specific cell type. These results provide essential support for 
the interpretation of experiments and test various microscopic 
pictures behind anomalous transport. 

A paradigm of anomalous transport is provided by the mo- 
tion of a small tracer in a disordered matrix of spherical ob- 
stacles, known as Lorentz models [79, 115, 126, 127]; see 
Section 3.3 for a thorough discussion and its connection to 
continuum percolation theory. As a localisation transition is 
approached by tuning the excluded volume, a growing win- 
dow of subdiffusion emerges at intermediate time scales. The 
tracer dynamics generically displays a double-crossover sce- 
nario from microscopic motion at very short time scales to 
subdiffusion and to normal diffusion at large times [115]. This 
phenomenology is preserved, whether the microscopic tracer 
motion is ballistic like in a porous medium or overdamped 
like in the cytoplasm [79]. The apparent exponent of the sub- 
diffusive regime decreases monotonically upon approaching 
the transition and converges to its universal value, a x 0.32 
or 0.42, depending on whether the tracer is restricted to the 
percolating space or not [127]. Subdiffusive transport over 6 
decades in times at the critical obstacle density and a suppres- 
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crowded environment 


tracer 


interactions 


subdiffusion exponent a 
and temporal window 


year 


reference 


10^ spherical obstacles 
(overlapping, equal size) 


point (ballistic and 
overdamped motion) 


hard core 


1 ^ 0.32 ^ 0, 
subdiffusion over up to 
6 decades in time 

1 — > U.4Z — > U 

(percolating cluster only) 


2006, 
2008 

ZUi 1 


refs. 79, 115 
rei. iz/ 


1000 spheres (equal size), 
fluid in obstacle matrix 


fluid particles 


hard core 


1/2 over 4 decades in time 


zooy, 

2011 


rets. lUU, lUl 


1000 spheres (bidisperse, 
partially quenched) 


mobile particles 


hard core 


0.3 over 4 decades in time 


2009, 
2010 


rets. 102, 103 


2000 spheres (size-disparate 
mixture) 


small particles 


soft repulsion 


0.63 and 0.5 over several 
decades 


2006, 
zuuy 


refs. 272, 273 


5000 spheres with weight 
distributed modelling HeLa 


small spheres 
(2-5 nm) 


hard core, soft repulsion 


1 — > 0.55, subdiffusion 
over 3 decades in time 


2004 


ref 240 


atomistic E. coli model: 
1000 macromolecules of 
50 different types 


macromolecules 
(10-lOOOkDa) 


repulsion, vdW attraction, 
electrostatic, hydrodynamic 


0.75-0.85 @ lOOps-lns 
(full set of interactions) 


2010 


ref 274 


E. co/i model: 15 molecule 
types, 1300 molecular-shaped 
or spherical particles 


macromolecules 


repulsion, vdW attraction, 
hydrodynamic (spheres 
only) 


£»o/Aho„: 3-10 (with 
HI), anomalous— >normal 
at 1 |is, Dshort/f ~ 2 (with 
HI) or 20 


2010 


ref. 271 


1000 spherical obstacles 
(non-overlapping, equal size) 
and 10* solvent particles 


globular polymer 
(400 beads) 


repulsion, vdW attraction, 
hydrodynamic (MPCD) 


1-0.7 depending on 
excluded volume and 
polymer size 


2010 


ref 275 


polymeric network 


free monomer 


repulsion, vdW attraction 


1/2 or 3/5 over up to 2 
decades in time 


2011 


ref 276 



Table IV. Overview of computer simulations on crowded fluids as guide to the discussion in Section 5.1.3. Empty fields repeat the entry 
above, arrows indicate a systematic change upon variation of some parameter. Abbreviations: vdW. . .van der Waals, HI. . .hydrodynamic 
interactions, MPCD. . . multi-particle collision dynamics. 



non-vanishing non-ergodicity parameters due to trapped parti- 
cles [100-103, 126, 273, 277] that would appear as immobile 
fraction in a FRAP experiment, see Section 4.3. 

Aiming at a more realistic modelling, the measured distri- 
bution of protein sizes in the cytoplasm of HeLa cells served 
as input for a simple cytoplasm model consisting of 5000 
spherical particles with exponentially distributed molecular 
weights [240]. The particles interact via a hard core plus a soft 
repulsive, parabolic potential; their overdamped motion was 
followed by Brownian dynamics simulations. The obtained 
MSDs display subdiffusion over 3 decades in time, where the 
exponent a decreases as function of the particle radius in a 
sigmoidal fashion from unity for the smallest particles (2 nm) 
to about 0.55 for the largest radii studied (5 nm). Thus, the 
excluded volume can, to a large extent, account for the exper- 
imentally observed subdiffusion in HeLa cells. 

A step towards a realistic modelling of the E. coli bac- 
terium was taken by McGuffee and Elcock [274], who de- 
veloped an atomistically detailed model of its cytoplasm in- 
cluding 50 different types of macromolecules at physiolog- 
ical concentrations. For systems of 1000 macromolecules 
in total, Brownian dynamics simulations including electro- 
static and hydrodynamic interactions yielded system trajec- 



sion of the time-dependent diffusion constant by 5 orders of 
magnitude was observed in the simulations [126, 127]. Even 
10% off the critical obstacle density, the motion is subdiffusive 
with exponent a » 0.5 over still 2 decades in time. Thus for the 
observation of anomalous transport over a finite time window, 
it is not essential that the system is fine-tuned to the critical 
point. The scenario of the localisation transition appears to be 
robust: subdiffusive motion with exponent 1/2 over 4 decades 
in time was found in simulations for a hard-sphere fluid ad- 
sorbed in a matrix of disordered, non-overlapping hard-sphere 
obstacles [100, 101]. Simulations for a similar system where 
the tracer fluid and the obstacles are correlated yielded sub- 
diffusion over 4 decades in time with an exponent 0.3 resem- 
bling the value from the Lorentz model [102, 103]. Anoma- 
lous transport was even found in simulations for transport of 
small particles in a glassy, slowly rearranging matrix of large 
particles [272, 273]; although the MSD does not display un- 
ambiguous power-law behaviour, the apparent exponent a(t), 
defined by ait) - dlog(5r^(f))/dlog(f), was found to be 0.5 
and less over several decades in time [273], and there is evi- 
dence that the localisation scenarios in the binary mixture and 
the Lorentz model are related [277]. At high volume frac- 
tions of the obstacle matrix, all of the mentioned models share 
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tories over 6 million steps or a time span of 15 [is. Thereby, 
molecular behaviour inside a cell is illustrated by vivid movies 
of educational value. Further, anomalous transport was ob- 
served over 4 decades in time with local MSD exponents a{t) 
between 0.75 and 0.85. Focusing on the role of excluded vol- 
ume by including only repulsive interactions in the model, the 
anomaly was markedly reduced. Hence, the effects of macro- 
molecular crowding extend beyond those of excluded volume, 
although a ramified excluded volume appears to be essential 
for pronounced subdiffusive motion. 

The role of hydrodynamic interactions (HI) on protein mo- 
tion was elucidated by Ando and Skolnick [271], who per- 
formed Brownian dynamics simulations of models for the 
E. coli cytoplasm comprised of 15 different macromolecule 
types: proteins were represented either by their molecular 
shape or by equivalent spheres of the same hydrodynamic ra- 
dius, the particles interact non-specifically via a soft repulsive 
parabolic potential or via attractive van der Waals forces, and 
hydrodynamic interactions were included for the equivalent- 
spheres model. It was found that HI yield a 3- to 10-fold 
reduction of the short-time diffusion constants, depending on 
protein radius and concentration. This finding is supported 
by previous results for monodisperse spheres [269, 270] and 
corroborates recent experiments [266, 267], see above. Long- 
time diffusion constants were determined as function of the 
protein radius at an observation time scale of 5 |as and display 
a further reduction by a factor of 2 (with HI) or up to 20 (repul- 
sion only). Notably, the equivalent-sphere model including HI 
reproduced the experimentally observed diffusion constant of 
GFP in vivo without any adjustable parameters. As a key re- 
sult, excluded volume effects and hydrodynamic interactions 
are likely the two major factors for the large reduction in dif- 
fusion of macromolecules observed in vivo. 

Echeverria and Kapral [275] studied the effect of crowd- 
ing on the conformation and transport of a globular polymer 
in an explicit, albeit mesoscopic, poor solvent. The crowded 
environment was formed by a frozen, random array of hard 
spherical obstacles, and the system was propagated using the 
multi-particle collision dynamics (MPCD) scheme (which in- 
trinsically includes hydrodynamic effects) for up to 1 million 
solvent particles. The simulation results show that the equilib- 
rium structure of long polymer chains is significantly altered 
towards a chain of polymer blobs being trapped in voids be- 
tween the obstacles. For increasing volume fraction of ob- 
stacles, subdiffusive motion of the centre of mass emerges 
with exponents a down to 0.7, and it appears that a decreases 
further as the polymer size increases. These findings are in 
accord with earlier lattice studies for polymers in disordered 
media [278-280], which showed that entropic trapping of the 
polymer in the voids leads to a stark slowing down of transport 
and to subdiffusive motion. 

Transient binding to polymeric networks, e.g., in hydro- 
gels, can render the transport of tracer molecules anoma- 



lous as well as demonstrated by molecular dynamics simula- 
tions [276]. Strongly bound tracers sliding along a relaxed 
polymer strand exhibit pronounced subdiffusive motion with 
exponents a = 1/2 or 3/5, depending on the fractal structure 
of the coiled strands. For a swollen hydrogel, the polymer 
coils are stretched and the window of subdiffusion was found 
to decrease quickly. Further, there is a competition between 
sliding along the polymer and free diffusion within the pores 
of the network. The observed anomalous transport depends 
sensitively on the attraction strength between the tracer and 
the polymer chain, which governs this competition. 

5.2. Crowded membranes 

5.2.1. Cellular membranes (in vivo) 

According to the "fluid mosaic" model by Singer and Nicolson 
[299], the plasma membrane of cells is thought of as an essen- 
tially homogeneous fluid bilayer of phospholipids with freely 
diffusing protein inclusions. Numerous in vivo measurements 
of protein motion in membranes during the past two decades, 
however, required several revisions of this simple picture. To- 
day, the plasma membrane is considered a patchy lipid bilayer 
densely packed with integral and peripheral proteins, where 
some of the lipids are organised into microdomains, some of 
the proteins form clusters, and others are tethered to the cy- 
toskeleton [19, 300, 301]. 

First hints on the heterogeneous landscape of the plasma 
membrane go back to early FRAP experiments [207, 295], in- 
dicating that a significant fraction of certain membrane pro- 
teins is immobile. These findings were corroborated by single- 
particle tracking techniques emerging in the 1990's, which 
enabled the observation of individual fluorescently labelled 
membrane proteins [223, 281-283]. The dynamics of mem- 
brane receptors in mouse keratinocytes was found to be highly 
heterogeneous within the same type of protein, and the trajec- 
tories were classified after their MSDs into four types of mo- 
tion: free diffusion, confined diffusion, immobile, and directed 
motion [281]. 

SPT experiments on IgE receptors by Feder et al. [223] 
yielded another class of MSDs, which required a subdiffusive 
power-law fit for 56% out of 241 trajectories (150 frames ev- 
ery 1.6 s) with an average exponent a - 0.64 ± 0.45; 27% 
of the receptors were immobile. The findings are consistent 
with corresponding FRAP measurements, see below. In sim- 
ilar experiments. Sheets et al. [282] preferred to account for 
the subdiffusive fraction by refining the classification of tra- 
jectories also according to their "shape" and by introducing 
the additional class of slow normal diffusion. For the Thy-1 
protein, the classes of fast and slow diffusion yielded diffusion 
constants as different as 0.081 and 0.0035 pm^/s, respectively. 

Significant enhancement of the temporal resolution to 
0.22 ms of the SPT experiments [283] provided a more de- 
tailed picture of band-3 protein transport in erythrocyte mem- 
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cell type 


probe 


experimental 

tpphmniip 


temporal and 


observation 


year 


reference 


keratinocytes 


E-cadherin and 
transferrin receptor 


SPT, FRAP 


33 ms-30s 


heterogeneous MSDs: free, 
confined, immobile, directed 


1993 


ref. 281 


RBL cells 


IgE receptor 


SPT, FRAP 


1-100 s 


56% of tracers: subdiffusion 

■.-1 f\ /' C , f\ A C 

With a = 0.65 ± 0.45; 
27% immobile 


1996 


ref. 223 


C3H lOTl/2 
fibroblasts 


Thy-1 protein 


SPT 


30ms-l s 


normal diffusion, fast and 
slow components 


1997 


ref. 282 


erythrocytes 


band- 3 protein 


SPT 


0.22 ms-1 s 


£>iOm.s/As« 80-90, 
1/3 immobile 


1998 


ref. 283 


rat kidney fibroblasts 


DOPE phospholipid 


SPT 


25 |JS-3 s 


3 subdiffusion regimes 

(double crossover): 

a = 0.74 — > 0.55 — > 0.79 


2002 


ref. 284 


various cell types 








subdiftusion @ 0.1 -10 ms 
(a = 0.53), short- and 
long-time diffusion: 
Do/D « 10 


2004 


ref. 285 


HeLa cells 


MHC class I protein 


SPT 


4-300 s 


subdiffusion: a = 0.49 ± 0.16 


1999 


ref. 286 


CHO cells 


H (TT T 1 TT i 

MHC class 11 protein 


SPT 


lOOms-3 s 


normal diffusion 


2002 


ref. 287 


C. crescentus 


PleC protein 


SPT 


1-9 s 


normal: D » 10"^ [im-/s 


2004 


ref. 288 


COS-7 fibroblasts, 
MDCK cells 


AQP-1 protein 


SPT 


lOms-6 min 


heterogeneous MSDs (free, 
restricted, immobile), normal 
case: D depends on protein 
content 


2008 


ref. 289 


HEK 293 cells 


Kv2.1 potassium 
channel 


TIRFM 


O.l-lOOs 


subdiffusion (a » 0.8), 
ageing, fractal space, 
stationary and non-stationary 
displacements 


2011 


ref. 171 


RBL cells 


diI-Ci2 lipid 


FCS 


Ti/2 a 30 ms 


subdiffusion a = 0.74 ± 0.08, 
or two components: 
Ui/Ui = iU 


1999 


ref. 290 


oligodendroglia 


MOG protein, 
sphingomyelin lipid 


FCS 


0.01-100 ms 


subdiffusion of MOG protein 
(ty = 0.59), two-component 
diffusion of lipid 


2005 


ref. 291 


HeLa cells (intracel- 
lular membranes) 


Golgi-resident 
enzymes 


FCS 


0.3-100 ms 


subdiffusion: a » 0.5-0.8 


2003 


ref. 292 


yeast cells 


Fus-Mid-GFP 
protein 


scanning FCS 


Ti/2 a 2s 


slow, normal diffusion: 
D « IQ-^ nmVs 


2006 


ref. 293 


COS-7 cells 


transferrin receptor 
(TfR), sphingolipids, 
anchored proteins 


variable-area 
FCS 


beam waist: 
190-390 nm 


w-dependent effective D: 
Tip_{w) = to+ wP-jADau 
/o > for lipids and anchored 
proteins, /o < for TfR 


2005, 
2006 


refs. 181, 294 






variable- area 
FCS 


beam waist! 
75-350 nm 


3— 5-fold reduction of /)(vv) 
for small w 


zuu / 


rei. loj 


PtK2 cells 


sphingomyelin, 
anchored protein 


STED-FCS 


beam waist: 
15-125 nm 


anomalous vv-dependence: 
Ti/2(w — » 0) « 10ms 


2009 


ref. 189 


rat myotubes 


ACh receptor 


FRAP 


1 min 


immobile fraction of ACh in 
dense patches 


1976 


ref. 207 


human embryo 
fibroblasts 


WGA receptor 


FRAP 


i"i/2 a 30 s 


only 75% mobile receptors 


1976 


ref. 209 


thymocytes, 
lymphoma cells, 
fibroblasts 


Thy-1 protein 


FRAP 


ri/2 < 1 s 


a; 50% immobile 


1987 


ref. 295 


COS-7 cells 


anchored, acylated, 

clliU Ll ullAllll^lliL'i OLVlXj 

proteins 


FRAP 


bleaching 

^tnnp' ^ 1 irn 

SLllLft^. Lllll 


normal, D = 0. 1-1.2 pm^/s for 

HiFTPrpnt nrotPiTiQ 


2004 


ref. 296 




4 membrane proteins 


FRAP 


Ti/2 « 2s, spot 

radius: 3.6 ^im 


normal, but D 10-fold smaller 
than for GPI-proteins 


2007 


ref. 297 


HEK 293T 


CD4 and chemokine 
receptors 


variable-radius 
FRAP 


spot radius: 
1.40-3.45 Mm 


D and recovery ratio depend 
on size of bleaching spot 


2007 


ref. 298 



Table V. Overview of in vivo experiments on cellular membranes as guide to the discussion in Section 5.2.1. Empty fields repeat the entry 
above. Abbreviations: TIRFM. . . total internal reflection fluorescence microscopy, STED. . . stimulated emission depletion. 
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branes: a fraction of proteins was immobile, the remaining two 
thirds exhibited normal diffusion at macroscopic time scales 
of a few seconds. The diffusion constant, however, was dras- 
tically reduced by a factor of 80-90 compared to its micro- 
scopic value at 10 ms. Refined experiments suggested a hop- 
ping mechanism due to interaction with the spectrin network, 
characterised by a mesh size of llOnm and an average hop 
frequency of 2.8 s 

Some years later, the same lab implemented a high- 
speed SPT setup with a temporal resolution of remarkable 
25|js [284], which permitted the detailed tracking of DOPE 
phospholipids in rat kidney fibroblasts. The obtained MSD 
covers five(!) decades in time and exhibits two crossovers at 
12 and 590 ms, connecting three regimes of motion (see sup- 
plement of Ref 284). All three regimes are compatible with 
subdiffusive motion with exponents a - 0.74 (short times), 
a = 0.55 (intermediate times), and a - 0.79 (long times, 
but regime covers less than a decade). Extrapolating to short 
time scales, the MSD appears to approach the free diffusion 
in a homogeneous lipid bilayer It was concluded that the 
plasma membrane is doubly compartmentalised into 750 nm 
and then into 230 nm compartments with regard to the lateral 
diffusion of DOPE lipids. Repetition of the measurements for 
many other cell types [285] corroborated these findings qual- 
itatively with hopping rates between 1 and 59 s"' and com- 
partment sizes between 30 and 230 nm. For FRSK cells (fetal 
rat skin keratinocytes), the regimes at short and long times ex- 
hibited essentially normal diffusion with diffusion constants 
separated by an order of magnitude; the intermediate regime 
between 0. 1 ms and 10 ms was found to be subdiffusive with 
exponent a = 0.53. 

On the basis of these observations, Kusumi et al. [19, 302] 
have proposed a revised model for the plasma membrane, 
which should be viewed as a compartmentalised fluid. The 
structure is provided by barriers or obstacles like the actin- 
based membrane skeleton ("fences") and anchored transmem- 
brane proteins ("pickets"). The motion of proteins and lipids 
is then thought of as a hopping process between differently 
sized membrane compartments, and transport is slowed down 
by such specific barriers. 

Deeper insight into the stochastic transport process is ex- 
pected from the investigation of displacement histograms for 
fixed lag times, i.e., the van Hove function. Such a route 
was taken by Smith et al. [286], who performed SPT experi- 
ments using fluorescence imaging of major histocompatibility 
complex (MHC) class I molecules on HeLa cells. Fitting the 
propagator of normal diffusion to the displacement histograms 
yielded the diffusion coefficients D(t) as function of the obser- 
vation time scale. A marked decrease of D{t) with increasing 
time interval was found, and the data were best described by 
subdiffusive motion with a = 0.49+0. 16 over the whole acces- 
sible temporal window between 4 s and 300 s. Unfortunately, 
the histograms suffered from large statistical noise prohibiting 



an answer to the question whether the distribution of displace- 
ments is gaussian or not. 

The translational diffusion of MHC class II membrane pro- 
teins in Chinese hamster ovary cells was studied by Vrljic et al. 
[287] to probe the plasma membrane for barriers from putative 
lipid microdomains. Using SPT at a rate of 100 ms, the cumu- 
lative probability distribution function of displacements was 
extracted at time scales up to only a few seconds (limited by 
photobleaching). In this time window, however, almost negli- 
gible deviations from simple Brownian motion and no signifi- 
cant confinement could be inferred. Most notably, the authors 
also studied correlations in the transport of close pairs of pro- 
teins, initially separated between 0.3 and 1 |jm. The idea was 
to detect a possible confinement by impermeable, but diffus- 
ing barriers. Again, the results followed the predictions of free 
diffusion and no evidence for a restriction to small, freely dif- 
fusing domains was found at the length scales under investiga- 
tion. 

Normal diffusion was also observed by tracking the trans- 
membrane histidine kinase PleC in the bacterium Caulobac- 
ter crescentus [288], which yielded the relatively low value 
D = 12 X 10"^ [im^/s. The experimental MSDs, however, cov- 
ered only a small time window between 1 and 9 s (due to lim- 
itations of the fluorescent label) and are subject to large statis- 
tical errors. Hence, the question whether transport is normal 
or anomalous could actually not be addressed. Further, trans- 
port was found to be spatially uniform across all positions in 
the elongated cell with exception of the cell poles and without 
evidence for active transport. 

The impact of putative lipid microdomains on the motion 
of integral membrane protein aquaporin-1 (AQP-1) in COS -7 
fibroblasts and Madin-Darby canine kidney cells was studied 
by long-time SPT experiments [289]. Labelling with quan- 
tum dots enabled the recording of long trajectories at 1 Hz 
over 6min. Trajectories were classified into free, restricted, 
and immobile according to their deviation from the extrapo- 
lated short-time MSDs at 91 Hz. Most AQP-1 proteins dif- 
fused freely over several microns on the scale of minutes, and 
diffusion was faster by a factor of 4 in protein-poor membrane 
blebs. It was concluded that AQP-1 is a largely non-interacting 
protein with a macroscopic diffusion constant determined by 
the concentration of obstructions in the membrane. 

Only recently, Weigel et al. [171, 303] investigated the dy- 
namics of Kv2.1 potassium channels in the plasma membrane 
of human embryonic kidney (HEK 293) cells. Labelling with 
quantum dots and tracking by means of total internal reflection 
fluorescence microscopy (TIRFM) permitted the collection of 
1000 trajectories each about lOmin long at temporal and spa- 
tial resolution of 50 ms and 8 nm, respectively. Scrutinising 
time- and time-ensemble-averaged MSDs, ageing behaviour, 
the propagator, and the distribution of waiting times indicated 
the coexistence of a non-stationary transport process compat- 
ible with a CTRW and a stationary walk in a fractal space. 
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The non-stationary process was attributed to binding of the 
potassium channels to the actin network, while the stationary 
process may have its origin in the molecular crowding of the 
membrane. 

During the past decade, FCS has routinely been applied 
to cellular membranes for in vivo measurements of lipids and 
proteins. Introducing fluorescent lipid probes in the plasma 
membranes of rat basophilic leukaemia (RBL) cells, Schwille 
et al. [290] obtained FCS correlations which clearly deviate 
from normal diffusion. The data can be described by the 
gaussian subdiffusion model with a - 0.74 + 0.08 and half- 
value times around 30 ms, or alternatively well, by the two- 
component diffusion model yielding diffusion times of 3 ms 
and 100 ms. A control experiment with giant unilamellar vesi- 
cles (GUVs) showed normal diffusion, but inducing phase sep- 
aration of the lipid mixture by adding cholesterol rendered the 
lipid transport anomalous again. 

For myelin oligodendrocyte glycoproteins bound to oligo- 
dendrocyte membranes, anomalous diffusion with a - 0.59 
was reported [291], while a two-component diffusion model 
appeared to be less likely. The results for sphingomyelin lipids 
in the same membrane suggest a two-component description 
with diffusion constants of 0.37 and 70 |am^/s; here, a possible 
explanation for the highly mobile fraction are dye molecules 
moving freely in solution. 

Intracellular membranes were addressed by Weiss et al. 
[292], who monitored three Golgi-resident enzymes both in 
the endoplasmic reticulum and in the Golgi apparatus of HeLa 
cells. Subdiffusive motion was found over more than two 
decades in time with half-value times of 10-30 ms and with 
exponents a that depend on the type of enzyme and range be- 
tween 0.5 and 0.8; for the same enzyme, similar values of a 
were obtained in the two different membrane structures. The 
data can be equally well fitted by the two-component diffusion 
model, but it was argued that such a fit seems unlikely. 

A modification of the FCS technique is scanning 
FCS [304], where the illumination laser repeatedly scans the 
probe, e.g., along a large circle. The same spot is visited 
with a low frequency, giving access to the fluorophore dy- 
namics at very long time scales. Combining this method with 
continuous wave excitation, very slow, albeit normal diffu- 
sion was observed for Fus-Mid-GFP proteins on yeast cell 
membranes [293], with diffusion constants on the order of 
lO^^* [im^/s and FCS diffusion times of up to 2 seconds. Com- 
pared to the motion of these slow proteins, the other membrane 
components rearrange quickly, effectively homogenising the 
originally heterogeneous environment of the protein. This may 
explain the finding of simple diffusion at long times, resolv- 
ing an apparent contradiction with the measurements quoted 
above. The example further highlights the importance of tak- 
ing into account the experimental time and length scales for 
the interpretation of transport in complex environments. 

Complementary to the previous FCS studies, the spatial as- 



pects of complex molecular transport in membranes were ex- 
plicitly addressed by systematic variation of the detection area. 
Using beam waists ranging from 190 to 390 run, Wawrezinieck 
et al. [181] have measured the motion of the GFP-labelled 
transmembrane protein TfR (transferrin receptor) and of a sph- 
ingolipid in the plasma membrane of COS-7 cells. The ob- 
tained two sets of diffusion times exhibit an affine dependence 
on motivating the FCS diffusion law, t^(w) = fo + w^/4Dgf{. 
Extrapolation to zero beam waist yields the axis intercepts fo, 
which were found to be -1-25 ms for the lipids and -20 ms in 
case of the TfR proteins, with uncertainty of 10%. The signs of 
fo can give a hint on the transport mechanism, and it was con- 
cluded that the motion of sphingolipids is hindered by isolated 
microdomains and those of the TfR proteins by a meshwork 
of barriers. Alternatively, the protein diffusion times may be 
interpreted in terms of fractal scaling law t^(w) ~ w^'^ accord- 
ing to Eq. (95), indicating subdiffusion with or = 2/dv, ~ 0.77 
[23]. The findings of Ref. 181 were substantiated by fur- 
ther experiments with a variety of molecules [294], provid- 
ing a broad range of intercept values: sphingolipids and GPI- 
anchored proteins yielded clearly positive fo, the transmem- 
brane TfR-GFP is characterised by a strongly negative value, 
and the axis intercept vanished for glycerophospholipid ana- 
logues. In conclusion, the slow dynamics of sphingolipids is 
attributed to cholesterol and sphingomyelin levels in the mem- 
brane, and the transferrin receptor protein is likely to be dy- 
namically confined by actin-based cytoskeleton barriers. 

It is clear from Eq. (90) that the diffusion time must vanish 
as the observation area approaches zero. Thus, the extrapo- 
lation w — > of the above diffusion law can only hold ap- 
proximately, and it is of interest to elucidate the behaviour at 
small w, below the diffraction limit. Observation areas with 
radii ranging down to 75 nm were achieved by equipping the 
microscope coverslips with circular nanometric apertures in a 
metallic film [185]. This approach permitted the elucidation of 
nanometric membrane heterogeneity and a considerable exten- 
sion of the measured diffusion laws t^(w) of Ref. 294, reveal- 
ing non-trivial crossovers. For ganglioside lipid analogues and 
GPI-anchored GFP in the plasma membrane of COS-7 cells, 
the apparent, w-dependent diffusion constants w^/4t<^(w) at 
short scales are 5- and 3-fold reduced compared to their large- 
scale values, with a crossover length scale of w » 100 nm. 
Even more interestingly, the diffusion law of TfR-GFP shows 
two transitions around 150 and 230 nm; the latter transition 
is explained by a meshwork of cytoskeletal barriers hindering 
transport of the transmembrane protein. 

Application of STED-FCS to the plasma membrane of liv- 
ing mammahan cells gives access to nanosized detection areas 
down to w = 15 nm in radius [189]. Phosphoethanolamine 
molecules, which are assumed not to form molecular com- 
plexes with membrane components, exhibited free diffusion, 
Td{w) oc w^, over the full range of investigated areas. The dif- 
fusion of sphingomyelin, however, was found to be strongly 
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heterogeneous: only for large radii, w > 80 nm, unobstructed, 
normal transport was found. For small detection areas, w < 
40 nm, the FCS correlations showed clear deviations from 
normal diffusion. Fits with the subdiffusion model yielded 
a < 0.7, but an explanation of the transport in terms of two 
diffusive components was favoured. The fast fraction resem- 
bled the unobstructed transport of phosphoethanolamine, both 
in space and time. The diffusion times of the slow fraction 
approached a minimum of about 10 ms for small areas. This 
offset was explained by a brief trapping event of that dura- 
tion by cholesterol-mediated complexes within an area of less 
than 20 nm in diameter Similar findings were made for GPI- 
anchored proteins. 

The first applications of fluorescence recovery after photo- 
bleaching (FRAP) on plasma membranes addressed the lateral 
motion of fluorescently marked receptor proteins in rat my- 
otubes [207] and human embryo fibroblasts [209]. In both 
experiments, only about 75% of the initial fluorescence was 
recovered indicating that a significant fraction of receptors is 
immobile. Similarly, the fluorescence recovery was found to 
be only about 50% for labelled Thy-1 antigen in the plasma 
membrane of lymphoid cells and different fibroblasts [295]. 
The obtained diffusion constant of the mobile fraction of Thy- 
1 proteins was comparable to those of the lipid analogues in 
the bilayer, which is consistent with a putative lipid anchoring 
to the plasma membrane. Feder et al. [223] studied the motion 
of IgE receptors in the membrane of leukocytes using SPT and 
FRAP. Motivated by the SPT measurements, they suggested 
a gaussian subdiffusion model for the analysis of the FRAP 
curves. The given data sets were fitted equally well by the 
subdiffusion model and the diffusion model with incomplete 
recovery, a fit to a combined model yielded a recovery of 90% 
and a subdiffusion exponent a ^ 0.6. It was shown how sub- 
diffusion could be discriminated from an immobile fraction if 
the recovery curves would span a much larger time window. 

In contrast, almost complete recovery of fluorescence was 
found in FRAP experiments on the plasma membrane of COS- 
7 cells for a variety of putative raft and non-raft proteins [296]. 
The recovery curves from bleaching a large, 4 |jm wide stripe 
were rationalised with normal diffusion, yielding diffusion 
constants between 0.1 and 1.2 |jm^/s for the different proteins. 
These values appear to be mainly determined by the type of 
anchorage, since they change uniformly upon perturbing the 
membrane by variation of the cholesterol level or temperature. 

In another FRAP study on four different membrane pro- 
teins expressed in COS-7 cells [297], it was found that trans- 
port over micron-scale distances is normal with 80-90% fluo- 
rescence recovery. The diffusion constants, obtained from fit- 
ting a single exponential to the recovery data, are in agreement 
with other in vivo measurements, but are an order of magni- 
tude smaller than those of GPI-linked proteins reconstituted 
into liposomes [297]. The diffusion constants of the mem- 
brane proteins were unaffected by manipulations of the cor- 



tical cytoskeleton, but they increased markedly (up to 3-fold) 
if the protein density in the membrane was reduced in semi- 
intact cells. In conclusion, the reduced mobility of membrane 
proteins is likely a consequence of the crowdedness of cellu- 
lar membranes leading to anomalous transport on smaller time 
and length scales than accessible by a typical FRAP setup. 

Baker et al. [298] employed FRAP experiments with a 
variable radius of the observation spot (vrFRAP) to address 
the distribution and transport of CD4 proteins and chemokine 
receptors on the plasma membrane of living HEK 293T cells. 
The obtained ratio of incomplete fluorescence recovery and the 
apparent diffusion constants depends sensitively on the radius, 
which was varied between 1.40 and 3.45 pm. Interpretation 
of this dependency allowed for the conclusion that proteins 
and receptors are compartmentalised into domains of about 
350 nm. 

5.2.2. Model membranes fin vitro) 

Model membranes reduce the complexity of cellular mem- 
branes and permit a detailed study of macromolecular trans- 
port under controlled lipid composition and crowding condi- 
tions, vital for a physical understanding of anomalous trans- 
port. Membranes are usually modelled by bilayers of phos- 
pholipids, either supported by a flat substrate or forming a gi- 
ant unilamellar vesicle. The structure of the lipid bilayer can 
be modified by addition of another lipid type, and crowding ef- 
fects can be mimicked by covering the bilayer with anchored 
proteins or by introducing protein inclusions. It is possible to 
fluorescently label lipids as well as proteins and to follow their 
motion individually. 

Schiitz et al. [169] performed SPT experiments on a sup- 
ported POPC lipid bilayer modelling a fluid membrane, where 
highly diluted fluorescent TRITC-DHPE lipids were intro- 
duced in the outer leaf It was found that 69% of the labelled 
lipids showed normal, free diffusion with D - 4.4 while 
the motion of the remaining part was still diffusive, although 
suppressed by a factor of 63. These results rely on the analysis 
of the cumulative probability distribution of the square dis- 
placements for a series of large times (from 20 ms to 300 ms), 
an approach that is well suited to identify and account for a 
mixture of fast and slowly diffusing molecules. The data for 
the slow fraction indicate further that motion is confined at 
short times below 100 ms, which the authors interpreted as ev- 
idence for corrals of about 130nm size after correcting for the 
tracking uncertainty. 

A systematic development of anomalous diffusion was ob- 
served in a supported bilayer using SOPC lipids where the mo- 
tion of the monitored TRITC-DHPE lipids was hindered by 
essentially immobile "obstacles" [170], the obstacles were re- 
alised by lipopolymers in the inner leaflet tethered to the sub- 
strate. Similarly as in Ref. 169, the distribution of squared 
displacements was analysed for a fixed lag time of 50 ms. 
At obstacle concentrations below 10mol%, the distribution 
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model membrane 
(composition) 


probe (size) 


experimental 
technique 


temporal and 
spatial scales 


observation 


year 


reference 


SLB (POPC) 


DHPE lipid 


SPT 


20-300 ms 


normal, two components: 
£>! = 4.4nm2/s, O2 = Oi/63 


1997 


ref. 169 


SLB (SOPC) + teth- 
ered lipopolymer 
("obstacles") 


DHPE lipid, 
bacteriorhodopsin 


SPT 


flag = 50 ms 


anomalous diffusion, 

D I linearly for [lipopolymer] f 


2005 


ref. 170 


SLB 

(DOPC/cholesterol) 


DHPE lipid 


PCS 


beam radius: 
0.25 [im 


D: 4.2 -» 0.5 [jmVs for 
[cholesterol] f, deviations from 
simple diffusion 


2003 


ref. 305 


multilamellar 
vesicles 


DPPE lipid 


PCS 


T,/,: 108-230 ms 


two-component diffusion for 
fluid-gel coexistence in bilayer, 

IICLCIU^CIICUUS SIIUCIUICS 


2005 


ref. 306 




C5-Bodipy-PC 


variable- 
waist PCS 


T|/9: 1 - 150ms 
w: 210-300 nm 


w-dependent diffusion constant: 

Tl/2(w) = fo + wV4Deff, 

fo < in gel phase 


2011 


ref. 307 


multibilayers 
(DMPC/cholesterol) 


phospholipid 


FRAP 


spot radius: 4 ptm 


normal diffusion, D 2-fold 
reduced for [cholesterol] f 


1992 


ref. 308 


GUV 

(DOPC/DOPG) 
+ proteins 


integral mem- 
brane proteins 
(0.5-4 nm) 


PCS 


Ti/2 a 2 ms 


D i for [protein] t, max. 2-fold 
suppression, subdiffusion: 
a > 0.88 


2009 


ref. 309 


SLB (SOPC) 
+ avidin 


DPPE-anchored 
avidin protein 


PCS 


flag: 0.3-300 ms, 
w w 0.17 ^im 


development of subdiffusion, 
a: 1 ^ 0.68 for [avidin] T, 
spatial heterogeneity, no ageing 


2010 


ref. 191 


SLB (SOPC) 
+ neutravidin 


lipid, DOPE- 
anchored neutra- 
vidin protein 


PRAP 


Ti/2 « 3 S 


normal diffusion, Dprotein up to 
7-fold suppressed, Dupn almost 
unchanged as [neutravidin] t 


2009 


ref. 310 



Table VL Overview of in vitro experiments on model membranes as guide to the discussion in Section 5.2.2. Empty fields repeat the entry 
above, up/down arrows indicate increase/decrease of some quantity. Abbreviations: SLB. . . supported lipid bilayer, GUV. . . giant unilamellar 
vesicle, [X]. . . concentration of X. 



histograms at this time scale coiTespond to normal diffusion 
and become increasingly anomalous for larger concentrations. 
The apparent diffusion constant (aR^^ /4fiag is linearly sup- 
pressed with increasing obstacle concentration and vanishes 
near 40mol%. The scenario is consistently described by ran- 
dom walks on two-dimensional lattice percolation clusters, see 
Refs. [50, 107] and Section 3.3. Similar results were obtained 
for the transport of diluted bacteriorhodopsin inclusions mim- 
icking membrane proteins. The work demonstrates the great 
potential of model systems systematically bridging between 
physics and biology. 

The influence of cholesterol content on the transport in pla- 
nar DOPC bilayers supported by an atomically flat mica sub- 
strate was studied by FCS of rhodamine-labelled DHPE trac- 
ers [305]. The diffusion constant of DHPE reduces from 4.2 
to 0.5 pm^/s as up to 60% of DOPC is replaced by cholesterol. 
Deviations from normal diffusion arise for a cholesterol con- 
tent of 30% and more. While the data are fitted by the model 
with two normally diffusing components, the diffusion con- 
stant of the fast component is much larger than without choles- 
terol. Thus, the two-component model may not be suitable for 
DOPC/cholesterol bilayers. 

Aiming at the reduction of coupling effects with the sub- 



strate, Hac et al. [306] studied multilamellar model mem- 
branes consisting of stacks of about 50 layers. The bilayers are 
composed of a mixture of DMPC and DSPC lipids, which ex- 
hibits phase separation and a fluid-gel transition in each com- 
ponent. The motion of fluorescent TRITC-DPPE tracer lipids 
was followed by FCS at different compositions and tempera- 
tures covering the full region of fluid-gel coexistence in the 
phase diagram. While transport is normal in the fluid regimes, 
complex transport emerges in the phase-separated regime with 
its heterogeneous structure. The FCS correlations can be fitted 
by a two-component model, but again, the obtained diffusion 
times deviate from the ones in the pure fluid and gel phases. In- 
stead of relying on a phenomenological fit model, the authors 
performed extensive Monte-Carlo simulations which quantita- 
tively describe the experimental FCS curves; a detailed discus- 
sion of the simulations will be given in Section 5.2.3. 

In a recent FCS experiment using the same DMPC/DSPC 
mixture [307], spatio-temporal information on the lipid motion 
was collected by variation of the detection area. FCS correla- 
tion functions were fitted by the model for normal diffusion 
yielding the half-value times as function of the area, which 
were analysed in terms of the so-called FCS diffusion law, i.e., 
Ti/iiw) = fo + w^/4Z)eff. In the fluid phase of the mixture, the 
dependence on is linear (fo = 0) as expected for simple dif- 



39 



fusion. Extrapolating to zero area, a negative offset, to < 0, de- 
velops as the fluid-gel coexistence is traversed by lowering the 
temperature. Such a negative offset typically arises for motion 
in a meshwork of permeable barriers [194], but there, transport 
is not subdiffusive in the time domain. For motion in a phase- 
separating lipid mixture, Monte-Carlo simulations show that 
the area 4Defffo correlates with the typical domain size of the 
spatially heterogeneous environment, demonstrating that a do- 
main size may be extracted from FCS experiments [307]. The 
connection between the half-value times of the FCS correla- 
tion and the detection area may alternatively be described by 
a power law (see supplement of Ref. 307), ti/2(w) ~ w^^", 
yielding exponents a between 0.5 and 1; the smallest values 
are realised for temperatures close to the phase boundaries. 

The lateral diffusion of a phospholipid probe was mea- 
sured with FRAP in multi-bilayers of DMPC/cholesterol mix- 
tures, exploring a large part of the phase diagram in the 
composition-temperature plane [308]. In all cases, the fluo- 
rescence recovery was complete and the uniform disc model 
for free diffusion, Eq. (1 14), fitted the data. A 2-fold reduction 
of the diffusion constant was observed as cholesterol content 
was increased from 0% to 50%, and the diffusivity increased 
up to 8-fold as the temperature was changed from 24 °C to 
58 °C. 

The lateral mobility of fluorescently labelled integral 
membrane proteins reconstituted in giant unilamellar vesi- 
cles (GUVs) of a DOPC/DOPG mixture was measured us- 
ing FCS [309]. The hydrodynamic radius of the investigated 
proteins ranged between 0.5 and 4nm. Consistent with the 
Saffman-Delbruck model, the diffusion coefficient at low pro- 
tein content displays a weak, i.e., logarithmic, dependence on 
the radius. The diffusion coefficient of proteins and lipids de- 
creases linearly with increasing the protein concentration in 
the membrane up to 3000 proteins per pm^. Such a concentra- 
tion is an order of magnitude smaller then in biological mem- 
branes and induces only small crowding effects. Hence, the 
observed reduction of mobility did not exceed about 20% for 
the lipids and was at most 2-fold for the proteins; deviations 
from simple diffusion were only moderate with the exponent 
of subdiffusion a not falling below 0.88. 

Crowding by peripheral membrane proteins was addressed 
by a study of avidin proteins irreversibly bound to biotiny- 
lated DPPE lipid anchors in a supported single SOPC lipid 
bilayer [191, 311]. The content of lipid anchors in the bi- 
layer dkectly controls its protein coverage. X-ray reflectivity 
experiments showed that the protein layer is separated from 
the lipids by a distinct water layer, and continuous bleach- 
ing experiments indicated that the fluidity of the lipid bilayer 
is retained even at high protein coverage [311]. Upon in- 
creasing the protein coverage, FCS measurements on dilute 
labelled avidin proteins revealed a systematic development of 
anomalous transport [191]. Converting the FCS correlation to 
the mean-square displacement via Eq. (91b), subdiffusive mo- 



tion was observed over a time window spanning from 0.3 to 
300 ms. The exponent a decreases gradually from 1 to 0.68 
in the most crowded regime with an excess of biotin anchors. 
The onset of anomalous transport occurs already at an anchor 
concentration of 0.1 mol% or an area coverage of 3-5%. This 
transition regime is characterised by a pronounced and long- 
lived spatial heterogeneity: different exponents a were repro- 
ducibly obtained at different spots of the sample. The distribu- 
tion of exponents becomes bimodal as crowding is increased; 
some spots exhibit simple diffusion of the proteins, while an 
increasing fraction displays subdiffusion. On the time scale of 
hours, ageing effects could not be detected. 

FRAP measurements on similar SOPC bilayers crowded 
with neutravidin [310] corroborate normal diffusion of the 
lipids with slightly reduced diffusion constants by up to 25% 
at the highest protein coverage, a fraction of about 10% of the 
lipids was detected to be immobile. Covering the bilayer with 
fluorescently labelled neutravidin, the fluorescent signal from 
the proteins did not recover on the time scales of the experi- 
ment, which was attributed to proteins being arrested in a kind 
of gel phase. Measuring protein transport by means of con- 
tinuous bleaching, a 7-fold slowing down was observed at the 
highest protein coverage, which may be compatible with the 
findings in Ref. 191 taking the different time and length scales 
of the experiments into consideration. 

5.2.3. Computer simulations (in silicoj 
Obstructed motion 

In his seminal works, Saxton [87, 312] introduced a min- 
imalist model for transport in cellular membranes. The posi- 
tions of lipids and proteins are restricted to a two-dimensional 
lattice for computational efficiency; tracer molecules perform 
a random walk which is hindered by randomly distributed ob- 
stacles occupying a single lattice site each. In the limiting case 
of immobile obstacles, the obtained diffusion constants are 
almost linearly suppressed with increasing obstacle area and 
vanish at the percolation threshold [312]. While the sharp lo- 
calisation transition disappears if the obstacles themselves are 
allowed to diffuse slowly, the systematic reduction of tracer 
mobility with increasing obstacle concentration is conserved. 
For immobile obstacles [87], the MSD displays subdiffusive 
motion over a growing time window and a crossover to nor- 
mal diffusion at large time scales. The apparent exponents a 
obtained from the intermediate subdiffusive regime depend 
on the obstacle concentration and approach their universal 
value only at the percolation threshold. The concentration- 
dependence of the exponent was reconciled with the dynamic 
scaling hypothesis only recently by including the leading cor- 
rections to scaling, which are again universal [116]; see Sec- 
tion 3.3.4 for a detailed discussion. 

The spatial discreteness of lattice models is overcome by 
continuum models, where the obstacles are usually imple- 
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space crowded medium tracer interactions observation year reference 



square 
lattice 


point obstacles 
(mobile, but slow) 


point 


excluded volume 


D i for [obstacle] T 


1987 


ref. 312 




point obstacles 
(fixed, random) 






subdiffusion on growing time window, 
D i linearly and a J, for [obstacle] T, at 
percolation transition: £) = 0, or ~ 0.68 

dynamic scaling for cluster-resolved 
MSDs, universality of exponents 


1994 
2008 


ref. 87 
ref. 116 


2D con- 
tinuum 


random disc obstacles 


discs 


hard core repulsion 


transient subdiffusion, a [ for 

[obstacle] f , highly fragmented void space 


2006, 
2008 


refs. 83, 
84 




random disc obstacles 
(overlapping) 


point 


hard core repulsion 


universal scaling of localisation 
phenomenology, subdiffusion over 
6 decades in time, cy: 1 — > 0.659 — » 


2010 


ref. 22 


square 
lattice 


extended obstacles from 
AFM images 


point 


excluded volume 


subdiffusion, a: 1 — > 0.59 as excluded 
area t, agreement with experiments 


2011 


ref. 313 




mobile "rafts", point and 
line obstacles ("fences") 


point 


excluded volume, low 
mobility inside rafts 


subdiffusion (a « 0.7) only for immobile 
rafts and raft-excluded tracers 


2007 


ref. 314 




regular network of line 
obstacles ("fences") 


point 


potential barrier 


high and dense barriers reduce Dit) by 
several orders of magnitude 


2008 


ref. 315 


triangular 
lattice 


DMPC/DSPC monolayer 
model (2 sites per lipid) 


lipid 


nearest neighbour, 
slow and fast hopping, 
fluid/gel state flips 


FCS correlation: subdiffusion [a > 0.63), 
agreement with experiment 

subdiffusion {a: 1-0.6) from first passage 
times, structural heterogeneity 

FCS with variable beam waist: 
Ti/2(w) ~ w^/", a: 1-0.5 


2005 

2005 
2011 


ref. 306 

ref. 316 
ref. 307 


square 
lattice 


DMPC/DSPC monolayer 
model (1 site per lipid) 


point 




transient subdiffusion, a » 0.85 close to 
critical point 


2011 


ref. 317 








-1- heterogeneous 
sticky network 


transient subdiffusion {a x 0.5), 
DolD up to 10 


2011 


ref. 318 


3D con- 
tinuum 


atomistic DOPC bilayer 
(288 fully hydrated lipids) 


lipid 


molecular force fields 


lipid transport at 1 ps is correlated over 
distances of 2.5 nm 


2009 


ref. 319 



Table VII. Overview of computer simulations on crowded membranes as guide to the discussion in Section 5.2.3. Empty fields repeat the 
entry above, up/down arrows indicate increase/decrease of some quantity. 



mented by a simple geometric shape. Sung and Yethiraj 
[83, 84] investigated the motion of tracer discs in a frozen 
obstacle matrix of randomly placed, non-overlapping, impen- 
etrable discs; the limit of dilute tracers is known as non- 
overlapping Lorentz model in the literature [320]. Similarly 
as for the lattice model, transport slows down with increas- 
ing obstacle concentration and subdiffusive motion emerges in 
an intermediate time window. A Voronoi map of the continu- 
ous space accessible to the tracer particles reveals that the free 
space is highly fragmented into voids of different sizes, which 
are connected by significantly correlated channels. Thus at 
high obstacle concentration, the primary mode of transport re- 
sembles a hopping motion between neighbouring voids; sim- 
ilar observations were reported from particle tracking experi- 
ments [285]. Transport on large scales, however, is the result 
of a large number of hopping events, and subdiffusive motion 
relies on the presence of a spacious and self-similar network 
of connected channels meandering the membrane. 

The full dynamic scaling picture of anomalous transport 



due to a localisation transition was established only recently 
for two-dimensional continuum models [22]. Extensive simu- 
lations for a Brownian tracer between uncorrelated and over- 
lapping circular obstacles confirmed the universal values of 
the dynamic exponents at leading and next-to-leading order, 
e.g., for the exponent of subdiffusion a -2jz~ 0.659 if also 
tracers in isolated regions ("finite clusters") are included. The 
subdiffusive motion was followed over six decades in time, 
and the crossover to normal diffusion away from the local- 
isation transition was found to extend over several decades 
as well. Further, the transport dynamics was rationalised in 
terms of crossover time and length scales, which diverge at the 
transition. Appropriate rescaling of the time-dependent dif- 
fusion constants for different obstacle concentrations yielded 
data collapse onto a universal scaling function that describes 
the crossover from subdiffusion to normal diffusion. The tran- 
sition to continuum models is a step towards a more realistic 
modelling, the universal phenomenology of the localisation 
transition, however, is preserved as expected at sufficiently 
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large scales. This is in contrast to obstructed transport in three 
dimensions, where lattice and continuum models belong to dif- 
ferent (dynamic) universality classes with different transport 
exponents [114, 115, 117]. 

An uncorrected distribution of obstacles, generated by a 
random process, is certainly a simplification. Realistic config- 
urations of the excluded area were obtained from atomic force 
microscopy (AFM) imaging of phase-separating lipid bilay- 
ers [313]. For mixtures of DSPC/DOPC lipids, the gel-like 
DSPC-rich phase effectively provides an excluded area to the 
motion of DMPE tracer lipids. The rastered AFM images then 
yield "obstacles" of irregular shape that extent over many lat- 
tice sites, and Monte-Carlo simulations can be used to inves- 
tigate the transport of the tracer lipids performing a random 
walk in the remaining space. This approach is somewhere in 
between the original lattice models and the continuum mod- 
els, and if it is justified that the macroscopic transport is not 
affected by the imaging resolution (cf. Ref 315) it can provide 
a quantitative bridge to the experiments. 

Apart from excluded volume, transport of membrane pro- 
teins may be hindered by other sources. Nicolau et al. [314] 
studied random walks of proteins on two-dimensional lattice 
models, where the proteins are hindered by either a regular lat- 
tice of cytoskeletal fence posts or by mobile lipid rafts that are 
either impenetrable for proteins or reduce their mobility. Ap- 
preciable subdiffusion similarly as for randomly distributed, 
fixed obstacles was only observed if the protein motion is hin- 
dered by almost immobile, impenetrable rafts. Depending on 
the raft mobility, subdiffusive motion was characterised by ex- 
ponents a between 0.65 and 0.75. The simulations indicated 
further that collisions with fence lines alone (without binding) 
do not explain the experimentally reported degrees of anoma- 
lous transport. Fences that impose very high barriers, however, 
can reduce protein diffusion at long times by several orders 
of magnitude [315], which was demonstrated by Monte-Carlo 
simulations on a lattice compartmentalised by a regular net- 
work of fences. It appears that here a plateau in the MSDs 
rather than a subdiffusive power-law growth is the manifesta- 
tion of slow transport. Since only a single length scale char- 
acterises the regular super-lattice of fences, the plateau can be 
understood as a reflection of the mesh size. 

Mixtures of DMPC/DSPC lipids 

Experiments on lipid transport in artificial lipid bilayers 
have shown that anomalous transport emerges as function of 
the lipid composition. Often binary mixtures of lipids like 
DMPC/DSPC or DOPC/cholesterol were used, which can 
phase separate into fluid- and gel-like domains. Sugar et al. 
[321] developed a model for DMPC/DSPC mixtures on a tri- 
angular lattice, where each of the two acyl chains of a lipid oc- 
cupy two adjacent sites and each site is in either the gel or the 
fluid state. The thermodynamics of the model is governed by 
a Hamiltonian with 10 relevant parameters, which were either 



directly inferred from calorimetric measurements or estimated 
by means of Monte Carlo simulations to reproduce the correct 
phase diagram. 

Relying on this model, Hac et al. [306] investigated FCS 
correlation functions of the lipids using Monte-Carlo simu- 
lations with a coarse-grained dynamics. Each lattice update 
allowed for the possibility of flipping the gel/fluid state of a 
lipid chain (Glauber dynamics) and of exchanging adjacent 
molecules mimicking diffusion (Kawasaki dynamics). The 
dynamics is controlled by a rate for the state flips and two 
hopping rates for fluid- and gel-like environment; the latter 
two were matched to the experimental diffusion constants for 
pure fluid or gel phases. The numerically obtained FCS curves 
show quantitative agreement with the experimental results if 
the flip rate is chosen close to the hopping rate in the fluid. 
In the regime of gel-fluid coexistence, the FCS curves deviate 
from normal diffusion and are fitted by the gaussian subdif- 
fusion model. The exponent a depends on the flip rate, and 
its smallest value of 0.63 was reported for pure Kawasaki dy- 
namics, i.e., in the absence of artificial state flips. Parameters 
for the hopping rates may be obtained by free volume theory, 
and correlating the structure and distribution of gel clusters 
with the lipid dynamics provides evidence that the anomalous 
transport arises from molecules exploring a heterogeneous en- 
vironment [316]. The model was further employed to address 
spatial aspects of the lipid transport by simulating FCS cor- 
relations for a wide range of beam waists, which allowed for 
the extraction of the mean domain size [307]; a discussion was 
already given in Section 5.2.2. 

The above lattice model for a lipid membrane of Refs. 
321 and 306 was further simplified by Ehrig et al. [317, 318]. 
They used a square lattice with each site representing a lipid 
molecule; the coarse-grained dynamics still comprises state 
flips of a single lipid and diffusion by exchange of two ad- 
jacent sites. Together with advances in computer hardware, 
this permitted the investigation of comparably large systems 
with a membrane area of (0.48 pm)^ or 600^ lattice points and 
processes on time scales of up to 1 s. The simulations demon- 
strate that close to the upper critical point of fluid-gel coexis- 
tence, the lipid motion exhibits anomalous transport over sev- 
eral decades in time. A double-crossover scenario is observed 
from simple diffusion at microscopically short time scales to 
transient subdiffusion and finally back to normal diffusion. 
The time window of subdiffusion increases as the critical point 
is approached, the local exponent ait) of the MSD, however, 
does not yet display a plateau and its smallest reported value 
is around 0.85 being still relatively close to 1. In consequence, 
only a 3-fold suppression of the macroscopic diffusion con- 
stant is observed compared to its microscopic value. The in- 
troduction of a frozen, spatial heterogeneity in form of a sticky 
network [317] renders the lipid transport significantly more 
anomalous in a limited time window and induces an up to 10- 
fold reduction of mobility and minimal values of the local ex- 
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ponent a(t) down to 0.5. 

Let us raise the question whether an asymptotic, univer- 
sal value for the subdiffusion exponent exists in a mixture 
as the critical point is approached further The analytical 
work by Kawasaki [322], employing a local-equilibrium ap- 
proximation, suggests that the self-diffusion constant is only 
mildly affected by the critical fluctuations; in particular, it is 
not singular at leading order and does not vanish. Similarly, 
recent microcanonical molecular dynamics simulations for a 
three-dimensional symmetric binary mixture near its conso- 
lute point revealed almost negligible anomalies in the single- 
particle transport [323, 324]. It appears that in two dimen- 
sions, the single-particle transport is coupled more strongly to 
the critical fluctuations. But it is unlikely that subdiffusion can 
exist over an arbitrarily large time window, which would imply 
a vanishing diffusion constant. 

With constantly growing computational resources, atom- 
istic molecular dynamics simulations may come into reach 
as an alternative source for the parameters of coarse-grained 
models. A detailed view of the structure of a DMPC/DSPC bi- 
layer was provided by a hybrid simulation technique, where an 
atomistic molecular dynamics simulation was supplemented 
by a semi-grandcanonical Monte-Carlo step converting a 
DMPC lipid into DSPC and vice versa [325]. Using 128 lipid 
molecules and 3200 explicit point charge water molecules, it 
was possible to generate system trajectories covering 50 ns in 
time, or two million integration steps. By means of a pure 
molecular dynamics simulation for a bilayer of 288 atomisti- 
cally modelled DOPC lipids, the experimental values of the 
short-time diffusion constants of the lipids could be quantita- 
tively reproduced by a simulation [319]. The study further 
showed that the motion of the lipids is correlated over dis- 
tances of about 2.5 nm at the time scale of 1 ps, which is a use- 
ful piece of information for the refinement of lattice models: 
it suggests that a single lattice site should represent a corre- 
lated area of several lipids rather than a single lipid molecule 
or chain. Molecular dynamics simulations of lipid transport 
over large distances seem not yet feasible as it would be ham- 
pered by the small diffusion constant of the lipids in the gel 
state (Z)gei ~ 0.05nm-/[is [306]). This emphasises the need of 
controlled coarse-graining and multiscale approaches for the 
study of slow, complex relaxation dynamics at the micro- and 
millisecond scale. 

5.3. Reaction kinetics 

So far, we have discussed the anomalous transport of individ- 
ual macromolecules in their crowded cellular environments. 
Physiological processes, however, are determined to a large 
extent by biochemical reactions. Adding to the complexity of 
anomalous transport, the reaction rates as well as the availabil- 
ity of reactants are significantly modified by molecular crowd- 
ing. The systematic investigation of this topic is relatively 
young, and we wiU only sketch some of the recent findings. 



Minton et al. [227, 326] considered the consequences of 
macromolecular crowding on the equilibrium rates of bio- 
chemical reactions from a thermodynamics point of view, fo- 
cusing on changes of the excluded volume. They found a non- 
specific enhancement of reaction rates if the total excluded vol- 
ume is reduced; examples for such reactions are the formation 
of macromolecular complexes in solution, binding of macro- 
molecules to surface sites, formation of insoluble aggregates, 
and compaction or folding of proteins. For protein association 
reactions, crowding is generally expected to increase the rate 
of slow, transition-state-limited reactions and to decrease the 
rate of fast, diffusion-limited ones. In the context of polymer 
systems, it was predicted that the subdiffusion exponent con- 
trols the crossover between reaction- and diffusion-limited ki- 
netic regimes [327]. A variation of the association rates under 
crowding conditions by an order of magnitude was reported in 
the experimental literature reviewed until 2008 in Refs. 227 
and 326. Melo and Martins [328] reviewed the kinetics of 
bimolecular reactions in model bilayers and biological mem- 
branes from an experimental perspective until 2006; their con- 
clusion "When going to crowded or phase-separated systems 
. . . the absence of experimental studies is almost complete." 
emphasises the need for systematic investigation of the reac- 
tion kinetics in crowded membranes. 

The kinetics of the reaction A H- B — > products on a perco- 
lating, two-dimensional lattice was studied using Monte-Carlo 
simulations [329]. Subdiffusive motion slows down the initial 
encounter of reactants, but concomitantly hinders them from 
escaping due to the increased return probability. The simu- 
lations showed that obstructions can slightly increase the re- 
action rate if the reaction probability on a single encounter is 
small. On the other hand, the reaction rate is reduced for high 
reaction probability, and taken together, subdiffusion of the re- 
actants narrows the range of reaction rates. The bimolecular 
reaction is similar to that of the reaction of a random walker 
with an immobile target. Realising anomalous transport of 
the walker by a hierarchy of non-reactive binding sites. Sax- 
ton [330] found from simulations that trapping may contribute 
significantly to the noise in reaction rates and further that the 
subdiffusive motion yields a power-law distribution of capture 
times. Schmit et al. [331] identified the concentration of reac- 
tants and the crowdedness of the medium as two antagonistic 
parameters for the rates of diffusion-limited bimolecular reac- 
tions. They showed that there is an optimal reaction rate (for a 
fixed ratio of the concentrations of reacting and inert particles) 
if the reactants are not too diluted and if the medium is not too 
crowded, i.e., diffusion of the reactants is not too slow. Con- 
tinuum models for the reaction of ligands with an immobile 
receptor realised the crowded environment by a dense hard- 
sphere fluid [332, 333]. Again, a non-monotonic dependence 
of the effective reaction rate on the volume fraction of the fluid 
particles was found for sufficiently small receptors. 

For a diffusion-hmited dimerisation reaction in a computer 
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model for the E. coli cytosol [334], it was observed that the 
reaction coefficients become time-dependent in the presence 
of anomalous transport; the obtained reaction rate was sup- 
pressed by two orders of magnitude and obeyed an approxi- 
mate power law over 4 decades in time. Hellmann et al. [335] 
found such a fractal reaction kinetics with a time-dependent, 
power-law reaction rate also in two-dimensional, off-lattice 
simulations for the reaction A + B — > with subdiffusive mo- 
tion of the reactants. Specifically, fractional Brownian motion 
was implemented by means of an overdamped Langevin equa- 
tion with a power-law-correlated noise. Key findings were 
that i) independently of the reaction probability the diffusion- 
limited formation of new reactant pairs dominantly determines 
the reaction rate and ii) subdiffusive motion enhances the seg- 
regation of reactants. Extending this approach to describe a 
multi-step enzymatic reaction, the double-phosphorylation of 
mitogen-activated protein kinase (MAPK), an increased per- 
formance due to anomalous diffusion was found when disso- 
ciation rates of the intermediate enzyme-substrate complexes 
are high while the irreversible catalytic step is slow [336]. 

Froemberg and Sokolov [337] set up a continuum descrip- 
tion in form of a reaction-iM^jdiffusion equation for the reac- 
tion A H- B — » with subdiffusive reactant motion governed by 
a CTRW. Supplying reactants of both types from both sides of 
the medium, stationary profiles of reactant concentrations and 
of reaction zones develop. These profiles exhibit accumulation 
and depletion zones close to the boundaries, and the reaction 
intensity was found to depend non-mono tonic ally on the reac- 
tant concentration at the boundaries. A different approach was 
taken by Boon et al. [338] for the n*-order annihilation reac- 
tion, A-I-A-I-----I-A— >0, who derived a nonlinear reaction- 
diffusion equation assuming concentration-dependent diffu- 
sion and reaction coefficients and imposing a scaling form 
of the solution, i.e., demanding the concentration to satisfy 
c(r, t) - r"I^C{rr"^^). Thereby, subdiffusive motion of the 
reactant is incorporated, cf. Eqs. (3), (36), (71) and (75). Then, 
the resulting steady-state profiles exhibit a long-range, power- 
law decay and fit experimental data for the Wg morphogen 
gradient in Drosophila wing disc. 

6. Conclusion 

The advancements of biophysical experimental tools in the 
past two decades have revolutionised the investigation of in- 
tracellular transport of proteins, nucleic acids, and artificial 
tracers as well as the motion of membrane receptors and lipids. 
The new methods allowed for a detailed view of the plethora of 
dynamic processes in vivo directly at the cellular level. Numer- 
ous new phenomena and features of the transport have been 
discovered as more and more data of ever increasing quality 
have been collected. Living cells, by nature, consist of many 
differently sized and shaped constituents, each with their own 
physiological function, leading to densely packed, highly het- 



erogeneous structures at the nanoscale — commonly referred 
to as macromolecular crowding. A widespread observation 
is that transport in living cells behaves rather different from 
the standard picture of Brownian motion, where the erratic 
trajectories are described by a probabilistic law as dictated 
by the central-limit theorem. Instead, the mean-square dis- 
placement, being the simplest quantitative measure of single- 
particle transport, displays complex and often also subdiffu- 
sive behaviour. Here, we have focused on anomalous trans- 
port, characterised by power-law increases on time windows 
covering several orders of magnitude. A unified physical pic- 
ture of this spectacular anomaly remains one of the grand chal- 
lenges of biophysics. 

From a theoretical perspective, a fundamental issue is to 
identify mechanisms that lead to a violation of the central-limit 
theorem on macroscopic time scales. To this end, several the- 
oretical models have been established that achieved this goal 
and yield a subdiffusive increase of the mean-square displace- 
ment. Here, we have summarised three of the most commonly 
used approaches based on different microscopic processes. 

Fractional Brownian motion as a representative of spatially 
gaussian models is distinguished by a strictly self-similar in- 
crease of the mean-square displacement. Its most prominent 
feature is a power-law correlated, gaussian random noise. As 
a consequence, the stochastic process is stationary, but not 
Markovian. The transport propagator, i.e., the probability dis- 
tribution of the displacements after a given lag time, displays 
a scaling property similar to simple Brownian motion. 

Second, continuous-time random walk models describe a 
microscopic hopping process, which asymptotically generates 
subdiffusive motion in the case of waiting-time distributions 
with sufficiently broad power-law tails. It contains a fractional 
Fokker-Planck description as limiting case where the propa- 
gator exhibits scaling for all time and length scales. The scal- 
ing function depends on the subdiffusion exponent a and is 
markedly different from the gaussian case. In particular, the 
non-gaussian parameter does not vanish even for long times 
and approaches a characteristic value depending only on a. 
By the very construction, the process is non-stationary if the 
mean-waiting time is literally infinite, which implies ageing 
phenomena and weak ergodicity breaking. 

The third class consists of Lorentz models, focusing on ob- 
structed motion in strongly heterogeneous landscapes. Crowd- 
ing agents lead to ramified structures delimiting the space 
accessible to the tracer, and anomalous transport emerges 
generically from the meandering motion of the tracer in such 
labyrinth-like structures. Subdiffusive motion occurs on larger 
and larger time windows as a localisation transition is ap- 
proached, until long-range transport would eventually cease. 
In the vicinity of the transition, scaling behaviour develops 
and can be rationalised in terms of a dynamic critical phe- 
nomenon. The properties upon approaching anomalous trans- 
port are controlled by three independent, universal exponents 
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governing the mean-square displacement, the crossover time 
scale, and the fractal spatial structure. The averaged dynam- 
ics is stationary, exhibits heterogeneous diffusion, and marked 
non-gaussian spatial transport. 

All three models are indistinguishable on the level of the 
mean-square displacement being subdiffusive. However, the 
predictions differ drastically if transport is considered on dif- 
ferent length scales, for example by analysing spatio-temporal 
information as encoded in the propagator. Then, the shape of 
the scaling functions serves, in principle, as a fingerprint to 
discriminate the various theories. Since the fuU scaling func- 
tion is often more difficult to measure than derived character- 
istics, we have exemplified the properties of the non-gaussian 
parameter as a potential clue in the quest for the microscopic 
origin of anomalous transport. 

For this endeavour, researchers have matured the experi- 
mental techniques that are suited to uncover the transport prop- 
erties of biophysical samples at mesoscopic scales. Single- 
particle tracking, fluorescence correlation spectroscopy, fluo- 
rescence recovery after photobleaching, among others, have 
widely been applied and yield complementary aspects of com- 
plex transport. Single-particle tracking provides access to the 
detailed particle trajectories and thereby to a complete statisti- 
cal characterisation, in particular to the full propagator. Imple- 
mentations based on video microscopy are typically restricted 
to Umited temporal windows at a resolution of the order of 
10 ms. Fluorescence correlation spectroscopy, on the other 
hand, probes shorter time scales and yields correlation func- 
tions over many decades in time. The shape of the correlation 
function serves as indicator for anomalous transport at the spa- 
tial scale of the illumination region. In recent developments, 
the beam waist has been made adjustable to probe the dynam- 
ics on a range of length scales, which then in principle allows 
for the discrimination of different propagators. Modern setups 
for fluorescence recovery after photobleaching are apt to de- 
tect slow transport on the time scale of several seconds and 
the length scale imposed by the bleaching spot; noteworthy, 
the method is sensitive to immobile particles as well. It moni- 
tors collective rather than single-particle transport, and anoma- 
lous transport becomes manifests in the recovery curve at long 
times. We see a need to evaluate the current theoretical ap- 
proaches for measurable quantities of each of the above tech- 
niques. On the other hand, the interpretation of experimen- 
tal data in the context of anomalous transport requires care- 
ful reconsideration of the underlying assumptions, e.g., spa- 
tially gaussian behaviour. Similarly, if the underlying process 
is (weakly) non-ergodic as is the case in CTRW or Lorentz 
models, time and ensemble averages generally do not coin- 
cide and caution is advised. Single molecule experiments are 
a promising tool to resolve these issues in the future. 

Experiments inside the cyto- and nucleoplasm of various 
eukaryotic cells and bacteria have provided ample evidence 
for subdiffusive motion in vivo over several decades in time. 



The overall investigated time scales covered the window from 
100 |js to 10 s, and the exponent a mostly ranged between 0.7 
and 0.85 depending on the specific experiment, but also on the 
size of the tracer; an even smaller exponent of 0.53 was re- 
ported for 5-nm gold beads in many different cell types. These 
findings are supported by in vitro measurements on crowded 
model solutions, yielding systematically decreasing subdif- 
fusion exponents a for increasingly crowded solution. The 
smallest values for a were obtained in dextran solution and 
reached down to 0.6 for dextran tracers and to 0.7 and 0.74 for 
globular proteins. A criticism on dextran as crowding agent 
is a possible interference of the tracer motion with the internal 
dynamics of polymeric dextran branches; for example, anoma- 
lous transport was regularly found to be less pronounced in 
solutions of globular proteins. 

Computer simulations have corroborated the experimental 
findings at least quaUtatively. SimpUfied model systems fo- 
cusing on the excluded volume describe a rich phenomenol- 
ogy of anomalous transport including the gradual emergence 
of subdiffusive motion, spatially non-gaussian transport, het- 
erogeneous diffusion with strongly suppressed diffusion con- 
stants, and a fraction of (on large scales) immobile particles. 
These phenomena have been related to a locaUsation transi- 
tion, and simulation data have been rationalised in terms of 
dynamic scahng. Subdiffusive motion over many decades in 
time appears already for a small tracer in a slowly rearranging 
matrix of equally sized particles with apparent exponents less 
than 0.5; it persists for infinitely long times for randomly dis- 
tributed obstacles at a critical value of the excluded volume, 
e.g., a universal value of the exponent a ^ 0.32 was reported 
from simulations following the subdiffusive transport over 6 
decades in time. More specific models for the cytoplasms of 
HeLa cells and E. coli use a distribution of particle sizes close 
to physiological conditions, yielding subdiffusive motion over 
several decades in time with exponents varying between 1 and 
0.55 as function of the tracer size. 

A number of experiments in vivo and in vitro, however, 
display simple, albeit slow diffusion with diffusion constants 
reduced by 1 to 3 orders of magnitude. For some experiments, 
data fitting required two normally diffusing components, ide- 
ally motivated by a physical picture hke compound forma- 
tion between tracer and crowding agent. For model solutions 
with adjustable degree of crowding, the reduction of diffusiv- 
ity is explained to some extent by an increased macroscopic 
viscosity. On the other hand, it has been observed that the 
Stokes-Einstein relation is violated for proteins as probe and 
that the visco-elastic response on the nanoscale strongly devi- 
ates from the macroscopic rheological behaviour. Eventually, 
transport was found to be less anomalous in simulations if hy- 
drodynamic interactions between the particles are taken into 
account; further, hydrodynamic interactions can already result 
in an up to 10-fold reduction of the short-time diffusion con- 
stant, in agreement with theory and recent scattering experi- 
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ments probing the nanosecond scale. 

Summarising, anomalous transport has ubiquitously been 
observed in experiments in vivo and can be recreated in 
crowded model solutions and in computer simulations, but not 
in all instances. The picture is not thus simple and the appear- 
ance of anomalous transport in experiments is far from being 
understood. At least, one can say for model systems in vitro 
and in silico that transport is slowed down upon increasing 
the concentration of crowding agent or increasing the size of 
the tracer. An unsolved problem is the importance of the dif- 
ferent contributions slowing down the motion, examples are 
excluded volume, polymer effects, compound formation, and 
hydrodynamic interactions. 

The situation seems even more controversial for the mo- 
tion of membrane proteins or phospholipids in cellular mem- 
branes. An central goal of such in vivo measurements is the 
elucidation of the membrane structure from the transport be- 
haviour by indirect evidence; such biophysical conclusions 
are beyond the scope of this work, and we kindly refer the 
reader to Refs. 19, 157, and 228 for pertinent reviews. In cel- 
lular membranes, essentially all kinds of transport have been 
reported ranging from normal diffusion with reduced diffu- 
sion constants over subdiffusion to confined motion and im- 
mobile tracers. Former measurements of membrane proteins 
displayed a heterogeneous dynamics in the sense that different 
tracers of the same type showed qualitatively different motion. 
During the past decade, advances in the experimental tech- 
niques have given access to the transport behaviour over large 
time windows. For example, high-speed tracking of phospho- 
lipids with 25 |js resolution probed the mean-square displace- 
ment over five decades in time and have provided evidence 
for a double-crossover scenario including a regime of subdif- 
fusive motion with exponent a » 0.53. Similar time scales 
can be covered by fluorescence correlation spectroscopy, but 
crossovers between different regimes have not been observed 
so far; the experimental data were typically compatible with 
either gaussian subdiffusion or normal diffusion of a fast and a 
slow component. For several membrane proteins, subdiffusion 
was reported at the millisecond scale with exponents typically 
between 0.5 and 0.8, while other experiments found normal, 
but very slow diffusion at the scale of seconds. Spatial aspects 
of transport have been addressed by FCS experiments with 
variable illumination area, revealing anomalous transport in 
the form of length-scale-dependent diffusion coefficients be- 
low about 100 nm. Thus, the diffusion coefficient may depend 
both on the time and length scale under investigation, and both 
should cover a sufficiently large window for a comprehensive 
description of macromolecular transport in membranes. 

As for cytoplasmic transport, the investigation of model 
systems is of paramount importance for a firm understanding 
of transport in membranes. In vitro studies of supported lipid 
bilayers or giant unilamellar vesicles agree that the lipid diffu- 
sion slows down in the presence of phase separation, which 



may be induced by increasing, e.g., the cholesterol content 
of the lipid bilayer It was found that gel-like regions effec- 
tively serve as a ramified, excluded area rendering the lipid 
motion anomalous. The experimental evidence has been cor- 
roborated by computer simulations for two-dimensional mod- 
els of lipid mixtures close to fluid-gel coexistence, yielding 
subdiffusive motion for molecules exploring a heterogeneous 
structure. Pronounced subdiffusion over three decades in time 
was detected for the motion of proteins anchored to a lipid bi- 
layer, with the exponent a developing from 1 to 0.68 as the 
bilayer was increasingly crowded by like proteins. Further, 
the transport of membrane proteins can become anomalous if 
a fraction of the membrane lipids is immobilised, e.g., by teth- 
ering to the substrate. This resembles obstacle models with a 
percolation threshold, where a ramified, heterogeneous land- 
scape of excluded area appears generically. As in three di- 
mensions, anomalous transport emerges in a finite temporal 
window upon approaching the percolation threshold, and sub- 
diffusive motion with a universal exponent a » 0.66 could be 
followed in simulations over 6 decades in time at the localisa- 
tion transition. 

As elaborated in this review, significant progress both on 
the experimental side as well as in theoretical modelling has 
been achieved in the puzzle of anomalous transport in crowded 
biological media. A lot of insight has been obtained, in par- 
ticular, to establish the phenomenon itself and to identify im- 
plications of crowding on macromolecular transport. Yet, the 
knowledge acquked also raises a series of new questions chal- 
lenging the frontier of biophysical research. Let us conclude 
this review by pointing at open problems and possible direc- 
tions of future work. 

So far, the interpretation of experimental data was mostly 
restricted to the mean-square displacement or the FCS corre- 
lation function assuming spatially gaussian transport. Only 
few studies discuss other statistical measures like distribu- 
tions of displacements or of apparent exponents or even two- 
particle correlations. We have emphasised that the charac- 
terisation of anomalous transport depends on the time and 
length scales probed, highlighting the importance of collect- 
ing spatio-temporal information. To this end, testable the- 
oretical predictions need to be elaborated and experimental 
methods to be perfected; a first step would be to focus on the 
non-gaussian behaviour, for which detailed predictions exist. 
Beyond that, multiple-time correlation functions contain in- 
formation on the stochastic process underlying the transport, 
which would quantify the non-Markovian behaviour in par- 
ticular. Here, a comparison between the different theoretical 
approaches remains to be done. 

While basic models have been established, they need to 
be refined to capture specific experimental conditions. For 
example, fractional Brownian motion should be generalised 
to account for crossover phenomena to normal diffusion, 
both at microscopically short and macroscopically long time 
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scales. Similarly, several experimental groups have excluded 
the continuous-time random walk model as a candidate for 
anomalous transport since ageing scenarios intrinsic to the 
model have not been observed. This apparent deficiency can 
be avoided by considering large, but finite mean-waiting times, 
thus allowing for a graceful exit to normal diffusion at long 
times. The Lorentz models dealing with obstructions should 
be extended to account for correlated, differently sized obsta- 
cles and slowly rearranging disordered environments. Further- 
more, the idealisation of a hard repulsive interaction should 
be relaxed, thus introducing barriers of finite rather than infi- 
nite height. Eventually, a significant ingredient to anomalous 
transport in the subcellular world are boundaries as imposed 
by (intra-)cellular membranes; the role of a finite space has 
been studied especially for the fractional Fokker-Planck equa- 
tion. 

Well-controlled experimental model systems offer a great 
potential to systematically investigate anomalous transport by 
examining adjustable crowding parameters over wide ranges. 
They provide the largely unexplored opportunity to address 
the same system by complementary techniques, allowing for 
a comprehensive description of the motion in crowded media. 
Further, they form a bridge from in vivo experiments to re- 
ductionist systems amenable to theory. Importantly, scaling 
behaviour could be tested here to discriminate between ap- 
parent subdiffusion as a mere phenomenological description 
of data in contrast to genuine anomalous transport over many 
orders of magnitude. As a central result of pushing experi- 
ments and theory further, it should become possible to predict 
the transport properties, at least qualitatively, from the outset 
of an experiment. On the other hand, it would be desirable 
to develop an experimental standard for subdiffusive trans- 
port, ideally covering sufficiently large windows in space and 
time, that is able to crosscalibrate the different experimental 
techniques [339]. Let us remark that the connection beyond 
the biophysical world, e.g., to transport phenomena in gels, 
porous media, and random heterogeneous materials, has so far 
remained largely elusive from an experimental point of view. 

In the broader context of cell biology, single-particle trans- 
port is only one aspect of macromolecular transport and one 
anticipates also fascinating collective phenomena associated 
with crowded environments. Similarly, living cells comprise 
an abundance of active, non-equilibrium processes, and their 
connection to anomalous transport needs to be unravelled in 
the future. The grand challenge, of course, is to clarify the 
physiological implications of anomalous transport, most im- 
portantly: does it constitutes merely a peculiarity or has it a 
biological benefit? For example, suggested scenarios on the 
acceleration of target search processes and bimolecular reac- 
tions need to be substantiated experimentally. More generally, 
important ingredients for systems biology Uke the kinetics of 
biochemical reactions, the dynamics of protein folding and un- 
folding, and intracellular signalling pathways may be strongly 



affected by macromolecular crowding. The development of 
a unified picture of anomalous transport and its physiological 
consequences will entail an immense effort of interdisciplinary 
research during the next decade at least. Or let us phrase it with 
Ellis and Minton [226]: "Join the crowd!" 
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